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Abstract 


An integral operator approach is used to derive solutions to static 
shape determination and control problems associated with large space 
structures. Problem assumptions Include a linear self-adjoint system 
model, observations and control forces at discrete points, and quadratic 
performance criteria for the comparison of estimates or control forces. 

Results are Illustrated by simulations. In the one dimensional case 
with a flexible beam model, and In the multidimensional case with a 
finite element mcdel of a large space antenna. 

Modal expansions for terms in the solution algorithms are presented , 
using modes from the static or associated dynamic model. These 
expansions provide approximate solutions in the event that a closed form 
analytical solution to the system boundary value problem is not available. 
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Static Shape Determination and Control for 
Large Space Structures 

Chapter 1. Introduction and Summary 

This report presents the results of the development and simulations of 
algorithms for the static shape determination and shape control of large space 
structures (LSS) • Observations of positions on the structure, and actuators 
for subsequent shape control, are assumed located at a relatively few 
discrete points along its surface. 

Quadratic performance criteria are defined to provide a means of 
determining **best** shape estimates and control forces. The resulting 
constrained optimization problems are solved using an integral operator 
approach, which proves ideal for the mixture of continuous and discrete 
problem elements. 

Results are illustrated in the one dimensional case with a flexible 
beam, and in the multidimensional case for a large space antenna. 

1.1 Background 

The development of the space shuttle has made it possible to design 
space structures larger than ever before, which may be carried into space 
and deployed or assembled tliere. Examples of such structures include the 
space platform, which would support experiments, laboratories, observation 
instruments and even habitation modules, and the solar power satellite, 
which would collect and transmit solar energy. 

Large space antennae, ranging in diameter frem 50 meters to one kilometer, 
are also being planned. They will assist in earth communications, radio 
and high energy astronomy, the deep space network as orbital relay antennae, 


and the remote sensing of soil moisture, salinity concentration and climatic 
conditions on the earth. The latter Information would assist agricultural 
productivity around the world. 

Satisfactory performance of these large space structures depends 
upon the competence of their control systems. Three kinds of control systems 
must be developed: shape, attitude, and orbit transfer and stationkeeping. 

In the past, the major deleterious Influence on shape was the Interaction 
between the control system, or systems, and the structural dynamics of the 
spacecraft. Such Interactions were minimized at the design stage, by 
guaranteeing a large separation between the modal frequencies of the struc- 
ture and the control system bandwidth. This is accomplished either by 
stiffening the structure, which Increases its natural frequency (and often 
Its weight)., or by reducing the control system bandwidth, which usually 
reduces the control system performance. 

however, in the case of the space structures now being designed, the 
enormous size, coupled with shuttle payload considerations, requires the 
use of lightweight, flexible materials. On the other hand, the performance 
criteria are extremely stringent. Furthermore, other influences, in 
particular gravity and temperature gradients, will exert significant torques 
on the structure. Thus design considerations are no longer adequate for the 
maintenance of appropriate shape. 

The sliape control problem is actually the dual problem of shape 
determination followed by shape control. Shape determination must be 
accomplished by the processing of possibly inaccurate observations of a 
number of predetermined positions along the structure. After the shape is 
estimated, shape control must be accomplished by means oi actuators (control 
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devices) placed at a finite number of discrete (Isolated) points » which 
produce forces or torques in one or more directions at these points. Since 
the sensing devices and actuators are likely to be both expensive and heavy, 
in comparison with ocher structural elements, they will be limited in 
number and in the choice of their positions. 

Thus we require methods for determining and controlling the sacipe of 
continuous structures by means of discrete or pointwlse observations and 
control devices. This is referred to as the continuous-discrete nature of 
the problem. . 

Within shape control four categories have been identified; dynamic | 

shape control (control of active vibrations), static shape control, model j 

verification, and engineering verification. This report deals with the 
problem of static shape control for large space structures. 

1.2 The Model 

In formulating the general system model it is helpful to consider the 
shape of the dish of a large space antenna. Its ideal or rest shape is 
a parabolic shape embedded in three dimensional space. If P is a point 
on the rest shape, the shape of a distorted antenna may be described by a 
three or six dimensional shape function U(P), which represents the transla- 

3 

Lional and/or rotational displacements in R of the distorted shape from 
the ideal shape. 

Thus we consider an n dimensional state function U(P) , defined on a 

I 

simply connected domain e R . We assume the state is governed by linear 
dynamics 

h U(P) * F(P) for P € Q, 0) 


where L is an n x n matrix of differential operators. 
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Associated with the dynamics <1) is a set of linear homogeneous 
boundary conditions 

Bj^(U) » 0 , i 1 i 1 

on r, the boundary of fl, which will determine the number of degrees of 
freedom of the antenna as a whole. The conditions (2) may represent portions 
of the boundary which are pinned « simply supported, or free. 

We will assume the system (1-2) is self-adjoint. 

The n dimensional vector function F(P) in (1) represents forces or 
torques acting on the system. In the shape estimation problem, F represents 
the unknown forces producing the shape distortion. F is to be determined, 
along with the shape itself, by means of a set of, possibly inaccurate, 
observations 

Yf » U(P^) + Z^, 1 ^ i ^m , (3) 

of the shape at the ra positions . 

In the shape control problem the vector F has the form 
m 

F(P) - [ OF 6(P-P.) . (4) 

i=l 

The representation (4) for F corresponds to the assumption that the forces 
are to be applied in one or more dimensions at the positions P^. A 
force applied to a rotational coordinate is a torque. 

To provide a measure of the optimal estimates of the shape and distur- 
bance functions, or alternatively the optimal set of control forces, we 
will define quadratic performance criteria. 


Thus the shape determination and shape control problems become 
constrained optimization problems, consisting of the following problem 
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elements: A continuous state which satisfies a self-adjoint linear boundary 

value problem, together with a set of m observations or forces applied at 
discrete points on the structure, and a quadratic performance criterion, 
which Includes both continuous and discrete components, and serves as a 
means ot comparison of estimates or control forces. 


1.3 Approach and Procedure 


Wc will apply an integral operator approach to the solution of both 
the static shape determination and shape control problems in the following 
manner: for a given forcing function F, the solution U of the boundary 

problem (1-2) may be expressed in terms of an Integral operator K: 


U(P) » KF « 


G(Pjij) F(Q) dQ 


(5) 


where the function G(p|q) is the Green's funotlot^, or influence coefficient, 
corresponding to the system (1-2). The integral operator K in (5) represents 
the inverse of the operator L on an appropriate space of functions. The 
use of the integral expression (5) in place ol the differential boundary 
value problem (1-2) eliminates some or all of the constraints in the 
optimisation problem, and proves particularly advantageous in the ease of 
a continuous-discrete problem mix. 


r roceduro 

We will begin by solving the static shape control and estimation 
problems for a one-dimensional shape function u, in Chapter 2. The results 
will be illustrated in Chapter 3 by simulations of a flexible beam, for 
both simply supported and pinned-free boundary conditions. 
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Consideration of the one dimensional case has several advantages: 

It is easier to use intuition about the results, and it is possible to be 
specific about the identity of the operator L and its inverse K* Thus 
exact solutions may be computed, and compared with solutions from modal 
approximations of the type which must be used in the multidimensional case* 

In Chapter 4 the results derived in Chapter 2 are applied to the case 
that L is a partial differential operator. The static shape distortion 
of a circular membrane and a rectangular plate are considered as examples. 

The analytical results are similar to those for an ordinary differential 
operator, but it is clear that even when the operator L is known, the 
specific Greenes function for a system governed by a partial differential 
equation may be difficult or impossible to compute. Approximate algorithms 
using the system modes (eigenfunctions) , which can still be computed 
analytically, are also presented. 

In Chapters 5 and 6 raultidimensional shapes, corresponding to most 
LSS models, are considered. In Chapter 5 the theory is developed. It 
parallels the theory for the one dimensional case, with some exceptions. 

The differential operator and the Greenes function are matrix operators. 
Observations and control forces may be applied to only some of the components 
of the state at each point. Furthermore, in most cases the differential 
operator L and the system modes are not explicitly known. Thus the modes 
must be computed experimentally, or by a modeling method such as the finite 
element method. Approximate solutions based on eigenfunction expansions 
corresponding to the static model are presented. 

In Chapter 6, in order to apply results to a finite element model of a 
large space antenna, the methods of Chapter 5 are adapted to the of 
eigenfunctions supplied by a dynamic (time-varying) model. A dc.^mpuion 
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of the finite element method is presented. The control problem is used to 
demonstrate the exact correspondence between solutions of the continuous 
static problem and the finite dimensional static model of the finite element 
method. Finally, results are Illustrated by simulations, using data from 
a finite element model of a large space antenna. 

Conclusions and future work are stated in Chapter 7 . 

The appendices include program listings and outputs for the simulations 

of the flexible beam (Appendix B) and the LSS antenna (Appendix C) . 

Appendix A contains a simplified sketch of distribution theory, the 

mathematical theory within which the use of the delta "function" may be 

considered legitimate. It also contains a proof of the identity of the 

free space solution of v\ 6(P-Q) , which is a part of the Green's function 
4 

for the operator V . 


1.4 A Comment on the Approach 

The integral operator approach Is ideally suited to the continuous- 
discrete problems of LSS shape control and determination. Physically the 
Green's function represents the response of the system to a unit impulsive 
force at one point. Thus, the shape control problem, for example, becomes 
merely the problem of determining the linear combination of Green's 
functions or responses at each point which produce the best approximation to 
the desired shape. 

The analytical problem of handling a continuous-discrete mathematical 
mixture can prove messy or awkward. The integral equation approach reduces 
the elements of the shape control and determination problems either to purely 
discrete or purely continuous problems which are more easily handled. 
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In addition, no approximations, other than the initial assumptions of 
linearity and polntwise application of forces or observations, which are 
common to most engineering approaches, are applied until the final computa- 
tion of the solution algorithms. This approach has value in both its 
simplicity and its generality. Intuition about the behavior of the system 
can be retained to the final computation stage. 

For example, it is easy to determine the additional constraints which 
must be applied in the case that the system has rigid body modes (eigenfunc- 
tions corresponding to zero frequencies), and to understand their physical 
interpretation. 

Furthermore, the shape control and estimation algorithms are not dependent 
on a particular model, since the only dynamical assumptions are that the system 
is linear and self-adjoint. A change in the model does not necessitate a change 
in the method, only a change in the eigenfunctions used to approximate elements 
in the algorithms. The eigenfunctions can be provided by lumped mass finite 
element models, which are themselves linear and self-adjoint. 

Finally, the use of integral operators rather than differential ones 
possesses these general advantages: 

(1) The expression of a solution as an integral equation automatically 
incorporates the boundary conditions, which must be stipulated separately 
if the problem is stated as a differential equation. 

( 2 ) The integral operator is usually bounded and often completely 
continuous, whereas differential operators are unbounded. Thus results 
concerning eigenfunction expansions, solutions of nonhomogeneous equations 
etc. are more easily obtained. 

(3) Numerical approximations and variational techniques which include 
several other methods ot solving problems with constraints are more 
easily applied to integral rather than differential equations. 
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Chapter 2. Static Shape Control in One Dimension 
2.1 Introduction 

In this chapter we present the general theory for a one dimensional shape, 
which will be illustrated by a flexible beam model in Chapter 3. While 
Che sK^pe of a large space structure is usually modeled as multidimensional , 
consideration of Che one dimensional case possesses several advantages: 

1) It is possible to be explicit about Che idenClCles of Che differential 
operator L and its inverse, the integral operator K. Thus exact solutions 

to the shape determination and control problems may be computed. 

2) Intuition about the physical meaning of results may be applied more 
easily to the one dimensional case. 

Procedure 

In section 2.2 we define the general linear boundary value problem 
(BVP) satisfied by a one dimensional shape function u, and discuss the 
existence of solutions. In section 2.3 we define the corresponding Green's 
function, and demonstrate its role in the solution of the BVP. We discover 
a mathematical distinction between the problem of shape control and those 
of attitude control and stationkeeping. 

We will state general shape control and determination problems for a 
one dimensional state in section 2.4 and 2.5, and use the Green's function 
to derive algorithms for their solution. 

In section 2.6 we will present eigenfunction expansions which may be 
truncated to provide approximations to elements of the shape control and 
estimation algorithms. Since in the multidimensional case approximations 
must be used, it is Interesting to compare them to the exact solutions 
available in the one dimensional case. 

Conclusions are stated in section 2.7. 


10 


2.2 The Boundary Value Problem 

j9 

Consider a surface which occupies a simple connected region 0 e R and 
is bounded by the curve F . 

Assume the surface is acted on at each point P e 0 by a force f(P), and 
that the static deformation u(P) of the surface satisfied the partial 
differential equation. 

Lu “ f (6) 

where L is a linear ordinary or partial differential operator, related to 
the stiffness of the structure, which also satisfies linear boundary 
conditions 

B^(u) » 0 , 1 ^ i ^ k^, for P e r. (7) 

Assume the bounda.'y conditions (7) are such that the operator L is self- 
adjoint. That is 

<Lu, v> = <u,Lv> (8) 

for any pair of functions (u,v) in an appropriate class which satisfy the 

boundary conditions. (The term "appropriate class" is purposely vague. 

See Appendix A.) The inner product <u,v> is defined to be the integral 
» 

<u,v> « u(Q) v(Q) dQ . (9) 

Jr. 

Solutions of boundary value problems do not always exist. Before the 
Green’s function can be defined and its role in the solution of (6-7) 
discussed, it is helpful to recall the following rule from linear differen- 
tial equations, which gives sufficient reasons for the existence of a solution: 
Consider the self-adjoint boundary value problem (6-7) and its corresponding 
homogeneous problem 

Lv * 0 , B^(v) « 0 , ^ 1 i 1 (10) 
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Then (a) The system (6-7) has a unique solution for each f If and only 

If the homogeneous system (10) has only the trivial solution. 

(b) If (10) has non-trlvlal solutions, the problem (6-7) has no 

solution unless the consistency condition 

<f,v> ■ I f(Q) v(Q) dQ « 0 Cll) 

■'fl 

is satisfied for every v(P) which Is a solution of (10) . This rule Is a 
simplification of Theorem 5.1 In Chapter 5. 

Remark 2.1: If a solution u(P) of (6-7) exists, and v, ,...,v_ are independent 

' ■ ^ g 

non-trlvlal solutions of (10), then u is not a unique solution, since 
s 

u + I c V. (12) 

i»l ^ 

is a solution of (6-7) for any set of constants c^. 

Remark 2.2 ; The consistency condition (11) becomes reasonable when we 
consider that seeking a solution to (6-7) for any function f in some space 
is equivalent to seeking the Inverse of the operator L on that space. If 
the null space of L is zero (l.e. the solution of (10) is only the trivial 
solution) then L is one to one and its Inverse may be defined. If (10) has 
non-trlvial solutions, L is not one to one and L ^ may be defined, if at 
ail, not uniquely on the range of L. The "consistency condition" guarantees 
that f has uo component in the null space of L, hence (with a little more 
work) that it is in the range of L. 

2.3 The Green's Function 

We first consider case (a) of the rule in the previous section. The 
corresponding homogeneous problem (10) has only the trivial solvit ion. Then 
the Green's function for the problem (6-7) satisfies 
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Lg(P|Q) - fi(P-Q) for P, Q e 11 , (13) 

B^(g) - 0 , 1 1 1 1 for Per. (14) 

It represents the response of the system at the point P to a unit impulsive 
force at Q. £(P<^) is the dirac delta function. 

Since L is self-adjoint, and both u and g satisfy the boundary conditions, 
we have 

<u,Lg> » <Lu,g> (15) 

which implies that 

u(P) » I u(Q) d(P-Q) dQ » I g(PlQ) f(Q) dQ . (16) 

•' •'ll 

Remark 2.3 ; Because the BVP (6-7) is self-adjoint, g(P|Q) is symmetric, 
that is g(P|Q) “ g(Q|P). [2] This is proved in the multidimensional case as 

Theorem 5.2 in Chapter 5. 

Remark 2. A ; The Green's function is the kernel of the compact Integral 
rator K such that 

Kf » f G(p|q) f(Q) dQ . (17) 

•’ll 

K is clearly the inverse of the operator L, where defined on the range of 
L , s ince KLu “ Kf ■ u and LKf » Lu = f . 

Remark 2.5 ; The solution of (13) is called a fundamental solution . The 
equation (13) is satisfied in a distributional rather than a polntwise 
sense. That is 

<Lg,v^> » <G,L*<|..> « <1,(0 (18) 

for all test functions (A test function is an infinitely differentiable 
function defined on R which has compact support. See Appendix A*) 
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The Modified Crcen*w Function 

We now consider esse (b). Suppose, the problem (10) has s independent 
solutions which we assume have been made orthonormal with respect 

to the inner product (9). We may not define the Green's function as in (13-14) 
bv'icausc 


<6(P-Q), Vj^> » j Vj^(Q) 6(P-Q) dQ ■ v^(P) 0 . 

(19) 

Thus the consistency condition (11) is not satisfied, 
the modified Green's function g(pjQ) wnich satisfies 

Therefore, we define 

Lg(P|Q) • «(P-Q) - 1 v.(P) v.(Q) 

1 

(20) 

B^(g) - 0 , 1 1 i • 

(21) 


We have subtracted the offending components of 6(P-Q) which lie in the 
nullspace of L. A solution to this system does exist. It is not unique, 
however, since the addition of any linear combination of the solutions 
Vi,...,v^ is also a solution of (20-21). We therefore Impose an additional 
constraint on g: 

<^8(p|Q) . v^> ■ 0 , 1 1 i 1 ® • ^22) 

The function which satisfies (20-22) is the unique Green's function of 
minimum norm, that is, the Green's function which Itself has no component 
in the nullapace of the operator L. 

We apply the relation (IS) to the modified Green's function. We 
note that 

<g, Lu> - I g(P|Q) f(Q) dQ 
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and 

<u, Lg> ■ f u(Q) (6(P-Q) - I v.(P) v.(Q)) 

i 

■ u(P) ■ I (f u(Q) V (Q) dQ) V. (P) 
ijl ^ 

• u(P) - I C. V. (P) . 
i-1 ^ 

Thus 

u(P) - f g(P|Q) f(Q) dQ + f c. v.(P) . (23) 

Jq 1-1 ^ 

The arbitrary constants are an expected consequence of Remark 2,1. For 
reasons given in the next segment, we may neglect the last term of (23). 

Rigid Body Modes 

As will be seen in the examples, the solutions of the homogeneous 
BVP(IO) are the rigid body modes, or degrees of freedom, of the system. 

They represent changes in position the structure may take as a rigid body. 

The pinnedofree beam in section 3.3 has one rigid body mode: it may rotate 

about the pinned endpoint. 

If a structure has free-free boundary conditions., which represent 
a structure floating freely in space, it may rotate or translate without a 
change in its shape. In three dimensions this Implies up to six rigid 
body modes. 

If the boundary is firmly fixed, the structure will have no rigid 
body modes. This is the case with the simply supported baom in Chapter 3, 
the distorted membrane and plate of Chapter 4, and the large space antenna 
with fixed hub in Chapter 6. 

Since shape distortion is measured with respect to the structure itself, 
it Is r?asonable to define a structure-centered coordinate system: the origin 

and axes are defined to be along the structure. To such a coordinate system 
the rigid body modes are invisible, and the solution of (6-7) for case (b) 
becomes (16), as for case (t>). Since the constants in (23) are arbitrary, 
no generality is lost in this assumption. 
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The consistency condition (11) will be seen to imply that no net 
forces or torques may be applied in the direction of any degree of freedom* 
Were this not so, an acceleration would result, contradicting the assumed 
boundary conditions* 

This condition (11) , coupled with condition (22) on g and the arbitrary 
constants in (23) , imply that the rigid body modes are both invisible to the 
shape control system and beyond its powers of influence* Translational 
and rotational motions must be controlled by the other control systems* 
Attitude control, orbit transfer and stationkeeping* This is the mathe- 
matical distinction between the systems mentioned in section 2*1. 

2*4 The iitatic Shape Control Problem 

In this section we define a general shape control problem for one 
dimensional shape functions* We first solve the control problem assuming 
case (a) of the rule in section 2*3. We then discuss the solution for case 
(b)^ which is slightly more complicated, due to extra constraints imposed by 
the consistency condition* 

We assume the control devices are located at the points P^, 1 < i < m, 
along the structure. The general model for the control problem is 


m 


Lu = ^ f 6(P-P ) 

(24) 

i=l 

(u) = 0 , 1 1 j 1 

(25) 


where u(P) is the shape, L is a linear differential operator as before, f^ 
is a force to be applied at the point , and (25) denotes an appropriate 
set of boundary conditions* 


Let \p be the desired shape of the space structure* Define the 


criterion 

J(F,u) = i f 2 ^ 1 (^(Q) _ u(o))2 dQ 

i=l ^ Jq 


( 26 ) 
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as a measure of performance. The constants are arbitrary weights and 

F ■ (f. . . .f )^. 

1 m 

The control problem is to determine the vector of forces F'fc which 
together with the corresponding solution u* of (24-25) minimizes J overall 
admissible sets (F,u). 


Solution of the. Control Problem 

There are two basic approaches to constrained optimization problems. 
One is to use Lagrange multiplier theory. We will use this method to solve 
part of the shape estimation problem. 


The other, perhaps more direct method, is to solve the constraints 
for an expression for some of the variables in terms of the others. This 
expression is substituted into the function of fewer variables, which cun 
be minimized without constraints. 

We will use the second approach in the control problem. We first 
assume the system has no rigid body modes ; 

The solution of (24-25) is given by 
r m 


u(F) 


rO’Iq) I I <S(l>,-g)l dQ 
hi 1=1 ^ ^ 


m 

1 

1=1 


g(F|l*j) 


(27) 


where g(PlQ) satisfies (13-14). Substitution of (27) into the criterion (2b) 


yield.- 

.1(F) = I f " r + V (ij.(g) - 1 f . g(g|p.))‘ dg . 
i=l ^ ‘ •'w i=l ^ ^ 


(28) 


The constrained optimization problem i24-2t') has become the simpler problem 
of ‘’linimtzing a tunctiun of m uukuowi^ constants without constraints. 
Simultaneous solution of the equations 
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X 

leads to the following necessary condition for an optimal solution 

F* - (f / ... 

I m 

(R + A) F* • B (30) 


The m X m matrices R and A have coefficients 
Ajj - Bd’^lO) gU’jlQ) JQ 


(31) 

(32) 


(33) 


and the m dimensional vector B has coefficients 

B. - I bU’JQ) ♦(Q) 00 . 

Once the optimal forces are determined^ the optimal shape is given by 
(27) s 


Solution of the Control Problem: Case (b) 

Wo assume that the homogeneous BVP corrospondii\g to (,24-25) has s 

indepondont solutions This is^ of course ^ oqulvalont to the 

assumption that the structure governed by (24-2.5) has s rigid body miuJes. 

In order for a solution to (24-25) to exist, the consistency condition 
v\ m 

0 - "V I f. 6 0M' )> - I t V, U’.) 04) 

•' i«l 1*1 ^ 

must Ih’ satisfied tor each function Vj. Thus the control problem Is to 
detxnmlne the sot of forces and shape I unction u which minimize the 

criterion v2n) subject to the constraints 124-25) and i34) . 

We will assume the coordinate system Is centered on the apaceoralt 
Ireeall the segment "Rigid Body Modus"). The solution ot (24*-25) Is given by 
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U(P) - I ^ g(P|P4) (- 

i-1 ^ ^ 

where g Is the modified Green's function which satisfies (20-22). 

We first solve the s constraints (3A) for the forces f^,...,fg in 
terms of the remaining forces 


^i “ I » 

j-s+1 J 


1 < 1 < s . 


It is clear that a necessary condition for any solution to exist is that the 
number m of forces applied must be at least as great as s, the number of 
rigid body modes. If we wish to obtain an optimal solution m must be 
greater than s, since for m = s the condition (34) determines the forces 
uniquely. 

Substitution of (36) into (35) yields 


m 

u(P) » I (g(p|p.) + I C g(p|p.)f.) 

i“s+l ^ i-1 J 


De f ine 


(P) » g(P|PJ + I c g(P|P.) 

j=l ^ ^ 


u(P) * I Yj(P) O'- 

i-s+1 

We substitute expressions (36) and (39) into the performance criterion, 
which results in 

i s m rt m 


^ six ' xLx 


+ (^(p) - I yAv) f,)*- dp . ( 40 ) 

*■ ■'0 s+1 

The criterion is now a function of the (m-s) constants f f , without 

S+1 m 

constraints. J is minimized by solving simultaneously the (m-s) conditions 


0 , i » s+i , . . . ,m . 


19 


f 

\ 

\ 


Let F and B be (m-s) dimensional vectors with components 

•k 


F = f 
1 ^1+s 




'l'(p) dP 


and the (m-s) square matrices R and A with components 


'ij 




^Ij “ ^ *^1 + L • 


(42) 

(43) 

(44) 

(45) 


n 

Then the optimal control law for the control problem (24-26) (34) Is 


(R + A) F » B . 


(46) 


Once the optimal forces » • • • determined from (46), the optimal 

* * 

forces fj^,...,f^ may be found from (36), and the resulting optimal shape 
is given by (35) . 

The non-constant terms in A and B are linear combinations of terms of 
the form (32) and (33) respectively. 


2.5 The General Estimation Problem 
For the estimation problem we assume the shape u(P) satisfies the 
boundary value problem 

Lu = f , B^(u) = 0 , 1 i i 1 . (47) 

where f(P) is an unknown function representing disturbances or inaccuracies 
in the model. Sensors placed at the positions P^» 1 i tn, yield the 
observations 

= u(pp + (48) 

where 4^ is an unknown constant representing inaccuracy in the observation 

at P , . Let Z = (4, ... 4 ), We define the performance criterion 
1 i m 
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I(Z,f) « y ? C ^ r + I [ f^(Q) dQ . 

^ i-1 ^ ^ ^ J a 

•I I (yi-«(Pi))^r.-^ + f ( 

^ i»l ^ 1 Jq 


(Q) dQ . 


(49) 


The estimation problem Is to determine the pair (u*,f*) which jointly 
satisfy (47-48) and minimize the criterion (49) over all admissible pairs 
(u,f). 

Solution of the Estimation Problem: Case (a) 

We assume there are no rigid body modes. Then the solution to (47) 
is given by 

u(P) = f g(P|Q) f(Q) dQ (50) 

Jsi 

where g(P|Q) again satisfies (13-14). Thus 

u(P.) = f g(PjQ) f(Q) dQ . (51) 


We substitute (51) into the criterion (49), which produces the criterion 


1 m 

'(f) ‘ J I (Vi - 
“ i=l 


g(P,|Q)f(Q) dQ)‘ r 


f 


r(Q) dQ . (52) 


The problem is now to minimize the functional J without constraints. A 
necessary condition for a minimum of J at is that the differential 

3J(f*,h) = 0 = J r7^y.,-f g(Pj(Q)f*(Q)dQ)(- f g(PjQ)h(Q)dQ) 

i=l ^ ^ ^ ^ 

+ I f*(Q)h(Q)dQ . (5 

•’i'2 

for all admissible variation h. (The unknown noise function f and variation 
h may be assumed to be in L 2 (*i), for example.) Thus it may be concluded 


(53) 


that 


m 


-1 


f*(P) = I r. g(P!P.)(y. - u*(p )) 
i=l 


(54) 


Substitution of this relation into (50) yields the optimal shape estimate 


m 


u*(p) “ I [r ^(y - u*(P )) { g(PlQ) g(P.|Q) dQ] . (55) 

i=l ■'W 



Note that u*(x) is expressed In terms of the unknovm discrete shape 


estimates u*(Pj^). Let 


and 


X - (u*(P.) . . . u* (P^))^ 
1 m 




Evaluation of (55) at x = , j * 1,.. 

condition for the vector X: 


(I + AR~^)X - Y 


(56) 

(57) 

m yields the following necessary 

(58) 


where A and R are the matrices of coefficients (31-32). 

Once the vector X has been determined the optimal shape estimate 
is given by (55) . 


Solution of the Estimation Problem: Case (b) 

We now assume the structure described by (47) has s rigid body modes 

V. , . . . ,v , which are orthonormal with respect to the inner product (9), 

X s 

The estimation problem is to determine the pair which minimizes the 

criterion (49) over all admissible pairs (u,f) which satisfy (47) and the set 
of consistency conditions 

1 < j < s . (59) 






f (Q) Vj(Q) dQ - 0 , 


We will show that the solution of the estimation problem for case (b) has the 
same form as tliat for case (a): 

The solution to (47) is given by (60), where g is the modified Green’s 
function (20-22). 


u(P) - 




g(PlQ) f(Q) dQ 


(60) 


We evaluate (60) at 1 i and substitute into the criterion (49) 

producing the criterion (52). Thus we have eliminated part of the constraints. 
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the boundary value problem (47). The estimation problem becomes the 

problem of minimizing the criterion (52) subject to the remainiug constraints 

(59). 

We will apply the Lagrange multiplier theorem [3] . We adjoin the 
constraints to the criterion by means of scalar multipliers 

J » 7 ! (V. - f g(PjQ) f(Q) dQ)2 r,"^ + j f f^(Q) dQ 
^ Jo ^ ^ Jo 


1»1 

s 


J-1 ^ 


f(Q) Vj(Q) dQ . 


(61) 


A necessary condition for a minimum of J at f* is that the differentials of 
J with respect to f and 1 j ^ s, are 0. We have 


3J 


3X 


= 0 


j 


f(Q) v^(Q) dQ , 


1 1 j 1 s , 


(62) 


and 


m . C 

3J(f,h) = I (y^ - u(P^))r^~-' (- g(PjQ)h(Q)dQ) 


i=l 




f (Q) h(Q) dQ + I X 

i=l ^ ■’ 


j= 

We factor the variation h to one side: 


h(Q) Vj(Q) dQ * 0 


(63) 


h(Q) dQ (-(y - u(P ))r fc(PjQ) + f(Q) + ! X v (Q)] - 0 

j*l J J 

Since this must be true for all admissible variations h, the bracketed 
term must be zero. We have 


f(Q) - 


1 

(y - u(P ))r g(P. I q) - I V (Q) . 

j=l J J 


(64) 


We apply the other necessary conditions (hil) . 

m 

r ^ t\f — 11^1 

i 


0 = ■ u(pj)( 


i«l 

s 


SI 


g(PilQ) Vj^(Q) dQ) 


- I , V. > , k » 1,..., s . 

j»l J J 


( 65 ) 



The set was chosen orthonormal, so ■ 6(J-k). 

Furthermore, from condition (22) on the modified Green's function we 

know that 
» 

g(PjQ) V. (Q) dQ « 0 , k « 1,..., s . (66) 

Jjj 1 

Thus we have ® 0 for k » s . 

We may conclude that 

m . 

f*(Q) “ g(PilQ) “(Pi)) (67) 

as .In case (a), and therefore that the optimal shape estimate u* Is also 
given by (55) . 

Remark 2.6 ; Note that because of condition (22) on g, the optimal shape 
estimate has no component In the direction of the rigid body modes. There 
may be components In the actual shape, but a shape control system has no 
means of determining them. 


2.6 Approximations 

If the Green's function Is known, the shape determination and shape 
control problems mny be solved exactly by the methods of this chapter, 
however, It will be seen In Chapter 4 that when L Is a partial differential 
operator It can be difficult to determine the Green's function. For large 
space structures, which are multidimensional, the determination of the matrix 
differential operator L, and consequently the Green's function. Is usually 
impossible. 

However, the Green's function, and the terms in the shape control and 
determination algorithms which involve the Green's function, may be expressed 
in terms of series expansions Involving eigenvalues and eigenfunctions 
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corresponding to the BVP (6-7). Truncations of those series can serve as 
approximations of the relevant terms. Even when L Is not known the eigen- 
functions and frequencies can be computed numerically, for example by the 
finite element method. 

Let normalized eigenfunctions of the boundary value 

problem (6-7), corresponding to the non-zero eigenvalues ^ 2 * **** 


{((ij} and {Xj} satisfy 
L (P) » Xj (P) 


for Pen, 


« 0 , 


1 < 1 < k for Per. 
— — 0 


Eigenfunctions corresponding to zero eigenvalues are rigid body modes. 
We have the following expansions: 

g(P|Q) - I f-*.(P)*.(Q) (70) 

i i ^ ^ 

and 

g(P|Q) f(Q) dQ - ^ f-^.(P) . (71) 

Jq j k ^ J 

Substitution of (70) for f in (71) yields 

g(P|Q) g(Ql») dQ = I -V • <72) 

Jn j x/ J J 


The expressions (71) and (72) provide approximations for the terms and 
defined by (33) and (32) in the control and estimation algorithms. 

The series expansions (70-72) are standard results of linear operator 
theory (2]. They are based on the assumptions that the integral operator 
K defined by 


g(P|Q) f(Q) dQ 
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Is a symmetric Hllbert-Sr i'tir’.t operator. The symmetry follows from the 
self -adjointness of the boundary value problem. Au operator K Is Hllbert- 
Schmldt If 

I|k|| - (f f |g(P|Q)|^dPdQ)^/2 < » . (73) 

Jn Jii 

In the case that L Is an ordinary linear differential operator, 
as in Chapter 3, the Green's function Is continuous on the compact domain 
ft, which Implies (73). If g is not known precisely, the property (73) 
must be assumed. 

2.7 Conclusions 

An integral operator approach to the continuous-discrete optimization 
problems of static shape estimation and control proves ideal for these 
problems. Solutions reduce to the solution of linear equations of dimension 
less than or equal to the number of observations, or control forces. 

A distinction must be drawn between the solutions for systems with 
rigid body modes and those without. The cotitrol law for a system with rigid 
body modes is more complicated, due to the imposition of extra constraints 
on the forces, which represent the requirement of zero net forces and/or 
torques in the directions of these modes. 

The estimation procedure for a system with rigid body modes is the same 
as for a system without them, but the resulting estimate has no component 
in the direction of the rigid body modes, because they are invisible to the 
shape estimator. The rigid body modes represent changes in attitude and 
translational movement, which must be the concern of the attitude control, 
orbit transfer and stationkeeping systems. 
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In the event that the Green's function cannot be precisely kno«m, 
approximations to the terms in the control and estimation algorithms may be 
computed from eigenfunction expansions available from linear operator theory. 
The eigenfunctions, often called modes or mode shapes, may be computed 
numerically even when the operator L is not known. 
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Chapter 3. Static Shape Control for the Flexible Beam 
3.1 Introduction 

A flexible beam provides a perfect illustration of static shape distor- 
tion and subsequent shape control. Consider a flexible beam which is 
supported at the end points, and is intender> to serve as a bookshelf. The 
desired shape, or rest shape in the absence of outside forces, is strictly 
horizontal. However, the forces of gravity act continuously along the beam, 
causing it to sag in the center. 

In order to achieve the desired horizontal shape, we apply a third 
support under the center of the beam. The natural stiffness of the beam 
together with the appllcatl^^^* of this additional force at the center approxi- 
mately counteract the effects of the gravitational force. Thus we observe 
static shape control by means of one polntwise force in the ordinary brick 
and board bookcase. 

In this chapter we will solve static shape control and estimation 
problems for a flexible beam of length 1, and boundary conditions representing 
simply supported, or pinneti-free endpoints. 

3.2 Shape Control for a Simply Supported Beam 

Consider the problem of controlling the static deflection of an elastic 
beam of length I, Define a coordinate system such that the x-axls passes 
through the endpoints of the beam, with one end at the origin and the other 
at X « £. Suppose control Is to be Implemented by means of transverse forces 
f , at positions x, , 1 < 1 < m, where 0<x, <x«...<x <t. See Fig. 3.1. 

At each point x t[0,£] denote the deflection by u(x)* Assuming no net 
tensile force on a cross-section, the shape of the beam is governed by the 


differential equation 


~ • I t. «(x-x ) (74) 

dx^ 1-1 ^ ^ 

The ends of the beam satisfy the boundary conditions 

u(0) - u'*(0) - 0 u(i) » u”(f) - 0 . (75) 



Figure 3.1 The Simply Supported Beam 

Let <j/(x) be the desired shape of the beam. As a measure of performance 
we define the criterion 

J(u,F) - f ’'l Ho 

T 

where F is the vector of forces (f, .♦* f ) and r. are non*-negative 

1 m X 

constant weights whose values are optional « 

The object is to determine the set of forces f^* t’hich together 
with the solution u*(x) of (74) minimizes (76) over all possible pairs 
(u,F). 

The existence and uniqueness of a solution to (74-75) follows irora 
the fact that the associated homogeneous system 
.4 

- 0 , v(0) » v"(0) - 0 , v(£) » v'*a) - 0 (77) 

dx 
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has only the trivial solution. Consequently the solution of (7A-75) Is 
given by 

m 

«(x) » I g(xlx.) f. (7f 

1-1 ^ ^ 

where g(x|^) Is the Green's function which satisfies 


= 6(x-^) 


g(ojo = g"(ojo = 0 , sale) - g"aU ' o • (8oj 

The Green's function represents the response of the beam shape to a unit 
impulsive force at x “ 

The solution of (79-80) is 


g(x|C) = 


a-SL)x ,1 


(x‘ - 21£; + 


(x2 - 21x + Kh 


0 < X < C 


C < X < £. 


Figure 3.2 displays the Green's functions which correspond to impulsive 

£ 

forces at positions C * n (g) , n = 1, 7. 

The solution of the control problem: Substitution of the solution 

into the criterion (76) yields 


H m 


19 II 9 

J(F) = ^ I fi i (I g(xix^)f^ - (P(x))^ dx ( 82 ) 

i-1 -'o i==l 


The problem of minimizing the criterion (76) subject to the constraints 
(74-75) has become the problem of minimizing a function of m unknown 
constants without constraints* A necessary condition for J to have a 
minimum at F* is 


(F^) ^ 0 


1 < i < m 
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This condition becomes 


m fl fi’ 

r, + I f.( g(xlx.) g(xjx ) dx « i|)(x) g(x|x.) dx . (84) 


If we define 


a 


ij 


g(xlx^) g(x|xj) dx , j im 


and 


ij;(x) g(x, dx , 


1 ^ i ^ m , 


( 85 ) 


( 86 ) 


then the necessary condition for a minimum of J at F* is that F* satisfy 


(R + A) F* » o (87) 

where R is the m x m diagonal matrix 



A is the m X m matrix with coefficients (85), and B is the m dimensional 
vector with coefficients (86). 


The Shape Control Algorithm for the Simply Supported Beam 


1) Compute the constants and defined by (85-86). 


Define R, A, B. 

2) Solve (87) to obtain F*. 

m 

3) The optimal shape u*(x) = ^ f *g(xlx. ) . 

i»l 

Figure 3.3 displays the optimal shape vs. the desired shape 
i|'(x) = sin — the second mode of the system (74-75), for two actuators 


at 1/4''. and 3/4H.. 


ACTUATOR POSITIONS MARKED BY O 



Figure 3.3 Optimal Shape vs. Desired Shape for the Simply Supported Beam with Two Actuators 
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3.3 The Control Problem for the Plnned-Free Beam 
A modification of the control algorithm Is necessary If Che system has 
rigid body modes, as Is the case with the plnned-free beam. 

The beam with one pinned and one free end point satisfies the differen- 
tial equation (74) with boundary conditions 

u(0) = u”(0) = 0 u"(*^) = u"Xl) * 0 . (89) 


We will again use the performance criterion (76) . The object is to 
determine the set of forces {f^} which together with the solution u(x) of 
(74) (89) minimizes (76) over all possible pairs ({f^}, u) . 

The system (74) (89) has the rigid body mode x (normalized). 

Physically this means the beam can have a non-zero slope or tilt as a rigid 
body. Mathematically it means that the corresponding homogeneous system 




v(0) = v”(0) = 0 v"(£) = v"\l) = 0 


(90) 


has the non-trivial solution v^^(x). Thus the system (74) (89) has a solution 
only if the inner product 


ra m 

Mx-Xi), Vj) • ° 


(91) 


The additional constraint (91) must be added to the problem of determining 
the optimal control forces. 

A solution to (79) with pinned-free boundary conditions does not exist 
because the inner product <6(x-5), v^> is not zero. The "modified'* Green's 
function which is appropriate to the system (74) (89) satisfies 


d g (x C) 
m 

dx'^““ 


= 6(x-£.) - X4 
I 


= 8„"(0|O 0 

m m 


(92) 




( 93 ) 
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We make the additional requirement that g^(x|4) have no component In the 
subspace spanned by the rigid body mode« 


(g„(«U),Vj) - I g^(«|{) 

* it ■'o 


X dx “ 0 


(94) 


r ■'o 

The modified Green's function which satisfies (92-94) is given by 


g^(xl?) » x^ 


33t 

''140 41 


401^ 


) - 



01^15 




(95) 


Condition (94) guarantees that g^(x|^) is symmetric and of minimum norm 
among all solutions of (92,93), Figure 3,4 displays g^(x|t) for impulsive 
forces at intervals of 1/8 

The Greenes function (95) represents the response of the pinned-free 
beam to one of a set of unit impulsive forces which satisfy (91) , 

Figure 3,4 displays the Green’s functions for impulsive forces at positions 

C. 

n (g) , n = 1,,,,,7. 

The solution of (74) (89) (91) is given by 
m 

u(x) = I f R (x|x ) (96) 

i^l 1 -Tn 1 


We solve (91) for f^ in terms of the outer :orces and substitute that 
expression together with (94) into the criterion (76) , which results in 


,<p, . -- i i -± f , + 5 1 f - , 

1=2 1 1=2 

2 ^ ^ ■ 'x“ “ V'(x))'^ dx 

•’ o i=2 1 


(97) 


T 

where F is the vector (f^,.,f ) , 

^ m 

Again, the optimization problem is reduced to one of minimizing a 
function of unknown constants. 
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Figure 3. A The Green’s Function for the Plnned-Free Beam 
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The necessary condition for a minimum at F* is 
(^) «0 

These conditions result in the following algorithm. 

(i) Compute the m dimensional vector B and m x m matrix A whose 
coordinates are 


**l “ 1 «m^*l*i^ ^ 

•ij ' J <** • 

•* J 0 


(99) 

( 100 ) 


(ii) Compute the (m-1) dimensional vector B and (m-1) x (m-1) 


matrix A whose coordinates are 

. X, 


*’141 *’l 


(101) 


*i4l*1-H 


•ij“’^l *"*11> 2 

*1 


, *141 • 

Vl,j4l " *1,141 “ 1.J41 X^ 


( 102 ) 


Let R be the (m-1) x (m-1) diagonal matrix 



(iii) The vector F* of optimal forces satisfies 

(R 4 A)F* - B . 


(103) 


(104) 


The optimal force f is found from (91) . 

m 

(iv) The optimal shape u*(x) = I fj*g (x|x. ) . 

i=l 1 m 1 
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Since the optimal shape u* Is a linear combination of Green's 
functions which satisfy (9A) , It will have no component In the subspace 
of the rigid body mode. If the desired shape tj)(x) does have such a component, 
that Is If (ip, v^) Is not zero, the optimal shape will approximate the 
shape 

ij)(x) - <)|»,Vj^> Vj^(x) . (105) 

That Is , It will approximate the desired shape minus Its component In the 
subspace spanned by v^(x). 

2 

As an example. Figure 3.5 displays the desired shape <|>(x) « £x - x , 

3 2 

the shape which approximates ^ £.x-x , and the optimal shape plus the missing 
rigid body mode component ^ lx. 

Those components of the desired shape In the subspace spanned by rigid 
body modes must be added by the attitude control system. A shape control 
system constrained to satisfy the boundary conditions cannot affect these 
components. 

3.4 The Shape Estimation Problem 

To illustrate the shape estimation algorithm we consider a simply 
supported beam oi length 1 and unknown sliape u(x) , which satisfies 
,4 

— ^ = f (x) on 0 ^ X ^ 1, (106) 

dx^ 

and 

u(0) = u"(0) = 0 u(l) » u"(l) = 0 (107) 

The function f(x) represents minor model inaccuracies or random disturbances 
acting on the beam. 

Assume sensors at positions x^ , 0 < x^^ < . . , < x^ < £ , produce 


observations 




3.5 Shape Control with Rigid Body Modes 
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1 < i < m . 


(108) 


As a measure of the accuracy of shape estimates we define the criterion 

m ... ft 


1 ™ ^ —1 1 9 

J(f ,u) ® 7 I (y^ ~ ^ + 2 * 

^ 1«1 ^ ^ ^ ^ Jq 


(109) 


The object Is to determine the function f* which together with the solution 
u* of (106-107) minimizes (109) c'-’er all possible pairs (f,u). 

The solution of (106-107) Is given by 

u(x) - g(x|5) f(0 d? (110) 

^ 0 

where g(x|c) is the Green's function (81). We substitute (110) Into the 
criterion (109) ; resulting in the criterion 

m . ,i . . rl 


m 9 1 

J(f) • I (y^ - J g(xJOf(0 dir + j 


tar d^ . (Ill) 


The estimation problem has reduced to one of minimizing (111) without 
constraints. A necessary condition for J to have a minimum at f* is that 
the Frechet differential 


m . r 

8J(f,h) = I T~^ (y, - g(x U)f*(OdO(- 
1=1 1 i Jq i J 

f*(C)h(^)d^ » 0 

for all admissible variations h. This implies 
"* -1 

f*(fj = I r^“ g(x^U)(y^ - u*(x^)) . 


9. 


g(x^U)h(5)df.) 


Then 


m , 

u*{x) = I r.“'‘'(y, - u*(x )) g(x j C)g(x . | C) dC 
i=l ^ ■’o 


( 112 ) 


(ll'’l 


Let 


and 


X = (u*(x^) ... u*(x )) 

X ui 


T = (yj^ ••• y^) • 



40 


Evaluation of (113) at x « and regrouping of terms yield the following 
necessary condition for the vector X: 

(I + AR“^) X « a r“^ Y (114) 

where A Is the matrix of coefficients (85) , and R ^ Is the diagonal matrix 
with diagonal entries 


The Shape Estimation Algorithm 

(I) Compute the elements of the matrix A given by (85) , and define 
X, R, Y. 

(II) Solve the system (114) for the vector X. 

(III) The optimal error estimates are g^ven by (112) and 
?i * Yf - u*(x^), 1 ^ 1 < m. 

(Iv) The optimal shape estimate is given by (113) . 

This algorithm is equally valid for the static beam with other boundary 

conditions, provided the appropriate Green's function is used. 

Figure 3.6 displays the optimal shape estimate versus the actual shape 

, . 1 TX. , 1 ,2irx, 

sin (~) + ^ (-^) , 

11 3 

for three exact observations at t and r 

4 ’ 2 ’ 4 


3.5 Approximations 

The approximations presented In section 2.6 take the following form 
on the domain [0,£] of the x axis: 

(115) 


«(x|0 » I Y 

k»l k ^ 


j-r " 1 

ft ^ 

b, • g(xlx^) !/<(x) dx » [ 

*'c k=l ^ 


(116) 


( 11 ?) 
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where {x^) ere the ectuetor or eensor poelcione, ere the non*zero 
eigeoveluee, end ere the correeponding normellted eigenfunctions of the 
essocleCed boundery velue problcas. 

Thus for the simply supported beem the approximations based on the 
first term of each expansion are given by 




«x 


a - 2 T ala(-^) slx.(^) 


bj • 2 ^ sin(-ji) (j ♦(x) eln(p) dx) . 


(118) 

(119) 


For the planed*free beam 
wx 


yXj WX WX WX 

. _ slii(-£i) elah(--i)^ slnh(-^) 

*1J ^ ^ cos w cosh ti ^ ^ cos w ^ cosh w 


where y ■ 3.927 satisfies tan y • tanh w. 


WX 

3 sln(--i) sinhC 

b • JL. t— -.*—1 .f M 

i ^4 * cos y cosh y 


yx 

^jf‘ 


sin ^ slah(fi) 
♦ (x)(rr^+ ‘ 


fWX. 


( 120 ) 


cos y 


cosh y 


O dx 


( 121 ) 


(the normalizations are approximate). 

Approximate algorithms constructed from the first term in the eigen- 
function expansions were included In the simulations of the examples in this 
chapter. The graphs of approximate vs. optimal results were indistinguishable. 
Numerical results are included in the program outputs in Appendix B. 

It is misleading to generalize from the approximations for the one- 
dimensional cace, for which satisfactory approximations result using only 
the first teiu. The expansions (118-121) telescope rapidly, because the 

magnitude the eigenvalues increases rapidly. The frequencies of 

2 

large space structures Increase relatively slowly (X^ * ^ 

observed in the output of the shape control program for a large space antenna. 
In ApperJix C. For multidimensional structures many more modes (eigen- 
functions) must be used. 
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Chapter 4. Shape Control of Structures Governed by 
Partial Differential Equations 

4.1 Introduction 

In Chapter 3 static shape control and estimation problems for one 

dimensional cases were solved using Green's function techniquv^s . In this 

chapter corresponding results for structures defined on multidimensional 

domains, governed by partial differential equations, are presented. It 

will be observed that the solutions are very similar to those for one 

dj-.iiensional domains. The major difference is that it becomes difficult to 

determine the analytical form of the Green's function, so that expressions 

in terms of eigenfunction expansions must be used . 

We consider as examples the shape distortion of membranes and plates 

which in equilibrium position lie in a plane. A membrane, such as a drumhead, 

or the mesh of an antenna, is distinguished from a plate by the absence of 

bending resistance. The restoring forces of a membrane are due exclusively 

to tension whereas plates have bending stiffness. Consequently, membranes 

2 

may be considered to be governed by the harmonic operator V , while plates 

/ 2 2 

are governed by the blharmonlc operator 7^ * V^(7^) . 

This distinction between second and fourth order dynamics is analogous 
to the model'ing distinction between a string and a flexible beam in the 
one dimensional case. 

For convenience, in this chapter we consider only systems without 
rigid body modes. 


4.2 The Boundary Value Problem and Green's 
Function for a Membrane 

Under suitable assumptions the shape distortion of a membrane is 
modeled by the differential equation 
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V u - £(P) , 


F e ft 


( 122 ) 


where V Is the Laplacian operator and certain known physical constants have 
been Incorporated Into the forcing function f. Eq. (122) Is known as 
Poisson's equation. 

We will choose the boundary conditions so that conditluus (6-8) are 

2 

satisfied for the operator 7 » L. We will then discuss the determination 

of the Green's function g(P|Q), and exhibit the solutions to the control 
and estimation problems for the unit disk. Finally we will exhibit approxi- 
mate solutions using the eigenfunction expansions (70-72). 

Green's theorem for the Laplacian operator takes the usual form 


I (v 7^0) - 0 ) 7^ v)dP = [ 
•'ft •'ft 


, 3u) 3vv , 

(V ~ - a, ^^-) ds . 


If we impot>e either of the boundary conditions u(P) - 0 or 


3n 


(123) 


0 for per. 


the right side of (123) will be zero for functions o) and v which satisfy 

2 

the boundary condition, and the operator V will be self-adjoint. 

For convenience we eliminate the latter boundary condition, since the 
homogeneous system 


.,2 ^ 3u « 

V u *= f , V = 0 
* 3n 


has the non-trivial solution u = C. 


(124) 


The Green Function 
The Green’s function for the system 

V^u = f, u(P)»0 Pel' 
satisfies 


(125) 


^ «(^*y*^,n) * 6(x-0 6(y-n) 


(126) 
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in rectangular coordinates P(x,y), Q(t,y) and 


7^ 6(r,e.».») . 


U27) 


i0 

in the polar coordinates P ■ re , Q » pe ^ . In both cases g 3 0 on F . 


The function y 


2it 


log R, where R is the distance QP, can be 


shown to satisfy V y “ <S(p|Q). It is called the free space solution 
since it is not required to satisfy the boundary conditions. 

Thus the Green’s function is given by 

g(PlQ) = log K + g(PlQ) 

where g(P|Q) satisfies 

V^g = 0 on a, g = - log R on r. 


(128) 


(129) 


The theory of analytic functions may be applied on convenient regions to 
determine g* hence also to determine g. For ^ equal to the unit circle 

|zl i 


2 

, ... 1 , ,r - 2rp cos(0-4') + p , 

g(r,0,p,(Ji) = log [ i— _-] 

1 - 2rp co8(0-(t') + r“V 
for \\Q in polar coordinates [7]* 


(130) 


Remark A.l: Through the use of conformal mapping it is possible to determine 

the Green’s function for some ot uer regions* but in general it is not 
possible to determine the exact function g* 


4.3 Tlie Control Problem for on the Unit Disk 
The control problem for the Laplaclan on the unit disk corresponds to 
the problem of controlling the shape of a circular net or drumhead to a 
desired shape by means of pointwise forces. Thus we desire to 

determine the set of forces {fj} at positions Pj » pj ^ 1. J ^ 
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which together with the solution u<r,6) of 


5 m «(r-o.) 

7^ u(r.6) •If, ^ ^ 


i-1 


J 


u(l,6) « 0 


(131) 

(132) 


minimizes the performance criterion 

r2it rl 


J(F,u) • 4 I [ l'*'(r,8) - u(r,e)l^r dr 

j«l ^ ^ -'0 ^0 


de 


(133) 


over all possible sets (u,F), where 


F - (f- ... f ) . 

X m 


(134) 


The optimal shape for the problem (131-132) is given by 


u*(r,8) 


‘i" 


2 2 
m r -2rp cos(6-(^.) + P. 

^ I t.- 1»8 I 1 — —J J -^1 




1 - 2rPj cos(6-4ij)+r pj 


(135) 


and the vector of optimal forces F* satisfies (R+A)F* * B, where 
R = ^^ij^ m X m matrices such that 

Rij * 6(i-j) 

and 

“II g(r.6. fdrde 

and B “ (b ^) is an m dimensional vector such that 

bi • f [ <Kr,e) g(r,8,p.,<i ) rdrd6. 

'o ^0 * ^ 


(136) 


(137) 


A7 


4.4 The Estimation Problem 

2 

The corresponding estimation problem for V on the unit disk is, 
given the shape observations 

+ 4^ . 1 li Im , (133) 

at positions to determine the error function f(r,6) and corres- 

ponding shape function u(r,6) which satisfy 


u(r,e) » f(r,e) , u(l,0) - 0 


(139) 


and minimize the criterion 
m 


1 ? 

J(F,u) = I (y^ - ur .,<^^)^ r^ 


2 - -1 + 1 
* 




i=l 


f^(r,e)r dr d0 . (140) 


The results of Section 2.4 yield the optimal error estimates 


Ci* = - u*(Pi»'l>^) 


m 


£‘<c.«)- iv I 

i=l 


-1 


log 


0 2 n 

r""-2rp . cos(0-4>*) + P* 

- 1 — i 


1 - 2rpj cos(0-4>j) -f 


where the vector X = {u*(P.) . . . u* (P )) satisfies 

1 tn 


(I + AR“^) X = AR~^ 


(142) 


(143) 


The matrices R and A are as in (136) and Y is the vector of observations 
T 

(>*i ... y^) . The corresponding optimal shape estimate is then given by 
m . fl 

u*(r,0) = I [r ?.* g{r,e,p,(t») g(r,8,p ,<J) ) pdpd4)] . (144) 

i-ri ^ Jq Jo 


4.5 Approximate Solutions 

For simplicity it may be desirable to compute approximations to the 
solution (135-137) and (141-143) using eigenfunction expansions. The 
eigenvalues and (normalized) eigenfunctions corresponding to 
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♦(r.e) - X»(r,6), ♦(1.6) - 0 


are 




1/2 

o on 

VT ~ 

h <v 


11 * 1 . 2 . • • • 


corresponding to the eigenvalues which satisfy 


•^o <^n> “ 0 * 


n*!.2. •** 


and 


♦„.(r 

anc 


♦ (r 
mns 


1/2 


• 4 

\Vi «*»>/ 


• e ) - yi - ^ 


1/2 

~UT 

< U >- 


cos b6 


sin a6 
1 < a, n < 


corresponding to the eigenvalues X which satisfy 

on 

1/2 

J (X ) - 0 , 
in nn 


where 0 ^ i < «> are. of course, the Bessel functions. 

Thus, a first approxiaatlon to the forces {f^} in the control law, 

2 

using the eigenvalue X^^ ■ (2.405) and eigenfunction 

(2.405 r) 


/? < 2‘«05 1 


satisfies 


(R + A)F « B , (145) 

where R is as before, F « (f^^ f^) is the vector of approximate forces, 

A A 

and A and B are the approximate matrix and vector with coefficients 
1 


‘Ij 


4 OO 1 OO J 


OO 


(146) 
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oo 


1 1 i» j 1 ® 


(147) 


The shape corresponding to the approximate forces In (145) is given by 

o « 


u*(r.6) - I f. g(r,e,p A ) 
1,1 1 11 


(148) 


since the Green's function still represents the response to a unit force 
at Pj, » p^e^^l. 

Using the same approximations (146) for the matrix A, the pointwise 
shape estimation vector X may be approximately computed from 

(I + AR~^)X e A R~^ (149) 

where X * (u(Qj^) ... '*(0^)) Is the approximation to X, and R~ and Y are 
as in (136) (138). The approximate estimates are then given by 


= Yi - u (Pj^) , 


1 i m , 


m 


f(r) 


A i ’i 00 00 1 

00 i=l 


(150) 

(151) 


m 


u(r) = I r C. ♦ (r) ♦ (pJ 

^ ^ i i 00 00 1 

00 


(152) 


Approximations of greater accuracy may be obtained by Including the next 
largest eigenvalues and their corresponding eigenfunctions. 


4.6 The Static Vibration of a Plate - The Boundary 
Value Problem and Green's Function 

The static vibrations of a plate may be modeled by the partial 

differential equation 

V^u = f(P) , PeO (153) 

4 9 2 

where V = V“^(V ) Is the biharmonic operator, and again certain physical 
constants have been included in the forcing function f for simplicity. 


50 


We wish again to choose the boundary conditions such that the boundary 

4 

value problem is self-adjoint. Green's theorem for the operator V takes 
the form 

[ - uv\)dP • f Iv It lu + (^^v)(^)-(V «)^] . (154) 

Jfl 

The problem of boundary conditions for plate vibrations is much more difficult 
than for the membrane. A useful discussion of boundary conditions is contained 
in [4] . 


The Simply Supported Rectangular Plate 

Consider a uniform rectangular plate on the domain 8 “ {(x,y) | 0 ^ •» 

0 < y < b). The boundary conditions for a siiiq>ly supported edge are 


" • 0 “"I-T + RsS’ 0 
3n 


(155) 


vhere n is the normal vector to the edge and R is the radius of curvature. 

For a straight edge R ■ •. Furthermore, since u is constant along the edge, 

■ 0. Tims 
3s 


„2 3^u ^ 1 . 3^u _ 3^u 

V u ■ — s- + B a;r + — T • — T 

an^ ^ as^ an^ 


(156) 


and the boundary conditions for a simply supported straight edge are 
2 

u « V u « 0 . (157) 

Clearly the conditions (157) make the right side of (154) equal to zero. 

4 

Thus V is self-adjoint for the simply supported rectangle. 


The Green’s Function 


The Green’s function for the simply supported rectangle should satisfy 

g(x,y,4,n) - «(x-0 6(y-n) (158) 


g M v^g « 0 


on X » 0 , a and y » 0, b . 


(159) • 
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The free space solution Y<pIq) which satisfies » <(p|q) Is 

y(p1q • 

where r represents the distance This Is proved In Appendix A. Thus, 

the Green's function 

g(P|Q) - Y(p|Q) + g(P|Q) (161) 

^ # 

where the function g satisfies V^g(p|q) • 0, plus boundary conditions such 
that g satisfies (157) . 

The function g in (161) Is no longer necessarily harmonic, as was the 
function in (129). It muse In addition satisfy two sets of boundary condl- 
tions. Thus It is much more difficult to determine the exact function g(P|Q) 
for a given set of boundary conditions. The Green's functions and solutions 
to the shape control and estimation problems will therefore be exhibited in 
terms of eigenfunction expansions. 

4 

4.7 Control Problem for the Operator V 
On the rectangle O^x^a, O^y^b, we desire to determine the set of 
forces {f^} at positions?^ - (x^.y^, 1 ^ 1 ^ m, which together with the 
solution u(x,y) of 
L ® 

V u ■ I 6(x-x.) 6(y-y.) (162) 

i«l 

2 

u » V u on X ■ 0, a and y ■ 0, b (163) 

minimize the performance criterion 

j(F,u) - ^ I '*'t[ f ('*'(*.y)“u(>‘»y))^ (i64) 

1-1 ^ ^ Jo Jo 

overall admissible pairs (F,u). 

The optimal shape for the problem (162-164) is given by 

* “a 

u (x.y) “If. g(x,y, X ,y ) 

1-1 ^ ^ ^ 


(165) 
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where the vector of optimal forces F* ■ (f * ... f *) satisfies 

X n 


(R + A)F - B 

The m X m matrices and A have coordinates 


R^j - r^ «(1-J) 


1 

ra rb 


‘ij 


f f gU.y.x^.yj) gU.y.Jt^iyJ dy dx 

Jo ^0 


( 166 ) 

(167) 

(168) 


and the vector B has coordinates 

bf - f [ iKx.y) g(x,y,x^,y^) dy dx. (169) 

* Jo ^o 

Since a complete analytical form for the Green's function is not known, we 
use the eigenfunctions 


. knx . luv 

a b 


♦,,t (x,y) - 


and corresponding eigenvalues 
‘kt • “I'' ♦ 


to represent the solutions (166-169). Thus 

kiTx 
m 


(170) 


(171) 


u (x.y) • I I « 


. kiTx 
sin sin 

a 


1 . iii , ‘’’i 

— .In .In ^ 


2 2 


1-1 k,i-i * ,« (A (i, ] 


(172) 


where the forces satisfy (166), and the coefficients of the matrix A and 


vector B in (167-168) are given by 


I 




/ kiix^ k” 

(sin sin — 

\ a a 


lTI> 


sin -r- 
a b 


illy 


sin 


lity 




^1 kj.i 

\ 

'^kf,*'*'^ “ I I '*'(*»y^ *** 


where 


(173) 


(174) 


( 175 ) 


Approximations are available by taking the first few terms in k and I* 
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4 

4.8 The EetlMtlon Problem for V 

The shape estimation problem for a rectangular plate is given 
the shape observations 

Vi - u(x^.y^) + 1 i i 1« . (176) 

to determine the error function f(x,y) and corresponding shape function 
o(x,y) which satisfy 

V^u(x.y) - f(x,y) , ( 177 ) 

2 

u • V u ■ 0 for X • 0, a and y ■ 0, b 


and minimize the criterion 

J(f,u) • y J (Y^ - u (x. ,y.))^r*^+ j [ [ f^(x,y) dy dx . (178) 

i ^o ^o 

The necessary condition for an optimal solution is that the vector 

X- (u*(x, ,yj . . . u*(*„.y«)) satisfy 
XX ID in 

(1 + A X) - a R"^ Y ( 179 ) 

where u*(x^,y^) is the optimal shape estimate at the point (x^.y^). 
the matrices A and R are defined by (167-168), and Y “ (Y^^ ... V^) . The 
optimal noise estimates are 


^ 1 *“ ^1 “ “ • 


1 ^ i ^ m , 


10 

f*(x.y) ■ I r,"^ G(x,y,x ,y ) 

i-1 1 1 

and the optimal shape estimate is 

a ® 1 r* 

u (x.y) - I Ir " cj gCx.y.C.n) g(x.,y.,4,n) dnd^. 

i»l ^ ^ ‘'o •'o 


(180) 

(181) 

(182) 
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To compute the vector X in (179) we use (173) for the elements of the matrix 
A. Then 

m • . . Xitx. , firy. 

f'*(x,y) •I I -wj — lain sin ^ sin sin -j—] . (183) 

1-1 k,t-l a a D D 

Finally, applying the expansion (72) to the optimal shape estimate (182), 

m • , , . kxx. - ley. 

u*(x,y) • I I 4^* 2 ®^® ®^" "b — ^ 

1-1 k,l-l abXj^j^ 

Again, approximations are obtained by taking the first few terms In k and 1. 

4.9 Conclusions 

Green's function techniques have been applied to the solution of shape 
control and estimation problems which have associated boundary value problems 
involving partial differential equations. In a manner analogous to those 
Involving ordinary differential equations. In the case of a multidimensional 
domain, however, precise knowledge of the analytical form of the Green's 
function is usually not available. Solutions may be expressed In terms of 
eigenfunction expansions. 

Although this chapter deals with systems which do not have rigid body 
modes, the techniques and solutions bear such a resemblance to those of the 
one dimensional case that an extension to systems with rigid body modes 


follows readily. 
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Chapter 5* Static Shape Control for Nultidiiaenslonal 
Large Space Structures 

5 « 1 Introduction 

This chapter addresses the problems of static shape control and shape 
determination for multidimensional structures* Chapters 2^ have addressed 
these problems for scalar shape functions, representing displacement in one 
direction, defined on one or multidimensional domains* However, large space 
structures are modeled as multidimensional states, representing translations 
and/or rotations in three dimensional space* 

We again use an integral operator approach based on assumptions of 
linear self-adjoint dynamics and boundary conditions* As might be expected, 
algorithms which are similar in appearance arise* 

However, there are important differences in interpretation and proce- 
dure. These include matrix, rather than scalar, differential and integral 
operators, controls and observations applied to only a part of the state, 
and the necessity for using approximate eigenfunctioiis provided by experi- 
mental or numerical methods, since the exact operators and corresponding 
eigenfunctions are usually not known. The algoritluns derived in this chapter 
will be adapted to the use of modes from a dynamic finite element model, 
and illustrated by simulated results, in Chapter 6. 

P rocedur e 

In section 5.2 we define the multidimensional linear boundary value 
problem fur a large space structure, and discuss the existence of solutions. 
We tlien define Green's functions for a multidimensional boundary value 
problem, both witli and without rigid body modes, and derive solutions to the 
boundary value problem for both cases. 
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In saction 5.3 ye define and solve the shape control problem for e .large 
space structure. discuss examples of the constraints Imposed on the 
control forces by the nresence of rigid body ir^t^des. In section 5. A we define 
and solve the shape determination problem. 

We present eigenfunction expansions for the more general multidimen- 
sional terms in the algorithms, vdiich involve Green's functions, in section 
5.5. A sumnary and conclusions are stated in Section 5,6. 


5.2 The Model and the Green's Function 
Consider a multidimensional system represented by the n dimensional 
state U(P), defined on a simply connected domain G e R^. Suppose the 
system Is governed by linear dynamics 

L U ■ F for F e n (185) 


where L Is an n x n matrix of differential operators. F(P) is an n dimen- 
sional vector function, or distribution, defined on G, which represents forces 
or torques acting on the system. 

Suppose the system satisfies linear boundary conditions 

B^(U) - 0 , 1 1 i i , for P e r (186) 

where r is the boundary of n. We will assume the boundary value problem 
(185-186) is self-adjoint, that is that L* ■ L and 

<LU, V> • <U, LV> (187) 


where U and V are any two admissible functions which satisfy the boundary 
conditions and <U,V^ is the inner product 


<U,V> 


U^(P) V(P) dP . 

h 


We will also require the usual vector inner product 


(188) 


<x,y> » xS’ » v^x . 


(189) 



We will use the norms induced by (188-189) and the weighted semlnoirm 

\\^{\^ » 

X and Y are vectors in the same space and R Is a symmetric square matrix 
of appropriate dimension such that R ^ 0. 

The reasons for the model formulation (185-186) become apparent ^hen 
one considers an LSS (large space structure) ant^>nna• The domain consists 
of the subset of three dimensional space occupied by the undistorted ideal 
shape, a perfect paraboloid. The state might be three dimensional also, 
representing vector displacements of points on the distorted antenna from 
their ideal positions. Boundary conditions represent a pinned antenna 
which may not rotate or translate as a rigid body in any direction, a free- 
free antenna which may rotate or translate along any of the three axes, 
or conditions between these two extremems. 

Other state representations are possible. It may be convenient to 
consider a six dimensional state which represents translations and rotations 
of a i 'Int about the three axes. This is the case if torques are to be as 
control mechanisms, in addition to translational forces, A torque can be 
considered an impulsive force applied to a rotational coordinate of the 
state , 


Solutions of Boundary Value Problems 
We consider under what circumstances solutions to boundary value 
problem (185-186) exist, and what form the solutions take if they do exist. 

We will apply the following alternative theorem for boundary value 
problems : 

Theorem 5,1: Consider the boundary value problem 

LIJ = F, B.(U) = 0 , 1 1 i 1 ’ 


(191) 
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Its corresponding homogeneous problem 

LU = 0 , B^(U) - 0 , 1 1 i 1 . (192) 

and the related homogeneous adjoint problem 

L*V = 0 , B *(V) » 0 , 1 < i < k . (193) 

1 ^ o 

L Is an n X n matrix of linear differential operators, L* Is Its adjoint, 

U and V are vector functions defined on the simply connected domain Q, and 
and are adjoint linear boundary operators defined on F, the boundary 
of 

Then; (a) if the problem (192) has only the trivial solution U = 0, so 

does the problem (193), and (191) has a unique solution. 

(b) if (192) has s Independent solutions U, ,...,U , then (193) has 

X s 

s independent solutions , and (191) has solutions if and only if 

X s 

<V ,F> = I v/(P) F(P) dP = 0, llils. (194) 

•'fl 

If the conditions (194) are satisfied, the general solution of (191) has 
the form 

U(P) = U(P) +l U^(P) (195) 

where U is a particular solution of (191), the c^ are constants and U^, 

1 ^ i ^ s are the solutions of (192). 

For discussions and proof of alternative theorems see [2). 

We have assumed the linear operator L and boundary conditions 
B^ = 0, 1 1 1 are such that the boundary value problem (191) is self- 

adjoint, that is that L* * L and B^* * B^, 1 1 1 1 k. Thus (192) and 
(193) are equivalent for our purposes. 

To observe the form the solutions actually take, we define Green's 
functions for the cases (a) and (b) of Theorem 5.1. 
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Green’s Functions 

We first consider case (a) of Theorem 5.1, that the homogeneous 
boundary value problem has only the solution U = 0. This Is equivalent to 
the physical assumption that the system has no rigid body modes. 

Define the n vector functions Gj(p}q), 1 j to satisfy 

L G^(P|q) - 6(P-Q) 

- 6(P%^) ... 6(P%S (196) 

B^(Gj) « 0 , 1 1 i 1 P ^ 

The unit vector e^ has zeros in all coordinates except the jth, where it 
has the value one. The points P(p ... p ) and Q(q ... q ) in (196) lie 

in S3. 

Gj (P 1 Q) represents the response of the system to a unit impulsive 
force applied to the jth coordinate of the structure at the point Q. 

Define G(p|Q) to be the nxn matrix function with columns , G(p|q) 

is the desired Green's function for the boundary value problem (191). The 
ijth coordinate (pjQ) represents the response of the ith coordinate 
of the state at P to a unit Impulsive force applied to the jth coordinate 
of the state at Q. We may write 

LG(pIq) - f(P-Q) (190) 

B^(G) « 0 , 1 1 1 1 1^0 ’ 

if it is understood that the boundary conditions in (199) are to be applied 
to each column of G Individually, 

The property derived in the next theorem will be useful when writing 
the solution of (191) in terms of the Green's function G, 

Theorem 5.2 Let G(p|Q) be the function defined by (196-197). Then G(p|Q)-G (QjP). 
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Proof ; For the noment we drop the assumption that the boundary value problem 
(191) Is self-adJolnt. Let Gj(p|q) and H^(P|r) be functions defined on Q such 
that 

L Cj(p|q) - ej 6(P-Q) , By(Gj) - 0. 1 iv Ik . 

L* H^(pIR) - e^ 6(P-R) - 0 , 1 iv <,k . 

Since Gj and satisfy adjoint boundary value problems, 

<Gj , L*H^> - <LGj . H^> , 1 1 i J i » • 

Thus 

j Gj^(P|Q) e^ 6(P-R)d? - ej*^ 6(P-Q)H^(P|R)dP . 

Evaluation of the integrals yields 

Gj’^(RIQ) Hj^(QIR) . (200) 

G^j(R|Q) » Hj^CqIr) 

But now we recall that L*«L and Thus Hji(QlR) • Gj^qIr). Substi- 

tution into (200) yields 


G^j(R|Q) - Gj^(q|r). Ill, jin.# 


We now seek the solution to the boundary value problem (191) , assuming 
case (a) of Theorem 5.1. Let U(P) be a solution of (191) . Then 


<U,L G 
where the 
<L 


.> - I U^(P) e 6(P-Q) dP - 
J Jjj J 

is the jth coordinate of U. 

U>-<G.,LU>-f G '*'(P|Q) 

J J 


«j(Q) 

By Green's theorem 
F(P) dP. 


Thus, 

U,(Q) - I g/(p|q) F(P) dP, 1 1 J 1“ . 

3 jjj 3 

If we apply this argument to all coordinates 1 1 j in, we have 
U(Q) - f G^(PIq) F(P) dP - j G(q|p) F(P) dP 


61 


by Theorem 5. 2. A change of variables yields the solution 
U(P) « f G(p1q) F(Q) dQ . 

The Modified Green's Function 


( 201 ) 


We now consider case (b) of Theorem 5.1. We assume the boundary value 
problem (191) has s Independent solutions , . . . , which we assume 
are orthonormal. If they are not, a Gram-Schmidt orthogonalizatlon process 
can be applied to generate an orthonormal set. 


Define the following vector functions: 

LG,(p|q) - (6(P-Q) - I V (P) v '*'(Q)l e 
J i«l 


B^(Gj) « 0 , 


1 ^ i <_ k , 


( 202 ) 

(203) 


where e^ is the jth column of the nxn identity matrix. 

Note that the right hand side of (202) has zero components in the 
space spanned by the functions that is that its inner product with 

these functions is zero. Thus by Theorem 3.1 a solution, in the distributional 
sense, to these problems exists. 

The solutions G^ which satisfy (202) are not uniquely determined, since 
the addition of any linear combination of solutions to the homogeneous 
problem (192) yields another solution. Thus we are free to impose 
another condition. We require that 


<G. , V^> - 0 , 


l<_l^s, lUin 


(20A) 


Mathematically this means we seek the solutions to (202-203) of minimum norm, 
those which lie in the orthogonal complement of the nullspace of the 
operation L - the space spanned by the solutions . Thus the functions 
Gj have no components in the direction of the rigid body modes. Solutions 
of (202-204) are unique . 
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Let G(p|q) be the n x n matrix function whose columns are the functions 
Gj which satisfy (202-204), that is G • (G^^ | . . . | G^j , Then G satisfies 


8 X 


LG - 6(P-Q) 

- 1 V (P)(V. (Q)) 
1 ‘ 

(205) 

B^(G) - 0 , 

1 ^ i ^ k , 

<G, V^> - 0 , 

1 ^ i ^ s . 



G(p|q) is called the modified Green *s function for the system (191), 
assuming case (b) of Theorem 5.1. The property derived in Theorem 5.2 may 
also be shown to be true for modified Green's functions. 

We seek a solution U to the boundary value problem (191) for case 
(b) of Theorem 5.1. We assume 


<F, V^> » 0 , 1 1 i 1 8 » 


(206) 


since without these conditions a solution does not exist. We will apply 
Green's theorem to the inner product <u, LG>. From (205) 

<U, LG> « f U^(P)II 6(P-Q) - I V (P) vJ(Q)) dP 

” i«l ^ ^ 

s . 

= u’^(Q) - i ( U^^(P) V (P) dP) V (Q)’^ . 

1 

But 

<U,V > - f u'^(P) V (P) dP 
* some constant c^^. 

Thus 

<U, LG> » u'^(Q) - I V^(Q)’^ (207) 

On the other hand, because the boundary value problem is self-adjoint 

<U, LG> « <LU, G> « I F^(P) G(p1q) dP . 

'ft 


( 208 ) 
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Equating (207) & (208) and taking the tranaposa, we have 
U(Q) - X c. V.(Q) + f C^(P|Q) F(P) dP . 

We apply Theorem 5.2 and a change of variables: 

U(P) - f C(pjQ) F(Q) dQ + X c. V (P) . (209) 

As one might expect from Theorem 5.1, the solution includes an arbitrary 
linear combination of rigid body modes, or solutions to the homogeneous 
problem. 

Remark 5.1 Naturally if a force is applied which does not satisfy the 
constraints (206)^ the system will still respond, but the boundary conditions 
will be violated. The conditions (206) usually translate physically into 
conditions that net forces or torques in one or more directions must be 
zero. 

Without loss of generality, we can define a coordinate system 
with respect to the space vehicle itself. In the case of the antenna we 
define the xy plane tangent to the hub of the antenna and the z axis along 
the axis of the paraboloid. We may fix the x axis along a particular rib. 
With the coordinate system so defined, we may ignore the rigid body modes, 
since rotations and translations of the antenna as a rigid body occur with 
respect to another coordinate system. 

We can then consider the solution of (191) to be 

U(P) - I G(Pl0) F(Q) dQ (210) 

‘'ft 

where G(p|q) Is the modified Green's function which satisfies (202-204). 

5.3 The Shape Control Problem 

Static shape control forces may be applied to some or all of the 
coordinates of the multidimensional state. Thus we define the following 
control problem: 
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Let <|)(P) be the desired shape of the space structure. Detemlne 
the set of control vectors of (predetenslned) dimension n(l) , 1 1 n* 

such that the resulting shape U(P), which satisfies the dynamics 


LU - I C F^ «(P-P^) 

and boundary conditions 

B^(U) - 0 . 1 1 i 1 

most closely approximates the desired shape 4^ on Q. The measure of best 
approximation is that the set (Fj^*,... ,F^*, U*) minimize the performance 
criterion 

J - 7 1 I iFjl I, ^ * Hn ' ■ “‘’'’I 'm(p) 

over all possible sets which satisfy (211-212) . 

The constant n x n (i) matrices distribute the control vector F^^ 
over the coordinates of the state U at P^. 5(P-P^) is the dlrac delta 

function for the multidimensional point P^. 

The n(i) x n(i) matrices are symmetric and 0. 

U(P) is a piecewise continuous symmetric positive definite matrix 
defined on fl. 

We first assume the homogeneous system 

LU - 0 » Bj^(U) - 0 , 1 1 i 1 <214) 

has only the solution U = 0. 

We apply the solution derived In section 5.2 to the boundary value 
problem (211-212) ; 


U(P) - f C(P|Q)lf C. F 6(Q-P.)1 dQ 
iQ 1 ^ ^ 


m 

- I 

!•! 


C(P|P^) 


(215) 
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where G(p|q) is the appropriate Green's function. We substitute (215) into 
the criterion (213) , ^ich becomes a functional depending solely on the 


discrete unknowns F^. 


(216) 


J - 1 1 I IfJ 1 I ;♦(?) - g(pIp^) dP . 

We seek the minimum of J with respect to the constant vectors F : 

Jfl " Jl G(p|Pj^)CJ^l^W(P) l-G(P|Pj)Cj] dP (217) 


0, lljlm 


Thus 


R.F. + I c/ ([ G(P |P)G(P|P.) dP) C F 
J J i-1 J 

- j G(P^|P)W(P) <KP) dP for 1 1 j . 

xn ^ 

Let N - I n(i) . 
i-1 

Let R be the block diagonal square matrix with diagonal blocks 
Let A be the N X H matrix of n(i) by n(j) blocks Aj^^, where 
Let D be the N dimensional vector 


( 218 ) 


( 219 ) 


( 220 ) 


(221) 


where 


Dj - C 


i L 


W(P) *(P) dP . 
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Let F be the N dimensional vector of unknown control forces: 

<^1^ • (222) 
Then the vector F* of the optimal control forces satlfies 

(R + A) F* • D . (223) 

Once the vector F* has been determined, the optimal shape U*(P) la 
given by 

U*(P) - I G(P|P.) C.F * . <224) 

i«l ^ ^ ^ 

The Shape Control Problem for Systems with Rigid Body Modes 

We now assume that the homogeneous boundary value problem has s 

solutions V, (P) . ,V (P) . We let G(pjQ) denote the modified Green’s 
X s 

function which satisfies (205). 

In order for a solution (211,212) to exist, the right hand side of 
(211) must satisfy the additional set of constraints (206) . That is 
m 

<V I C F 6(P-P.)> » 0 , llils, 

j-1 ^ J 

which by definition is the set of conditions 

I v/(P) I C F 6(P-P.) dP»0, llils. 

J ^ J J J 

Evaluating the integral yields 

! vJ(P ) C F • 0 . llils- (225) 

j-1 ^ ^ 

The shape control problem is now to find the set of forces {F^> and shape 
U(P) which satisfy (211-212) (225), and minimize the criterion (213) ever 
all possible sets. 
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Exanples.of Constralntg 

Example 5.1: A three dimensional static with three rigid body modes. 

Suppose the state U(P) is three-dimensional, representing the displacement 
vector of the actual position of the structure corresponding to the point P 
from its ideal position, and the antenna has three rigid body modes representing 
translations along any of the three axes. An orthogonal basis for the space 
spanned by these rigid body modes is 



Note that if U(P) is a three-dimensional state then 

U(P) + f c.V 
1-1 ^ ^ 

does represent a translation of that state. 

The constraints (225) become 
m . 

I (C F - 0 <226) 

j«l 3 J 

where coordinate of This is equivalent to the 

condition that the net force applied in any direction of the state U over 
all the points P^ is ^ero , If the sum of the forces in any direction is 
zero, no net acceleration is applied to the structure as a idiole, which is in 
keeping with the free boundary conditions. 

Example 5 .2; A six-dimensional state. 

If torques are to be applied as part of the control scheme it may be 
convenient to consider a six-dimensional state, the first three components 
of which represent displacements as before, and the second three components 
of which represent rotations-. 
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A torque is en Inpulsive force applied to e rotetionel coordinate. 

Suppose that the systea has six rigid body stodes, representing constant 

translations or rotations from an ideal position. A basis for the space of 

rigid body modes is (1 0 0 0 0 0)^, (0 1 0 0 0 0)^. (0 0 1 0 0 0)^. 

(0 0 0 1 0 0)^, (0 0 0 0 1 0)^ and (0 0 0 0 0 1)^. The later three vectors 

represent unit rotations about the three axes. 

The constraints (225) again become 
m . 

I (C F - 0 , 1 < i < 6 . (227) 

j.l J J ~ 

These constraints represent the fact that the net sum of forces or torques 
applied to any coordinate of the state must be zero, a requirement which guarantee) 
zero translational or rotational acceleration applied to the state. 


Example 5.3: A three-dimensional state with six rigid body siodes. 

Suppose for computational convenience we wish to consider a three- 
dimensional state, but the vehicle is allowed to both rotate and translate 
along three axes as a rigid body. One basis for the six rigid body modes 
is 


- (10 0) V2 * (0 1 0) V3 - (0 0 1) 


where 


Vj(P) - T2P 

V^(P) - T3P 



0 

0 


1 COS0 

0 

sine 

CO80 -sine ' 


T2- 0 

1 

0 

sinO cosd j 


\ -sine 

0 

cos6 


and 



T^. T2 and T3 represent rotations by an angle 6 about the x. y, and z axes 
respectively . 
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The first three constraints yield the same conditions as In example 5.1. 


The last three constraints yield 

I P/ T/ C.F, - 0 
j-1 


T « T 

j 1 ''J'J 


1 <1 <3 . 
For the rigid body mode V^(P) ■ T^P this la 


(228) 


I CO80 - alneX (Pj^ nine + p^^ co80)l C^F^ - 0. (229) 

1 


This expression Is the requirement that the sum of the forces applied 
times the displacements at the points vhere the forces are applied must 
be zero. But this Is just one of the constraints which resulted from 
rotational rigid body modes In example 5.2: the sum of the torques must 
be zero. 

It Is easily seen that the condition that the sum of the torques be 
zero for each coordinate Is satisfied if the constraints (228) are satisfied. 
Thus, the constraints for six rigid body modes are the same, however the 
state vector is defined. 

The procedure tor finding the set of optimal control vectors for systems 
with rigid body modes Is as follows: 

1) Substitute the solution (215) Into the criterion (213) 
li) Solve the constraints (225) for some of the control vectors In terms 
of the others. 

iii) Substitute tho expressions derived in (11) into the criterion J, 
which now becomes a function of fewer control vectors, 
iv) Minimize J with respect to this smaller set of control vectors. 

The minimization process will result In a system of linear 
equations which, when solved, yield the identity of the optimal 
set of these vectors. The other control vectors may then be 
determined from (11) . 
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The plnned»fr«e b««a In section 3.3 trss sn sxsmpls of this procedure in the 
esse of s one-diaensionel state. 

5.4 The Shape Deteraination Problea 

The desired shape f in the control problea in the last section vill be 
based on the difference bet«reen the catlaated shape and the ideal parabolic 
shape. The estiaated shape aust be coaputed froa observations of soae or 
all of the components of the state, taken at a number of predeterained 
points along the structure. 

Thus we seek to deteraine the estimates of the noise vector F(P) and 
shape function U(P), defined on Q, based on the observations 

Yf - U(P^) + 2^ . 1 1 i 1 a . (230) 

which minimise the performance criterion 

j • X Ji • 1^1 ' ^1 ^ L ’ * *'^^^**^ **** 

over all admissible sets (U,F) which satisfy 

LU - F, ? c n and B^(U) - 0 PcT. (232) 

Th. ' jtant matrices are n(i) x n, the n(i) dimensional vectors 
repr..8vint nols? cr inaccuracies in the observations Y^, W(P) is a 
continuous poi-itive definite matrix on 0, and are n( x n(i) constant 
positive definite matrices. 

We will assume the boundary value problem (232) has no rigid body 
modes. The estimation algorithm for systems with rigid body siodes is the 
sane, with the exception of the fact that the rigid body modes themselves 
cannot be estimated. The derivation of this fact follows as in section 2.5. 
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We will evuiaute the solution (210) of the boundary value problem (232) 
at the points 1 ^ 1 ^ m, and substitute into the criterion J. 

u(P.) = ( G(P Iq) F(Q) dQ . (233) 

•J - tI ll^i - C f G(P IQ) F(Q) dQllJ.l+i f Nf(Q) 1|2 d(. . (234) 


The criterion is now solely a function of the continuous unknown vector 
function F(Q)« To minimize J with respect to F we find the Frechet derivative 
3J(F,H)» where H is any admissible variation^ and set it equal to zero* 

3J(F,H) = i 'Y - C f G(P 1q) F(Q) dQ]^ rA-C. f G(P Jq) H(Q) dQ] 

1 •'ft ^ 

+ f F(Q)’^ W~^(Q) H(Q) dQ * 0 . 


If we transpose the equation, factor out H and recall (233), we have 

[ r(Q) w'Aq) IF(Q) + I G(q[p,) c V'Ac u(P ) - Y )] = 0 
JQ 1=1 

Since thio wst be true for all ;nr'ssible variations H, we have 

F(Q) = W(P) I G(o! ' C.^rT^ (Y. - C U(P )) . (235) 

1=1 ■ 11111 

Ve still do not know the optimal estimate of F at this point, because the 

estimates C. U(P.) are still not known. We substitute (235) into (233). 

X i 

U(P ) = f G(P jQ) W(Q) [ f G(QlP^)cARi (Y^ - G(P^)) • (236) 

^ ^ i-1 


Then we have , for 1 j m » 

U(P.) + I (f G(p Iq) W(Q) G(q|p ) dQ) cJr, C U(P ) 

J i=l ■'ft ^ X A - 

= ! f G(P Iq) Q(Q) G(q1p ) dQ) c\ Y . (237) 

i=l •'ft J 

We will solve this set of m matrix equations for the vectors C. U(P^), 

J J m 

m 

i - i < m. MuZtiplv both sides on the left by C , . Again define N = T ti(i). 

' J i*l 


(Recall that n(i) is the dimension of Yj,.) 



72 


Let A be the matrix of m blocks by m blocks, where the Jth block 

is an n(J) by n(i) matrix. Thus A is M x M. 

Let r”^ be the N x N block diagonal matrix with blocks 


Let U*(Pj^) be the optimal estimate of the shape function U at and let 
C be the N dimensional vector formed by "stacking" the n(l) dimensional 
vectors (240) 


Let Y be the N dimensional vector 

(y ^ y 

'*1 •** m ' * 

Then the vector 0 satisfies the system of linear equations 
(Ijj + AR”^) D « a R“^ Y . 


(241) 

(242) 


Once the vector C is known ^ the optimal estimate F* of the noise vector 
F is given by 

5 T 1 

F*(P) « W(P) G(pjp^) R^ (Y^ - C^u*(p^)) . (243) 

The optimal shape estimate U*(P) is then given by 

ID f 

U*(P) * (J^ G(?1q) W(Q) C(QlP^)dQ) (244) 


5 . 5 Appr oxlmat ions 

In this section approximations will be presented, which Involve 
eigenfunctions corresponding to the static boundary value problems (211-'2l2) 
(232) which parallel those in section 1.6. 

However , most finite element models for large space structures are 
dynamic, rather than time-invariant. Therefore, in the next chapter, approxi 
matlons will be developed for the use of eigenfunctions from the dynamic 
model corresponding to 0.85-186) , 



I 

I 

I 
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It was demonstrated In Theorem 5.2, section 5.2, that G(p{q) was 
symmetric. We will also assume that G(p|q) Is a Hllbert-Schmldt kernel. 


that is that 


f 11 g(p1q) 1|2 dPdQ < « . 
J o 


(245) 


Let K be the integral operator with the Green's function as kernels Then 


for F(P) In the domain of K, 


Kf = f G(p!q) f(Q) dQ . 
Jo 


(246) 


Let ^ ^2 — — •*. be the non-zero eigenvalues of K, and ^ 2 ^ 

be the associated eigenfunctions, such that 


(247) 


The non-zerc eigeu alues of K are the inverses of the non-zero eigen- 
values of L. end the eigenfunctions are also the corresponding eigen- 

functions of L* 


We will assume the eigenfunctions {((>^} have been normalized with 
respect to the inner product (188). 

From integral operator theory we have the following expansion for the 


Green’s function: 


G(i'lQ) ^ I u. 4>AP) 4>/(Q) 

1-1 


(248) 


If we assume, es ir Chapter 2, that W is the identity (matrix) on we 


have the following expansions: 


f G(p|q) G(q1r) dQ = I <{)^(P) 


(249) 


f G(p1q) i;;(Q) dQ = I <^^{P) . 


(250) 
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These expressions are generalizations to the multidimensional case 
of the approximations offered in Chapter 2. 

If W is not the identity matrix, (249-250) become 

f G(P|q) W(Q) G(q1r) dQ - X ^ y U ♦ (P) 4'/(Q) <♦. , *.> (251) 

Jjj 

and 

f G(P|q) W(Q) «KQ) dQ « [ y ♦ (P) «|>. . (252) 

J n i-1 

:>*6 Summary and Conclusions 

Procedures for static shape control and determination of multidimen- 
sional large space structures were derived in this chapter, under the assump- 
tions that the structures were continuous, governed by linear self-adjoint 
boundary value problemt^, and that the control forces are applied and observa- 
tions taken at a number of predetermined points along the structure. 

Approximate optimal control functions and shape estimates, in terms of 
eigenfunctions corresponding to the static model, were presented. 

As one would expect, the problem formulations and solutions for multi- 
dimensional states bear a strong resemblance to those for scalar state 
formulations derived in Chapter 2. This is due to the commonality among linear 
self-adjoint systems. 

However, there are significant differences in interpretation and 
procedure. The differential and integral operators become matrix operators 
rather than scalar. Observations and control forces may now be applied to 
parts of the state, on to linear combinations of state components, rather 
than to all of the state. The additional constraints imposed in the case 
of rigid body modes must be interpreted and handled with ;nore care. 

Finally, it is now nearly Impossible to know the differential and 
integral operators, or their eigenfunctions, with analytical precision. 
Approximations must be supplied using eigenfunctions computed experimentally 


or by a numerical method such as the finite element method. 
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Chapter 6. Finite Element Models: 

A Large Space Antenna 

6.1 Introduction 

In Chapter 5 static shape determination and shape control algorithms 
were derived for a multidimensional model defined on a multidimensional 
domain, the situation most likely to correspond to large space structures. 

It was assumed the structural models satisfied static self-adjoint linear 
boundary value problems ol the form 

L U(P) - F(P) . U(P) - 0 , 1 1 i , (253) 

where U(P) represents an n dimensional state vector of displacements at 
the point PeJi^Lisannxn matrix of differential operators and 
1 ^ i are linear boundary operators defined on the boundary F of 

Terms in the solution algorithms for the stati": shape estimation and 
control problems involved the Green*s function^ or impulse coefficient, 
of the associated boundary value problem* Since it is highly unlikely 
that the precise Green’s function for such a problem is known, approximations 
to these terms by means of expansions involving the eigenvalues and eigen- 
functions which satisfy the corresponding eigenvalue problem 


L 




®i^j 


0 . 


1 < 




(254) 


were presented. 

However, it is likely that the most convenient eigenfunctions will be 
those supplied by a finite element model, which approximate those for the 
dynamical boundary value problem 
2 

M(P) ^ + L U(P,t) » F(P,t), P c ft, (255) 

3t^ 


1^1^ Rq . Per 


U(P,t) - 0 


» 


(256) 
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associated with the static problem (253) . These eigenfunctions satisfy 


L ij(P) - Xj M(P) }j(P) - 0. Pen, 

(257) 


(258) 


and arc orthononnal with respect to the norm induced by the weighted 
inner product 

<U,V>„ - f ’/(P) M(P) U(P) dP (259) 

“ J(2 

rather than the usual inner product 

<U,V> » I U^(P) V(P) dP . (260) 

This chapter investigates the modifications necessary for the use of the 
eigenfunctions {<Pj} which satisfy (257-258), rather than those for the static 
problem. 

The finite element method is outlined in section 6.2. In section 6.3 eigen- 
function approximations are derived for terms which involved the static Green’s 
function, using eigenfunctions for the dynamic problem. In comparison, 
we solve the discrete static control problem in section 6.4 in order to 
demonstrate the remarkable consistency between the discrete and continuous 
solutions . 

Finally in section 6.5 we present specific examples of algorithms for 
multidimensional shape determination and control, which are illustrated 
by simulations using an available finite element model of a large space 
antenna. Tables . nJ plots of results are included at the end of the 
chapter . 

For convenience, only the case that there are no rigid body modes, 
or non-trivial soltitions of the unforced (homogeneous) boundary value problem, 
will be considered. T’ :? extension of these results to the case of system 
does have rigid body modes is obvious. 
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6.2 The Finite Element Model 

The finite element method Is a modification of the Raylelgh-Ritz 
procedure for solving self-adJolnt boundary value problems. The Raylelgh- 
Rltz method will be described briefly first. It Is based on two principles: 

1) The unique solution of the self-adjoint boundary value problem %fhlch 

governs a system Is equivalent to the unique function In a certain 
class which minimizes an integral, or functional, which usually 
represents the energy of the system. 

Examples of such equivalences are the following: 

Example 6.1: The solution of the system of linear aquations A x* « b, where 

A Is a symmetric matrix, is equivalent to the unique vector x* which 
minimizes the functional J(x) ■ <Ax,x> - <b,x>. 

This equivalence is equally applicable If A Is a self-adjoint 
linear operator. {10] 

Example 6.2: The function y e C (0,1] is the unique solution of the boundary 

value problem 

- ^ y * 0 <261) 

y(0) - y(l) - 0 (262) 

2 

if and only if y is the unique function in (0,1) which minimizes the 
integral 

J(u) - I {pU)(u’(x)I^ + q(x)(u(x)]^ - 2f(x) u(x)} dx . (263) 

1 0 

(Ref. 111]). 

2) The second principle of the Raylelgh-Ritz method is that the functional 

J is not minimized over all appropriate functions (in example 6.2. for 

2 

example, J is minimized over those functions 1*^ (0,1) which satisfy 

(262)). It is minimized ever a smaller set consisting of linear 
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coabinatlons of certain basis functions referred to as 

coordinate functions, which are defined on the region and satisfy 
the boundary conditions. 

Thus, the s<^lution of the linear boundary value problem becomes the 

question of determining the set of constants such that the function 

n 

£ . ^ c f minimizes J overall such sets, a finite dimensional problem. 

i-1 ^ ^ 

In effect we are finding the best approximation of the solution to the 

original problem in terms of the functions The trick in the Rayleigh- 

Ritz method is to find a sequence of suitable functions such that 

n 

as n goes to Infinity the functions f ■ c converge t.’ 'he solution 

“ i-1 ^ 

of the boundary value problem. Frequently used sets are piacew^s^ 

linear polynomials and cubic splines. 

The finite element method is a modification of the Rayleigh-Ritz 
method for more complicated structures, which cannot be described accurately 
by as simple an equation as (261). The- domaiis of the strsidur^ is divided 
into smaller regions, or elements, which are Ir.crj’covtaected at a discrete 
number of nodal points. 

The displacements of the . ructure at the nodal points form the unknown 
constants. The displacemen • at one node represent translations, rotations 
or higher order terms in one or several dimensions. Within an element, 
a set of displacement functions is chosen to define displacements between 
the nodal points in terms of the displacements at them. These functions 
correspond to the coordinate functions of the Raylelgh-Rltz method. 

A state vector X representing the displacements at all the nodal 
points is formed. The usual order in the vector is that all displacements 
for the first node are first, followed by all displacements for the second 
node, and so on. 
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The object of the finite element method is to determine the state vector, 
or displacements, which will yield the closest approximation to the actual 
displacement pattern of the structure. The derivation of the equation that 
this vector satisfies is an application of the following principle which 
Is analogous to the first principle of the Raylelgh-Rltz method. 


Hamilton’s Principle; Let L » T-V be the Lagranglan of a system, where T 
is the total kinetic energy and V Is the potential energy. Then the actual 
p^th of the system in time, X(t), renders the Integral 


f 2 

L(X,X,t) dt 

^1 

stationary with respect to all possible neighboring paths the system may 
take between times tj^ and t 2 * Therefore the Frechet differential 


3J(X,H) 



L(X + oH, X + aH, t)dt 


o«0 


0 


for all admissible variations U. This Is a classical problem in the 
Calculus of Variations, which leads to the Euler -Lagrange equations for the 
system: 

L^(X,X,t) - ^ L. (X,X,t) - 0 (264) 

(Ref. 13] , p. 181). 

For dynamic finite element models 

T » I X^Mi and V - i x'^KX . (265) 

M and K are square symmetric matrices and M is positive definite. 

The mass matrix M arises out of an analysis of the inertial forces 
acting at the nodes. The coefficients of M are referred to as mass 
influence coefficients, which relate the accelerations at the nodes to the 
resulting inertial forces. is the force at coordinate i due 

to a unit acceleration at coordinate j. The total inertial forces acting 
on the system may be expressed in vector form by F^ » MX. 
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The stiffness matrix K arises out of an analysis of the elastic 
force relationships at the nodes. The stiffness influence coefficient 
represents the force at the coordinate i due to a unit displacement 
of coordinate J . In vector form the elastic forces acting on the system 
X may be written ■ K X. The stiffness matrix K in the discrete systnm 
corresponds to the linear operator L in the continuous systems (253) or (255) . 

The coefficients of M and K are computed by integrations over each 
element using the coordinate functions. 

If the Euler-Lagrange equations (264) are evaluated for the finite 
element model the following equations result: 

for a conservative system: MX + KX = 0 (256) 

and, if a vector of nonconservative (outside) forces F(t) is 
acting on the system: MX + KX » F(t) , (267) 

In a static system X » 0, which yields a system of linear equations 
as a necessary condition for the state X: 

KX = F . (268) 


The final step in the finite element method is to solve (267) or (268) 
for X, the vector of nodal displacements, given a known force vector F. 

The system (267) Is self-adjoint if and only if the weighted inner 
produce 

<X,Y>j^ - X^MY - Y^MX (269) 

is used. Consequently there exists a corplete set of eigenvectors (modes) 
{ 41 ^}, 1^1 where is the dimension of the state X, and corresponding 

eigenvalues (A^), such that 

M « K , 1 1 i 1 Np* (270 

* 2 

The eigenvalue A^ > , where is the frequency corresponding to the mode, 

A 

or eigenvector, 
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Eigenvectors corresponding to different eigenvalues are orthogonal 

under the norm (269). We assume they have been normalized with respect 

to that norm. The solution of (267) is given by 
N 

XCt) - I c.(t) ♦. (271) 

i-1 ^ ^ 

** O A 

where Cj^(t) satisfied 

Thus, given a known vector F of non-conservative forces the solution 

A 

of (269) is expressed In terms of the eigenvectors and frequencies which 
satisfy (270) • These are the modes and frequencies supplied by the finite 
element method, which must be used to approximate the static shape control 
and determination algorithms a 

Because of computational limitations, only a fraction of the total 
number of modes are actually computed « 

The solution of (268) is discussed in section 6.4. 

In summary, the basic steps of the finite element method are as 
follows: 

Summary of the Finite Element Method 
i) The domain is divided into a number of elements, which are inter- 
connected at a discrete number of nodal points. 

li) A state vector X is formed, representing the displacements of which 
knowledge is desired at each node. Displacements within an element are 
expre^^ed in terms of coordinate functions. The unknown constants in the 
displacement functions are the displacements at surrounding nodes. 

The mass matrix M and stiffness matrix K are computed. The state vector 
X there satisfies MX *f KX » F, for dynamical systems, or KX » F, for static 
systems, where F is a vector representing outside forces acting on the 


system. 
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iv) The nodes and frequencies which satisfy M - K 
ere then computed « Solutions to the nodel nay be expressed in terns of 
these modes. 


The Lumped Mass Method 

A simplification of the finite element method, the lumped mass method, 
is frequently used for models of large space structures at JPL. The entire 
mass of Che structure is assumed to be concentrated at the nodal points, 
which are interconnected by massless segments. Thus no coordinate functions 
need be defined. The mass matrix is a diagonal matrix, with identical 
entries for all translations corresponding to the same node, and zeros for 
rotations or higher order terms U2]. 
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6.3 Approximations from the Dynamic Model 

A 

Given the eigenfunctions which satisfy (257->S) and are orthonormal 
with respect to the weighted norm(259\we wish to generate approximations 
to terms in the shape control and determination algorithms. The eigen- 
functions (270) supplied by the finite element method are discrete approxi- 
mations to those of (257-8) . 

If the Green’s function C(?jQ) is not known « we require approximations 
for Che following quantities: 

G(p|Q) (272) 

jj G(p1q) W(Q) G(Q|r) dQ (273) 

jj G(p1q) VI(Q) *(Q) dQ (274) 

where i|»(Q) is a known function and W(Q) is a symmetric positive definite 
matrix. 

We will first assume that we have available the continuous eigen- 
functions for which the finite element method provides approximations. For 
convenience from this point forward we drop the hats ou these eigenfunctions. 


which satisfy the following properties: 

L ^j(P) - Xj M(P) ? t a (275) 

(P) - 0 , 1 1 i i ^ ® ^ 

<♦,, ♦,> - f ♦/(?) M(P) ♦.(P) dP “ • (277) 

J i ® Jj, J ^ 

Properties (275) and (277) easily yield the following property: 

<^j, - I ♦j'*'(P) I t^(P) dP - Xj 6(i-j) . (273) 


The application of the Green's function (198-9) to solve the boundary value 
problem (275-6) yields 


♦j(P)-Xj 1 


(279) 
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If there ere no eigenfunctions corresponding to the eigenvalue X « 0, 
that Is, if the nullspace of the operator !• Is only the zero vector, the 
functions fonn a complete aet for all functions in « suitable class 
which satisfy the boundary conditions. 

If there are eigenfunctions corresponding to X • 0, the modified Green's 
Function defined in Chanter 5 has no component in the nullsnace which is 


spanned by these functions. Therefore in either case the column vector 
Gj (P Iq) can be expanded in terms of the eigenfunctions corresponding to 
non-zero eigenvalues: 

(280) 


c.j(p|q) • I ♦j^(P) Yj^(Q) 


where continuous scalar functions defined on n* If we define 

Yj(Q) “ (Yjj^CQ) ••• 


G(P|Q) « I ♦.(P) Y.CQ) 

i ^ ^ 


(281) 


In order to determine y j * we multiply both sides of (281) on the left by 
4i^^(P) M(P) and integrate over fl. 


dP . 


f v/(P) M(P) G(P|Q) dP - J f 4/(?) M(P) 4.(P) y.(Q) 
h ^ j ■'G j j 

If we apply the orthogonality relationship (277) : 

I 4j(P) M(P) G(P|Q) dP • y^(Q) . (282) 

J T 

If (202) is compared with (279) it is observed that Yj^(Q) * “ ^i • 
and 


G(P|Q) - if" • 

1 Ai 1 


(283) 


where the sum is over the non-zero eigenvalues and elge.tfu tctions of the 
system (255-6). 
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We use the expression (283) to find expansions for (273-4) : 

f G(p[q) W(Q) ♦(Q) dQ » X ^ I ♦.(?) ♦/(Q) W(Q) *(Q) dQ 

» I ^ ♦^(P) j ♦i'^(Q) W(Q) <KQ) dQ 

- I ^4i(P) .<!»>„ . (284) 

Finally we evaluate expression (273) : 
f g(p1q) W(Q) G(Q|R) dQ 

* f (I f- ♦/(R)) dQ 
Joi^i^ X j'^iJ J 

■ H ♦!<'■> ♦/<“> ‘♦l- *j’w 

In the event that the matrix W(P) Is chosen to be the mass matrix M(P), 
the relation (285) becomes 

f G(p|Q) M(Q) G(q1r) dQ » X (^.(P) 4./(Q) . (286) 

•’Q ^ V 

The expressions for (272-4) in terms of eigenfunctions for the dynamic 
problem are very similar to those in terms of eigenfunctions for the 
static problem. The major differenr.^ is the loss of orthogonality with 
respect to an unweighted inner product. 
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Section 6.4 The Discrete Control Problem 

It is satisfying t although not unexpected, to note the resemblance 
bettreen the solutions of the shape control and determination problems for 
continuous and discrete models. For exsmple, the discrete control probl^ u 
analogous to the problem (211-212) is as follows: 

Let X be an dimensional state vector representing the displacements 
for a sequence of nodal points along a structure. Let Y represent the vector 
of desired displacements. Suppose m scalar forces Fj are to be applied 
to coordinates r(J) of the vector X in order to achieve the desired 
"shape" Y. Then the control problem is to determine the control vector 
X which is the solution of 

KX - CF (287) 

and minimizes the criterion 

J • ^ 11^11^ + I 1|X-Y||^ (288) 

over all pairs (F,X) which satisfy (287). 

C is an X m matrix with entries ■ 6(i - r(j)). R is a symmetric 

constant m x m matrix, and M is the mass matrix of the corresponding 
dynamical model. Since we are considering systems without rigid body modes, 
there are no nontrivial solutions of KX » 0. Thus K is non-singular, and 
the solution of (287) is given by 

X - CF (289) 

when F is known* 

Finding K ^ is analogous to finding the inverse of the operator L, 
that is to finding the Green's function such that the solution of LU « F 
plus boundary conditions may be expressed as 

U - L~^ F - I G(p1q) F(Q) dQ . 

•'n 
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As In the continuous cese, while It is eesy to refer to in theory, 
in practice the ayaten dlnenaion ia on the order of 10^, so it is 
desirable to find a neans of approximating rather than actually 
computing it * 

We substitute (289) into the criterion J : 

J - y (F^RF) + j(k"^CF - Y)^ M(k“^CF - Y) . (290) 

We minimize (290) with respect to the unknown vector F: 

|j • f\ + (K"^CF - Y)^ M K“^ C - 0 . 

This equation results in the following necessary condition for F: 

(R + c\“^1K”^C) F • K“^ M Y . (291) 

Once F is known from this m dimensional system of equations, the optimal 
shape X is given by (289) . 

Since it is awkward to compute we seek eigenfunction expansions 

for it, and the terms K ^ M K ^ and K ^ M Y« We assume we have available 

the eigenfunctions and eigenvalues of the corresponding dynamical system 
•• 

MX + KX » f , which satisfy (270), together with the orthogonality conditions 

<«i,^j>jj • «(i-j) (292) 

and 

’ *i ’ ^i 

Let ♦ be the by matrix | ♦ , . ] ] . Then 

o 

K ♦ - A (294) 

where A is the diagonal matrix with diagonal entries 1 ^ i Thus 

X » (*“^) A 
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and 


a a 

ftp. 


i»i "1 


Furthnnoore 


N 


-1 


• ‘i X7 *1 ♦ J) “<[ ^ *P 


N N 
o o 


I I rr- *i 

iti i J i J “ 


and 


N 

1 . . T 

I —J ♦! ♦i . 
i« .l 

N 


K"^ M Y - [ i- Y> . 

1-1 ^1 1 1 


(295) 


(296) 

(297) 


Note the marked resenblsnce between the discrete expressions (295-7) 
the analogous expressions (283-4) (286) in the continuous problem. 
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6.5 Applications to a Large Space Antenna 

In this section ve present actual algorithms for the shape determination 
and control of a large space structure, subject to the choice of certain 
constants, which may be varied in simulations. The algorithms are illus- 
trated using eigenfunctions and frequencies provided by a finite element 
model of a large space antenna, which has been developed at JPL. 

The model is constructed by the lumped mass method described at the 
end of section 6.2. It assumes 18 ribs and 882 nodal points locatea on 14 
concentric circular cross-sections of the mesh. The ribs are assumed to 
be very stiff in comparison with the mesh. The hub of the antenna is 
assumed to be firmly fixed to the bus of a more massive spacecraft, so that 
there are no rigid body modes . 

3 

Available data on the model includes the rest coordinates in R , which 
represent the positions of the nodes on the ideal shape U® , the masses 
at each node, and 33 modes and frequencies. 

We will restate the problems and algorithms to incorporate two subtle 
refinements necessary for the application to a large space antenna. 

The first arises from the fact that the mode shapes and Green's function 
represent displacements of the antenna from its ideal shape. The actual 
antenna shape is the sum of its ideal, or rest shape U° , a perfect paraboloid, 
and its displacement. Thus the Green's function represents the displacement 
of the antenna from its ideal shape due to a unit Impulsive force at one 
point. 

The second refinement is that shape estimation is accumpllshed first, 
and the resulting shape estimate U* is used as the desired shape in the 
control problem. Once the forces necessary to control the ideal shape to 
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the shape estimate are determined, the negative of those forces will bring 
the estimated shape to the optimal corrected shape. 

After the full algorithms are stated, we state the corresponding approxi- 
mations used in the simulations, which are based on the expansions developed 
In section 6.3. 

The results of the simulations Include tables representing comparisons 
of results for varying choices of control and observation positions, number 
of modes In the approximations, weighting matrices and choices of actual 
distorted shapes. Plots of the first eleven mode shapes, and the actual 
distorted shape, estimated shape and corrected shape for various Initially 
distorted antenna. 

The computer program listing and output for the shape control of a 
large space antenna are found in Appendix C. 

The Shape Estimation Problem 

Consider an n dimensional space structure, the shape U(P) of which 
satisfies the following linear self-adjoint boundary value problem on the 
L dimensional domain n with boundary T: 

LU(P) • F(P), Pen (298) 

Bj U(P) - 0 , 1 1 j Per. <299) 

L is an nxn matrix of linear differential operators, which Is related to 
the stiffness of the structure. , 1 ^ j a*’’* linear homogeneous 

boundary conditions. F is a vector function of unknown disturbances. 

The shape estimation problem is to determine the unknown disturbance 
function F* and shape function U*, based on the m observation vectors 

1 ^ 1 ^ m, 


Vj . U(Pj) + Zj. 


(300) 



s 
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which satisfy (298-9) and minimise the performance criterion (301) over all 
possible pairs (F,U) which satisfy (298-9) . The vectors 2^ represent, noise 
in the observations. 

J Ha ‘ I I 

and V are symmetric positive definite weighting matrices of appropriate 
dimensions. 

The Static Shape Control Problem 

Given the optimal shape estimate U*(P), the shape control problem is 
to determine the set of m control forces applied at the positions P^, 

A 

which together with the resulting shape U(P) which satisfies 

LU(P) • f C F. «(P-P.), 

U(P) - 0. 1 1 J Ik, 

minimizes the criterion 

j • I I I I ^ ll|^ + i i i 1 lw(p) dP 

over all possible sets {U,{F^}) which satisfy (302-3). The matrices are 
positive semidefinite and the matrix V is positive definite. 

The forces F^, 1 ^1 ^m, when applied to the positions P^ of the 
ideal shape U*, will produce the closest approximation to U* obtainable 
by the polntvise application of forces at those positions. Consequently, 
because the system is linear, the application of the negatives of the 
forces F^ to positions P^ on the estimated shape U* will produce the optimal 
shape correction of U* to the desired shape U*. 


P c fi 
Per 


(302) 

(303) 
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The Shape Detemlnatlon Algorithm 
Assume the positions observations and their directions 
determined by are known. Choose the weighting matrices and W In 
the criterion (301) . Then 


i) ...Compute the block matrices 1 1 given by 

Aij • cj (j^ G(P^jp) W(P) G(P|pj) dP) C^ 
where G(p|q) is the associated Green's function for the system. 


(305) 


11 ) Form the matrix A whose block coordinates are A^^ and the diagonal 
block matrix R ^ whoae diagonal blocks are R*^. Form the vector Y by 
"stacking" the observations Y^. 


Ill) Compute the solution U of the system 

tl AR"^1 0 « a r“^ Y ^306) 

The vector G contains the optimal pointwise shape estimates C^ U*(P^). 

Iv) The estimate of the continuous optimal shape distortion AU* is given by 

AU*(P) » I ([ G(P|Q) W(Q) G(q|p )dQ)c/ R?^ (Y. - C, U*(P.)) (307) 

1-1 ■'G X 1 1 1 1 1 

The optimal shape estimate Is U* - 4- AU*. 


The Optimal Shape Control Algorithm 

A A 

Assume again that the positions P^ and matrices C^, R^ and W have 
been chosen. Assume also that the desired shape or optimal shape estimate 
U* is available. Then 


i) 

D 


J 


Compute the block matrices A 
given by 


ij 


given by (305) and the vector elements 




(308) 



I 
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li) Form the block matrix who.<e block components are A^j, 1 1* j 

A A 

Form the block diagonal matrix R whose diagonal elements are R^, and the 
vector D by "stacking" the vectors 

Ill) Solve the system (309) of linear equations for the vector, of optimal 

A 

forces F. 


(R + A) F - D 


(309) 


Iv) Ttte optimal shape correction resulting from t>ic: application of these 


forces at the points is 

dU - f G(P|P.) C.F/ . ( 

1-1 1 1 1 

If the negative of the forces is applied to the shape estimate U*, 

A 

the resulting shape is U* - AU, the optimal corrected shape. 


(310) 


Approximate Algorithms 


Me assume the weighting matrices W and W are chosen to be the mass 
matrix of the dynamical model which corresponds to the static model (253) . 
The eigenfunctions and frequencies for that model satisfy 


“k “ \ \ 


Then an approximate Green's function, based on the first n^^ modes, is 


given by 


■ J, -T *k<f> • 

k.l 


(311) 


Furthermore, the elements A. 


and D., 


determination algorithms are given by 


in the shape control and 


(312) 




013) 
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Substitution of Oil) into the sxprsssion (307) for ths optimal shape 

sstimata U* yields 

a 

■ Ji kh 

k 

Thus the coefficient of the mode in the approxiaste shape estimate 

is 

*i«i - v<'i» • «i5) 

k 


These computed estimsted sK>dal coefficients may be coapared^to the actual 
coefficients of the known distorted shape. Representative cos^arisons 
may be found in the tables at the end of this chapter. 

Substitution of expression (311) into the expresalon (310) of the 
optimal shape correction dr yields 

* “ 1 T 

dU = I I ^ V<^l) C.F, 

1-1 k-1 « 2 * » i i i 


(316) 


Thus the coefficient of the mode iQ <^he optimal shape correction 

dU is 




(317) 


u, 


Comparisons of these terms with the actual coefficients are also found in 
the tables. 


Results of the Simulations 

The tables 6. 1-6. 3 at the end of this chapter exhibit representative 
results for the following choice? of variables. Figures 6.4-6.10 Illustrate 
the results of shape determination and control simulation for selected 
distorted shapes. 

Control/Observation Positions ; The control and observation points were 


chosen colocated both In position and direction. Since conventional 
stability <)uestious do not arise in static problems, colocation serves 
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the convenience of the progremer , but le not neceeeery for eccurecy. 

Either nine or eighteen points were chosen on e given circle* Thus 
they were loceted on every rib or every other rib on the circle, The 
second, fifth, eighth eleventh circles were tried. (Teble 6.1] 

The forces/observstlons were chosen to be ell In the x direction 
■ (1 0 0)), the y direction ■ (010)), the t direction (C^ ■ (001)), 
or both in the x end y directions et esch point (C^ ■ (1 1 0)). Table 6.2 
compares results for the sane shape and varying numbers of points and 
directions. The results for the t direction arc not Included (sec remarks 
below) . 

Modes: The number of modes used in the approximations was either 7 or 11. 

Plots of the first eleven modes are contained In Figure 6.1 - 6.3. 

Welghtina Matrices : The weighting matrix V(P) was chosen to be the mass 

matrix M of the finite element model. This is a natural choice when 
using modes from the same model, since the Inner product for the space 
spanned by the modes Is weighted by M. 

A 

The weights and are scalars in these simulations. They arc 
chosen to be the same nuisber R in both the control and estimation problems. 

The criteria was that R be as small as possible, while large enough that 
the matrix (RI + A) Is invertible. The correct choice of R varies from 
circle to circle, but appears to be half->way in order of magnitude from 
the minimum and maximum elements of the matrix A. 

Observations ; A good test of an estimation alg'trithm is its performance 
when given exact observatloi\s of a known shape distortion. This provides 
a means of comparison of the accuracy of the results. The program was 
provided with the modal coefficients of several known distorted shapes, 
from which it computed exact observations. It uses the exact observations 
in the shape estimation algorithm. 
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It should ho roMMborod vhon oboorvlnt rooulto that tho aedo ohapoo 
roproaont diaplacoaaota of tho antonna froa Ita natural or idoal ahapo 
U*. Thus if 'M ropraacats tho coabioation c£ aodoa la tho dlstcttod 
ahapo of tho antonna, tho actual ahapo la U* dU • U. 

Roaulta 

1) Aa long as tho valuo of tho volgbting factor R la choaon aaall onough. 

It doos not appoar to aoutor on idtich circle the obsorvations aro choaon. 
[Tablo 6.1] There is one oxcoptloa: tho innormost circle aay not bo used, 

hocauao of the assuaptlon that tho hwb la fixed, tho values of all the 
modes on this circle aro aero. 

2} Good results arc obtained fron obscrvatloo/control forces applied 
only in the x direction, or oquivalontly only in the y direction. Thus 
if. observations and/or control forces may be applied in those, or in radial 
direct ior.8, satisfactory roaulta can be obtained. [Table 6.2] 

On the other hand, when observations/control forces vore applied 
It the t direction , results verc very poor (and are not included in the 
tables). Exsoination of the modes reveals two reasons: The first is that 
in the lower order nodes there is very little displacement in the < 
direction. This Is due to the assumption that the ribs are very stiff In com 
parison with the mesh, so the lower order modes consist of ribs being 
pinched together at some points and spaced apart In others. (Figures 6.1- 
6.3). The second reason Is that the changes In the z direction do not 
vary much on the same circle. Control /observation points on two circles 
simultaneously were tried, but results, although better, were still poor. 

For a fixed number of observations, slightly better results are 
obtained If they arc taken at different points In one dlrcctlct', rather 
than in several directions at fewer points [Tables 6.2 an<' 6.3 j. 
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3) Hor« control/obstrvation points than modaa should ba uaad. Aslda from 
cht £acb 'chat this is aasily obsarvad iron tha data, it is a sattar of 
coBBon Sanaa* Both problav^ invol\*a tha dataraination of tha coafficiants 
of aach of tha modaa. Ont must hava at laast as many piacaa of Indapandant 
data as one has unknowns* 

Ht>wavar, It is astimatad that thara will ba from 50 to 150 obaarva- 
tions takan of LSS antanna* Sinca it is unlikaly that 150 >odes will ba, 
or could ba, used in tha modeling, this restriction does not actually 
pose a problem. 


Table Symbols 

The ith mode. 

The rest shape, or ideal shape, of the antenna. 

Ab The modal displacements of the actual distorted shape. 

1) The actual distorted shape: U • U” -f AU. 

AU* The modal displacemettts of the shape estimate. 

U* The estimated shape; 0* ■ U* + AD*. 

AU The modal displacements of the shape resulting from the application of 

the control forces. 

t' The antenna, shape resulting from the applicaticn of control forces: 

0 • U* + AU. 



Table 6.1 
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Table 6,3 

Measureneuts aod Controls applied to both x and y directions at 9 points 
vs. X direction only at 18 points, on the fifth circle 

Actual Shape: U“ + lO^j^ + 

11 nodes used in approximations. R ■ 


Actual Coefficient 

9 X . 9 V 

18 X 

U“ 

U* 

A 

u 

U* 

U 

♦i 

10. 

10.010 

10.010 

10.0 

9.994 

*2 

b. 

.000 

.000 

-.001 

-.001 

*3 

0. 

.ob4 

.004 

-.000 

-.001 


10. 

10.011 

10.010 

9.996 

9.994 


0. 

-.006 

-.006 

.004 

.005 

♦e 

0. 

-.003 

-.003 

.001 

.001 

H| 

0. 

.007 

.007 

.001 

.002 

^8 

5. 

5.255 

5.111 

4.996 

4.994 

4> 

9 

0. 

-.551 

-.559 

.003 

.003 

10 

5. 

4.644 

4.791 

5.003 

5.003 

11 

0. 

.488 

.530 

.000 

mm 


For a fixed number of observations, it appears better to take them at 
different points In the sane direction, rather than to take observations 
of several directions at fewer points. 
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Ideal Shape. First, Second, and Third liodes 
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Figure 6.2 Fourth, Fifth, Sixth, and Seventh Modes 
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ESTIMATED SHAPE 


CORRECTED SHAPE 


Fij^ure 6, 
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CORREaED SHAPE 


ACTUAL SHAPE U® + 20^^ 


ESTIMATED SHAPE 


CORRECTED SHAPE 
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Figure 6.10 O'- + lo + 5 4^. + 3 ^ ^ ^ 
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Chapter 7. Conclusiona and Future Work 

It is possible to accurately determine and control the static shape 
of a large space structure by means of a number of control devices and sensor 
measurements at discrete points along the structure. 

An integral operator approach to the continuous-discrete optimization 
problems of static shape estimation and control proves ideal for these 
problems. Solutions reduce to the solution of linear equations of dimension 
less than or equal to the number of observations, or control forces. 

Elements of the linear equations involve the Green's function, or influ- 
ence coefficient, of the structure, which represents the recponse of the 
structure to a force at one point. In the event that the Green's function 
cannot be computed analytically, approximations based on modal expanaions 
have been presented, involving modes either from the static or associated 
dynamics model, which may be computed experimentally, or numerically. 

The distinction between the shape control system and the attitude 
control, orbit and stationkeeping system arises In connection with the 
rigid body modes of the structure. The rigid body modes represent transla- 
tions and/or rotations lii space of the structure as a whole, clearly a 
concern of the attitude control, orbit and statlonkcepir.g systems. 

Un the other l;:and, the rigid bodv modes are Indetcctable to the sliape 
control system. Furthermore, a shape control system may not apply a net 
force in the direction of a rigid body mode, tc correct It, since this 
would violate the boiutdary assumptions upon which shape control forces are 
computed. The latter restriction places addlMonal cor.iitralnts on the 
stiape control forces in the case that rigid body modes are possible. 

The use of modal expansions for terms In Che sliape control and 
determination algorltlims invites the inevitable trade-off between accuracy 
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and computational difficulty. If a few modes are used and the structural 
distortion involves significant components in higher order modes, the 
shape control and determination schemes will not be accurate. On the other 
hand, the use of many modes increases the necessary storage, time and 
expense of computation. A compensating factor is that while dynamic 
shape control must be accomplished on board the spacecraft, and within a 
short response time, static shape control may be accomplished by ground 
computers over a much longer period of time. Thus, the use of modal 
approximations may not present a difficulty. 

Future Work 

The solutions of both the shape determination and control problems 
depend on the solutions of linear systems which have dimensions on the order of 
the number of observations or control forces applied. It is estimated that 
actual large space antennae will require from 50 to 150 observation points 
for static control* It is therefore desirable to develop a geometric 
scanning algorithm, which would successively process data sets of antenna 
sections in an adaptive manner* 

Despite the fact that linearity and self-adjointness are common 
engineering assumptions, it is probable that large space structures will 
not always have these characteristics. It is anticipated that the integral 
equation techniques used here will be applied to an iterative technique for 
the solution of non-linear problems, and that it will be adapted for 
the solution of non-self-adjoint problems. 
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Appendix A. Some Mathematical Background 
A.l A Little Dletrlbutlon Theory 


We should give some consideration to what is meant mathematically 
by a solution to (13-14) or (24-25) . 

A classical or strict solution to an nth order differential equation 
Lu ■ f is an n times differentiable function y which ’'satisfies'* 
the differential equation: Ly • f on (a,b] . 

Clearly it is not possible for a function to be both n times differen- 
tiable and to exhibit delta function behavior in a combination of its 
derivatives. 

A rigorous development of the theory of solutions of equations of the 
type (13) may be found In distribution theory: 

Distribution theory was developed to provide a rigorous framework 
for "functions" such as the delta function. One cannot deduce from the 
definition 


6(x) - 



X i< 0 

X • 0 



6(x) dx - 1 . 
6(x) 4(x) dx ■ 


♦ ( 0 ) , 


(318) 

(319) 


or even that such expressions are meaningful. Thus a polntwlse definition 
of the 6 function does not characterize It. 

On the other hand. If the 6 function Is defined by (319), the other 
Information about It can be deduced, laus 6(x) Is defined b y its act i on 
o n other functions through the Inner prod uct 

<«,♦> - P ♦(*) 4(«) dx • 4(0) . 
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In distribution theory this concept is extended to an entire collection 

of g eneralized functions , or distributions. Rather than characterizing 

distributions by poiatwlse values, they are defined by their "action" on a 

specific class of functions, called test functions. Test functions are 

I 

infinitely differentiable function'^ on R which vanish outside of some 
bounded domain. Eligible test functions for boundary value prob!lems on tlie 
interval (a,b] must vanish outside of [a,b] . For problems defined on the 
test functions must vanish outside of 

On one dimensional domains, a distribution t "acts" on a test function 
through the inner product 


<t 




t(x) 4>(x) dx . 

Two distributions and t 2 are equal if 
eligible test functions 4* • 

The d er ivati ve of a distribution t is defined by <t’,(^> = <t, 

The nth derivative is defined by - <t- (-1)^ • Note that 

1 j n 

dx 

again the definition describes actions on test functions rather than soitie 
pointwise behavior. 

If U we denote a partial differential operator on U by 

k . 4* . . . k . 




1... x^ 


where K is the vector (k^^...,kj,) and !k! = k^-t . . .+kj . As example of 

3 

this notation, if ^, = 3, a point in R is denoted by (x^,X 2 »X 2 ), and 


K = (2, 0. 3) , then 
7 


D 


K 


3 

“ “ ■ V ■ "5 

3x^ 3x^ 
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In a distribution T acts on a test function t through the inner, product 

<T,*> - I T(P) ♦(?) dP » 

Again, two distribution Tj^ and Tj are equal if for all eligible 

test functions and the derivatives of T are defined by 

<dS:, ♦> - (-1)^*^^ <T, dS> . 

By this new definition of the derivative, since test functions are 
infinitely differentiable, distributions are Infinitely differentiable. 

Finally, the distribution T is a genp. alized solution of LU ■ F if 
<LT,t> ■ <F,4» for all test functions ♦. This removes the problem with 
finding solutions to (13), that is, how a function may be n times differen- 
tiable and yet have delta function behavior in a combination of its 
derivatives. 

If T corresponds to a pointwise defined function which satisfies 
LT » F but is not sufficiently differentiable it is called a weak 
sclution. If T corresponds to a function which is sutficlently differentiable 
sc that the differential operations in LT « F may be performed in the classical 
sense, T is a c lassical solution, or stric t solution. Classical solutions 
are easily shown to be generalized (distributional) solutions, so none of 
these solutions is lost by appealing to distribution theory. 

Ex amp les 

A.l) ^ ^ ® claseical solution t ■ C. It also has the weak 

solution C ■ H(x) (the heavy side step function). 

^2) * ^ ® generelited or distributional solatioa t • d(x), 

which is neither a weak solution nor a strict solution. 

A. 3) Green's functions, which are solutions of Lt ■ d(x-0 are weak 

solutions, since they stay bt defined pointwise but lack sufficient 
diffarentleblllty to be strict solutions. 
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The use of the alternative theorem 6.1, and the assumption of the 
existence of complete orthonormal eigenfunction expansions which are the 
basis of the approximations, depend on the assumption that the operators 
L and K be defined in Hilbert spaces. The Hilbert spaces which can accommodate 
members such as the delta function are known as Sobolev spaces. An excellent 
treatment of Sobolev spaces is contained in [9] . 

A.2 The Free Space Solution of « - ^(pIq) 

The equation 

7^ Y - - fi(PlQ) (320) 

represents the response of a plate in free space at the point P to a unit 
negative impulsive force at Q. 

Theorem ; A fundamental solution of (320) is given by 

Y(x,y,C,n) = log R (321) 

where R is the distance PQ. 

Proof ; We wish to show that (321) defines a solution in the distributional 
sense. Thus it is necessary to show that 

<yS,4» • <Y,(V^)*^> ■ - ♦(Q) 

for all test functions 4i, where the inner produce <u,v> in free space is 
<u,v> - 1 u(P) v(P) dP . 

Let R be a circle of radius e about Q. 
e 




r 

i 
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The tunction (321) is continuous except for a removable singularity at 
R • 0. Thus it is locally integrable and 

Y(P) (vS(P))dP • li® I - Y(P) (A(P))dP. 

J2 C-K) JR'^-R 

Iv V 

We apply Green's theorem, making use of the fact that ^ vanishes for suffi- 
ciently large R to eliminate the surface integral at Infinity. Thus 


y(P) (V (KP))dP 


' 1 

R-^-R 


•^r2-r 


V y(P) ♦(P) dP 


- f ^ (V^)]ds 


J 'an 
3R 


On R^-R . v\ « 0. 
e * * 

the boundary of ds “ cdS and 


1(A) -A(|i)ids. 

The first integral on the right is zero. On 


Therefore 
r2n 


^ r2-r 


y(P) V^^(P)dP 


" 1|r Ir 1 de 

♦ n 


+ e 


.1 


0 

2n 


(|i) -v2* (|l)ldc 


(322) 


S “ 4^ r + |) , 


and 


v\ - ^ (log r 4 1) , 


l-f * 2^ • 
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Furthermore, the test function 4 hes continuous derivatives of all orders 

2 

which have compact support In R . Hence ^ and any of the derivatives are 

2 

bounded on all of R . Thus 


!|^ I 1 IM2 and |||-1 1M3 In . 


We apply these relations to the elements of (322) : 

1 e f Y 1^ (V^4»)de llM, E^log e (2ir) - o(e). 

^0 

I ^ f ^Ir^ I -liT E + 1) M3 (2n) ■ o(c) 

^0 

1 e f^’' vh (|J) de I 1M2 (log € + j) (2it) - 0(E) . 

^0 


Finally, 


-ef 

* n 


♦ (R) 1^ (^^)de 




♦(6) de. 


Taking the limit as E 0, only the last term provides a contribution. 
We can conclude 


y(vS) dQ = - 4.(Q) . # 

■'R 
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Appendix B. The Flexible Beam Program Listings and Output 

B.l The Simply Supported Beam Control Program 

B.2 The Plnned-Free Beam Control Program 

B.3 The Simply Supported Beam Estimation Program 



The Simply Supported Beam Control Program 
Listing 


2 « 

3* 

b* 
6 1 
?♦ 
6 * 
9 * 

n* 

12 * 

13 * 

16 « 
17 * 
I 6 « 
19 * 
2 D» 
21 « 
22 « 
23* 
29 * 
25 * 
26* 
27* 
2 d* 
2 9* 
^3* 
31* 
32* 
33* 
39* 
33* 
3o • 
37* 
3d* 
39* 
43* 
41 * 
* 2 * 
43* 

44 • 
4 *^* 
46 * 
47 * 
4rt* 
4 9 * 
53 * 


CO^^ON STAkTf hSTAI«»HSIM#H^AXf&RP^AX»XLtK£y »X2 au 1 *KS 
ft CAL AU3«13) « bU3)f 40ftKC53) « A9 (13*131 

RLAL 2(50)»U(3Ct) •FSKSOvuLl 
ftcAL AAU0il3) 9d3(10)tUA(30) 

RiAL YZUOtYYUCO 
RtAL QU3) 

CATA YZ/X0*d./ 

DATA 0/13*UeS/ 

DATA Nb»NOA*RC0/ 1«10«I0/ 

DATA ZU) tucn »PSI a>*UAUI/H*34 / 

c 

c 


c **••*• INTROtUCT 10 A ••*••* 

c 

C This PftOGRAr: COPPULS Th£ OPUPAL DlSCfttTi. fCftCCS 

C FOR THE SHAPE CONTROL PKObLEH FOR THE Sl-R^LY >JP=»0RTt3 

C oEAftt Af<C GftAPhS Thi ftESUtUAG SHAPE WS ThE DESIRED SHAPE. 

C THE OUAuRATIC COST IS ALSO C3PPjrE3« 

C 

C 

C PLEASE DEFINE Thi fOLlOylNG 9AP1A3LLS. 

C XL IS THE LlNGTh D? THE SEA^* 

C NH IS THE MURDER CF ACTUATORS* 

C NP IS The NjPdER 3F POINTS ALONj ThE oEAN AT I2HICH Y3J 

L LlSh THE graphs TO toE PLOTTED. 

C NP IS LESS THAN OR EOUAL TO 53. 

C xzn>f l = ARE The actuator PCSITICAS. 

C at CERTAIN X2(I) IS >c.i-JEEN 3 . AND XL. 

C 0(I)| I:]t«..tNH APl the UEIChTS CK TnE rCRCES Fd) SCUAPLO 

C IN THE 3UAORATU JD ST CRITERION. 

C IF C(1);0»ALL If Tfl matrix C*a PAY cE SU6ULAF* RESlLTiNL 

C IN MO SOLJTION. 

c RtCCPHE^D Oz l.E-S • XL**7 

C 
C 

C PLEASE v^hOOSE ONE OF THl FOL.OwiNa OPTIONS. 

C 

C J0PT=1 Only exact d^ti^al forces JILL 5E COYSIDEREO. 

c JOFTr2 t**CT OFTIxAl FORCES AND FIRST A F P A CX * ; I CM 5 

C bASEO 04 lUENFUNCT 104 EXPANSIONS ARE TO 5t w04S1DlPl 


c 

c 

C ThE DESiRED SHAPE IS T«t PARmSOLA Y = LEA«Th*x - X*X 
C 

c 

C PLEASE u^OOSt ONE Of THl :"*OLLOJ14G OPTIONS. 

C KOFI OITERHINES WHAT GRmPh.S APE GEAtRATcw. 

C 

C XOPTjI UL SRAPhS. THl OPTiHAL FCPcESt ShA Ft AD cCST 
L WILL bE ^RINTtO. 
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51* 

C 

K0PT:2 eXACT 0?T 19AL SHAP& VS. SesMCo SHAP£. 

52* 

c 

KOPT-3 OeStRLC VS. APPROXIPATC ShAPCS 

53* 

c 

KOPTra eXACTf APPROX IPATCf A90 0eSIR£0 SHAP£S OM QHi 3VAPH 

54* 

c 

A0PTr5 ttOTH 2 AAU 3. 

55* 

c 

KOPTso 60TK 9 A90 3. 

56* 

c 


57* 

c 


sa* 

59^ 



60^ 


K0PT=9 



J0PT=2 

62^ 


IF(K0PT.GT.2> jOPtr2 

63* 


PU3.U159 

a%* 


6P:20 

63* 


Xtrl03« 

66* 


C£L=XL/RP 

67* 


NP=NP*l 

66* 


URlTe<6*i» XL 

69* 

X 

fORRAT (////lX»23rt7 H£ LENGTH OG TH£ 3 cAR IS*fl3.2l 

70* 


bRUL(6*2» 

71* 

2 

forhat<////j 

72* 

c 


73* 

c 

THl following VARlAdltf ARE VECiSSARY FOR T m£ JPL OJAjRAfJRES 

7H* 

c 

SUbROUriNE. 

75* 

t 


76* 


STARtiJ. 

77* 


HSTARr*a| •XL 

78* 


HRIN^.001*XL 

79* 


hPAX:* )5*XL 

6U* 


LRPAXsi.E-9 

ai* 


KEY53 

62* 

c 


63* 

c 


84* 

c 

These constants are necessary FCR the PLOTTlKii SL6RCUTUE5. 

85* 

c 


86* 


NG^l 

67* 


ncsc'x* 

6 6* 


NT5:-1 

69* 


NTir2 

VO* 


NT2 = 2 

91* 


NT3-2 

V2* 


TlClr«o» 

93* 


T U2-« •• 

94* 


t;l3 = v* 

95* 


XLENre . 

96* 


YLEN=6, 

97* 


00 495 1N-It3 

96* 


NPlIN 

99* 


00 7 NXZ=ltIN 

iOu* 


X2<SXZ»tNX2*XL/ {lk*U 

IJl* 

7 

CONTINJE 

K2* 


CALL V0UT(O*NR» Ua 7hCThE UUG^'TS C(ln 

1 J5* 


CALL VOJT UZ«NP*33 #33HUVhE VECTOR OF ACTJAT3R POSITlOViJ 

104* 

C 


X05* 

c 

HERE XHl exact a matrix ANO 6 VECTOR ARE ODR^JTtO. 

Ic 6* 

c 


1 J7* 

0 



X22 



ioe 

109 

ui 

na 

U3 

I19< 

Xl$ 

116 

117 

116 

119 

120 
121 
122 
123 
129 
123 
126 
127 
126 

129 

130 

131 

132 

133 
139 

135 

136 

137 
136 
13C 

190 

191 

192 

193 
199 

195 

196 
I**? 
196 
199 
15J 

151 

152 

153 
159' 

155 

156 

157 
156 
159 
loO 
161 
1 62 
U3 

169' 


c thc SU6IIOJT1NC 6V£c C09PUTCS th£ cxact 9ecro2 a« 

c 

CALL Bvecidl 
DO 50 izXfhn 

xsxam 

CISX9(X-2.«XL1 
00 53 J=ltN9 
r^xziJ) 
c 

C2-r*Y-2*9XL*t 

C3-X*X^Y*Y 

C9iX»X»Y*Y 

c 

A( 1*J1:^ ( <X*XL |9 <Y«XLl/i36#9XL9XL> 1 « nx«#7 1/7.4 (X # *5 I 9 (C |9C2 1^5*9 

1 IX9«3)9C1«C2/3*) 9 U9Y/06.9XL9XU) •nXL997-Y9 97 1/7.9XL«Y 996-XL«*7 

2 9.2 9 (XL99 5-Y«9 51«< l3«9XL9XL9C31-CXL996-Y9#9l9i3.9XL9939 Xl9C31 

3 9( XL9 93*Y99 31 9U«/3«)9 <5«9XL9XL«C399.9XL*999C9l - tX L •< L -Y ♦ Y 1 ♦( XL • 

9 C9^ (XL9«3MC3l9lXL-Yl9(XL9XL9C9n ♦ ( X/ ( 36 • «XL « «2 1 1 • ( Y -Xt 1 • C ( Y 

5 *97-X«97>77«-XL9*59 (Y9*6-X9 9 61 9.2 •< y 9*5- X 995 I •< 02^2 • I *XL ♦ X*X I • 

6 .259 (Y9*9*X«99)9 13 .9 XL* C2 ^XL^X^X > 9 ( Y •«3^X «93 1«C2 «(X*X«2 *9Xt 9XL 3 / 
73«*( Y9Y-X9X)9.59(XL9X«X«C2II 

50 CONTINUE 

00 51 J-2«NX 
JJrJ-l 

00 51 UltJJ 
A(JtIUA(l»J) 

51 CONTIMJL 
WRm(o»21 

CALL Aour (AiN0Ai 99 t9^tl9*l9iQT9£ £X>CT A *96771x1 
CALL VOUT(o»NM*13«13hOIh£ » VLCTORl 
00 63 1 = 1 »I9A 
00 60 0=ltKH 
A0( I«J >=A(1«J1 
6U CONTINUE 

UO 61 UltN*9 
A(i( It n = AO (ii I) 9 0 m 

0l CONTINUE 

CALL rOUT(6QtNcA«AAtNNt2lt2lHOTH£ EXACT FATXIX C9A) 

C 

C SOR IS A J^L llNcAA LOUATION SOIVUO SUoACLTINE« 

C 

CALL SCR (AC tNCA *NF«6t NOB t Nb t i 3 C tO CR K ) 

CALL rfOUT (aiNN t25# 25H3VLCT09 0- 0?TI«AL •37:£i) 

CO TO 9C 

33 yRITE( 6t3ll 

31 fCRN AT( 15X*26huH ATR IX IS NEARLY SINOOLAR) 

00 TO 5)J 

90 lflJOPT.EU.il 00 TO 170 

C 

c 

0 MERt T.9 l APPROXIRATE VALUES DP A ANJ d ARE CDRPJTtD* 

C 

21=2.9(XL9#7)/ (PI9961 
C 

00 153 UlfNN 
Xrxr ( 11 
UO 153 j:itNR 


f 





169* 


Y:x2(J» 

166* 


AA( l«Jl:Zl«(SlN(Pl«X/Un«SlN(Pt*V/XU 

167* 

I5G 

CONTlMUe 

US* 


WftlTCUtJ) 

169* 


CALL ttOUtUA«.YJA«Y H«N>*«32t32H>)FIRST AA>X)X I XAf 13N TO A XAT4IX) 

170* 

C 


171* 

C 

THC SUS.'iOUTtNe aA> QNAUTiV AN AMYOXINAT£ V-CfOY d« 

I7i* 

c 


173* 


CALL BAP(Att) 

179» 


CALL VOUT(oeiNMt3it32HGFlftST AFPROXIFAUCA TO « VCCTOAI 

175* 


00 163 I=ttNP 

176* 


00 160 J:l>NP 

177* 


A0( ItJI=AA(ltJI 

178* 

160 

CONTINUE 

179* 


00 16$ I:itNN 

160* 


A0( It 1)::AQ (1* 1) * 0 ID 

Idl* 

16S 

CONTINUE 

U2« 


CALL NUUT(AO»NOAiNlltNllt27i27HOThE AFFRCXDATE FATNIX A«C1 

1S3* 

c 


16«* 

c 

SCR IS * JPL LINEAF lOUATION SOLVUo SUbXCUTlNEt 

las* 

c 


166* 


CALL SCRIAOtNOAthPtcB tN&etNet«30tfaCAK) 

167* 


CALL VOUT (»dtNNt2')t26NJVECT3X OF APPX0X1»AT: POTcial 

16b* 

170 

CALL COSTIetOtNPtO 

169* 


WRlTL(6tl7D C 

19C* 

171 

F0RPAT(///t7X»17hTHC CXACT COST ISt EIS.S i 

191* 


IFiJOPTttO.D 00 TO 17$ 

19?* 


CALL COSTUBtOiNP.tC) 

193* 


URlTei6tl72> C 

196* 

172 

F0RPAT|/,VtlX*23hThE APPROXIFaTC CCSI IS * ElStS) 

195* 


URITtI otl9«) 

196* 

175 

IFIJOPT .eO .11 WRDL(6*l$$i 

197* 


WRlTc(6t21 

1 rb» 

C 


199* 

C 

HERL T.'IL shapes A4I CONPUTcu* 

2v0* 

c 


2Jl* 

c 

2(11 IS the X YALJi OF T Ht ITh »01MT OY A .*,IA»H. 

• 

c 

li(l> IS The V VALLE OF The ITh FOUT CN THE GRAFH CF THE 

2w5* 

c 

OPT DAL SHAPE. 

20b* 

c 

L»U> IS THE V VALLE OF THE ITH FOIM CA THC ORAFH. CF 

205* 

c 

APPftOXtHATE SHAPE. 

206* 

c 

psmi iS The V VALUE of the ITh FCUT CA The CifAPH CF 

207* 

c 

OESIRED SHAPE. 

20 6* 


CC EOO A:2 iKP 

2 J9* 


Z(Kl-<*-ll • 3EL 

2K* 


X = ZIK 1 

211 * 


PSIlKl :XL«A-X*X 

212* 


UIKtrC. 

213* 


UA(Ki:j. 

2 1 b* 


LO IVC UltNP 

215* 


IF( X.GT.XZI III 03 TO IcO 

2 1 6* 


r.:(XZ(;i-XL)«X» IX»X-2.*X2UI*Xl»XZII)*»2 1 

217* 


uC TO 13'. 

2 1 n* 

lov 

u:<X-XLI*XZ(ll.|XZ(il»>i:>!.«XL*X*X.V) 

21 V* 

US 

L'(K|:0(AI*j.^(i>/(o.*XL> 

22 j* 


IF ( JOF T .Nt .2 I 00 TO DO 

221 * 


UAlAUjAKI ♦ C 1/ (».*XLI 
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dZ2* 

IVO 

ZZ5* 


22<i* 

&91 

ZZU 



192 

ZZl* 

199 

22S* 


229* 

195 

232* 

196 

231* 

200 

232* 


233* 


23*» 


235* 


23^* 

215 

237* 


232^ 

C 

23V# 

250 

2V2# 


2V1# 


2«2# 


2V3# 


2H%# 


2<i3* 


246# 


2%7# 


2 # 6 ^ 


2««V# 


257* 


251* 


232# 


233* 


25»»* 

276 

25b* 


^53* 

273 

257* 


256* 

271 

25V* 


26) * 

272 

26 1* 


264 • 

273 

263* 


2 6** 

274 

265* 


266* 

275 

267* 

2bC 

26$ * 


2v V* 

3vO 

2 73* 


271* 


2 72* 


27 3* 


27** 


2 7 5* 


2 76* 


2 7 7* 


2 7d* 



CONTINue 

WRirC(6»t«»t XtPSl (K)«U<Nt 
to TO iCO 

W01TCU*1T»> XtPSt (TtfUUOt'JAIK) 

FOMAT|lhO»SX*OHPC& lTlCI«*«Xti;»OeSlACC ShAAe*07X*SftSHAfC»2x» 

I tIHAPPROX. SMAPCI 

FOMATllHOOXtAHPCS IT tON»AX tUtOCSl ACO SHAPt*0TX«9M$MAPC I 

POPPATi^ tlXtPD.atSElS.Ol 

COMT IMUC 

NP2:2*NP 

bO 219 l:ltNP 

TT( I> :p$im*l .2J 

TVI (•>TP>:U(t)*1.2! 

cotTiKiue 

to TO (900t300»3iC*29C«900«290l«KOPT 
CALL liuNPLT 

CALL 5‘LFOAPT*; IN. I • *XLCX. YLtNI 
CALL 'aLSCAL (2»AP-AL>*TV.NP2*A6> 

CALL ?LAaet(*THt INPLV SUPPORTCti 0CA9* .29t« LctaTN IN P'TEaS'tl^* 

1 ‘blSPLACCPENT* »li» 

CALL P.6RAF 

CALL P>.AXl!>t~2<XLLt<0.t 
CALL PLCJAYtZtUiT* tNTltTlCl) 

CALL PLCURV(2>P9I«APtNT2*TIC2> 

CALL PLCJRV4 Z«UA»1P«NT3tTlC3t 
CALL PLCUAV (X2*Y2«NH«KT9tTlC9> 

call PlT^XT (2.».>.5« .ltO.>33HA:TUAT3A AOilflOVS PATXf.3 SI <*39tll 
CALL PLTEX'; (t .3»5 4>«.lJ»0.«4RhLESlR£3<*> NS OPTlPAL<Ot V! APPtCXlP 
1 ATc(/I S.4A>ES»4««1 > 
tb TO l2TOi27l«<T2« V73>214«2T9»2T6I*AP 
CALL PLTCXT (3.2i7. Ui .1»0.« AS.-A-Sl'/cN ACT JAT 3TS • I S* 11 
to TO ^60 

CALL PLTEXf O.Utl.Jt.ltO.f 12H04t ACT'J AT OTt 12 « 1 1 
to TO < dC 

CALL >LT-:xr<1.3t7. 3*.ltO.« llhTdO ALU* AT OTS • 1 3t 1 1 
to TO 2o0 

CALL ALTcXT(3.2«T.0t.l.0.> onTlitEE ACT d AT STS > 1 I i 
Ob TO xdb 

call PlTlXI (3.3.T. J* .1iO.«1 -HF3UA ACT J AT 3 H • 1 A 1 1 T 
bO T C < oO 

CALL ALT£XT(1.3«r.3«.l<0.«l4irIVC ACTUkT STS • 14« 1 ) 
to TO kdO 

call AlT^XT ( 3. j«T. Jt.l>3.«13MStX ACTU AT OT>t 1 3« 1 > 

IMKOAT.Cb.Ht itb TO ATu 
CALL AvYPLT 
CALL fcuT.HT 

call >Lf 3R.HT •. U.l 4* .XLCN. 'LtNt 
Call PL S L AL ( Z ♦**P * Ab ♦ T T * N a 2 • Lb t 

CALL PLAdC. I*ThC SIHPLY SUP*0TTl 3 ot AT* *Z>. ♦ *.i VJ TT IN AiTcTS'tUt 
I •blSPLACt>ENT‘.li» 
call PLvjTAA 

call pl axis <-:• xlca ao . > 
call »l CJPWCZ.uM* .nti.tuii 
call PL COA V TZ .PK I t». a .KT2.T Il2 » 

CALL »LCJ.TVTXZ.YZ»NH*NI J«riCSt 
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a79« 

afto* 

2tS* 
29^* 
2«»« 
2«6* 
2»7* 
266* 
269* 
290* 
291* 
292* 
29J* 
29H* 
295* 
296* 
297* 
296* 
299* 
300* 
301* 
302* 
303* 
3U6* 
3J5* 
306* 
3JT* 
306* 
309* 
310* 
311* 
u e* 

3i3* 

31 *• 
31S* 
316* 
317* 
3U* 
319* 
3*0* 
321* 
32i* 
323* 
329* 
32b* 

32 6* 
327* 
32 «• 
329* 
330* 


CALL rLr£xri«.Si$.5« .Itw.tJlHACTUATOR AOlirtOML HA*Kr.3 4I 
CALL fLTCXt(2^«94>«.>O»0.»ltHCCStflC0(*> «S OPflXALtC) SHAPCS01« 
in 

60 to tiaot J 2 t. ) 2 <» 92302A»};9«32A>*Af 

326 CALL ALT£xr<3*2'>7.4* .1*0.»1)HS£VC« ACT'J At O^t • 1 »• 1) 

CO to 33C 

320 CALL ALT£Xr(3.4«7.O*.l«O.<12M0N£ ACtU At 0t« 1 2« 1 > 

60 to 330 

321 CALL ^LtCXT(9.3*7.0«.l*O.«13MTW0 ACTUAT OtS 1 1 3* 1> 

60 to 330 

322 CALL PLtLXt<3.2*7.0* •nO.tlSHtHXEI ACtU AT 0A> 1 1 3»m 
60 to 330 

323 CALL ALT£Xt(3«3«7.0«,ltO.<lt‘<'4uf( ACT J AT 32 j « 14 • 1) 

60 to 330 

324 CALL rLt£XT(3.3«7,Ot .l*0.«14Hf 1V£ ACTU At 3A3 1 1 4 • I ) 

60 to 330 

329 CALL >LtcXn3.3»7.0»-nO«»13M91X ACTUAt Oil « 1 3t 1> 

330 If (K0PT.C6.2> :«0 10 49v 
CALL AwVALt 

390 CALL eOLPLT 

CALL PLfOR«(*L IV.U»tSLt.N*rLENJ 
CALL PLSCAt (?;4f*66*VVtNP2«Lvl 

CALI PLA3£L(«TH£ S IVLt SUP^OStEO a£A X* »2 9« • L£'< 3 14 19 XitLXSSloi 
1 *Ul9PLACeP£Nt*«12l 
call PL6XAF 

CALL PLAXlSI-2«XltN »0.> 

CALL PLCJPt<2*PSI«NP.ST2*TlC2» 

CALL PLCURV(Z»OA»LP .NT3»tIC3) 

CALL PLCJPt<XZ.y£.N*fVI9.Ti:5l 

CALL PLTCXTJ2 .5*9.5. .1.0..3CMACTUATCX fCSITlCLi Pi«KtC oY X.3C.1! 
CALL PLr£Xr<2.3«9. a* .lo«3.*39Hj£St t£J( «i Vi APP40X 1 XAt • I / 1 SHA>r>. 
1 35*n 

60 TO i37J.371.372.3T3.374.3V5. 376l *S1 
376 CALL PLTCXt(3.2*7.0*.1.0..l9Hi£y£N ACtLfe T Cf S t la « 1 ) 
bO to 4«3 

370 CALI PL1CXr<3.4.7.u*.t«O..l<H0t.i ACtlAT Cftll • 1 ) 

60 to 4«3 

371 CALL PLtexI (3 .3. 7 .j« .1.0.«13 hTLO ACtU A T C fS . 1 3 . 1 I 
63 to 4tJ 

372 C«LL PLT£Xt<3.2*7.0*.1.0..1iHTV2££ ACU A t Cf * * 1 a . 1 » 
bC TO <>43 

373 CALL PL text (3 .3. 7 .0*. 1.0. .i«Hf CUR ACT U AT C A 5 . 1 4 * 1 > 
bO TO 443 

374 CALL PltLXtt3.3.7AJ«.1.0..1«Hf ive ACTL At C Ai • 1 4 . 17 
iO TO '.43 

375 CALI PLTeXT(3*3*7.b. .1.0..t!MSlX AUUAtCAi.l’ .M 
4V0 COLT 19. 4 

CALL Ai,VfL7 
44a C09U9JC 

CALL tLCPl 7 

t OObt fOSutT TO lit*. ACE 44J CALI tVUPLT 
OUO Step 


J3i* 




V3 Of 1.34AUATJ09; 40 

732 i,tp;.j»l SUPilT.btC 


6U3935TiC5 



1« SUbROlTlNE bAPtO) 

2* C0MJ10N STAkTtflSTA^ tE^MAXtXLt AEYtXZU J) « Vi 

3* fiiAL 8(ld> 


H • 
$♦ 
6* 
7* 

a* 

10* 

12 * 

1 3* 
14* 
IS* 
l6* 
17* 
13- 
IV* 
ii3* 

2 1 * 
22 * 
23* 


C 

t This sudkoutine ccrpuies an approxipatc t vector for 
T rtt SIMPLY SJPPORTiU oiAM. 

C 

C INTEGRATIONS ARt P£RFOR.iEO 6Y THE JPL OJADRAtJRES SJaROiTlNE* 
t ftOPbS AND R0P2 ARE PART OF THAT SUbRCUTlNE* 

C 

L PSl IS The LcSIREL shape* 

c 

PI-3.14159 
00 53 1=1 fVi 

C v^L ROM bS (ST ART *K »X *f QFX t PST ARt HPl NthPAX tc RPaX« AAS« R * KEY ) 
10 PSI=XL*X-X*< 

F0FX = PS1*S1N tPl*X/XL) 

CALL R0M2 

IF(K,eO.U go to 10 

d( 1) = A.\S*SiN tP I*<2 a )/XL) *2.* (XL** 3J /(?I**4 I 
50 CONT INuE 

RETURN 
ENL 


1* 


SOoROJT IN£ dY£C(i) 

2* 


CCPPCN STAKTfHSTAF* PM I,'. • pPA X * £ RPA X t XL t XE Y • XZ ( lo JtKF 

3* 


REAL cJ(U) 

«* • 

c 


5* 

L 

This S’JoROUTINE C34PJT£S THE EXACT 6 VELT39 ?0R T HE 

6 • 

C 

SIMPLY SOPPORTEu giMM. 

7* 

c 


b * 

L 

iMtbRATlCNS ARt FcRfCRfEL oT T r-£ wPL CLALRATURlS SLcRClTUL 

9* 

L 

RO*oS A;.0 R042 kki PART OF THAT SuS^^aUTiNE* 

1 U* 

0 


11 • 

c 

PSl IS I He. jESIRE3 SriA^t* 

12* 

L 


13* 


uO 5j 

1 4* 


2 = a2 ( 1) 

15* 


CmLL Pu'ljS(STAKT»XL.X>r'OFX*HSTAR#H*lVf •iR^AX#A'^S*<*<EY) 

1 0* 

lo 

IF ( , GT .2) oO TO !:» 

17* 


0=(2-XL^*X*(X*X-2.*2*XL*2*2) 

1 0* 


GO TO 2^ 

19* 

l5 

ur(X*XL)*2*^2*2-2.*yL*X*X»X) 

C • 

2 C 

o=G/ (6**XL ) 

2l • 


P:=3.1^15V 

< i * 


PSI=X*XL-X.X 

23* 


fOrX-P;>I*.j 

2 4* 


CALL R-.*^2 



IF ( ^ . 1 ) GO ro i i 

2c* 


b( i ) rk .t 

27. 

? J 

COMIN.t. 

2 ^ * 


Ft Tl RN 



1* 


SUttROUUNe C0STUt«tN)>»O 

2 ^ 


CO!<NON STAi(Tt(4ST«A thillNt HMAX *£^!<AX *XLt KCVtXZ U J ) 

3* 

U m 


REAL BUO)tQ(ID) 

5* 

C 

c 

THIS SUllROUTIMC COPPUTCS tHl QUADRATIC COST FUNuTICAAL 


c 

FOR A SLT OF NR FORCES »tl) AT ROSIUQXS XZ(I) AIQN<> 


c 

A SIMPLY SUPPORTED BEAM 


c 


9* 

c 

C0SX=U/2)« (SUMtC (It *0 Ui**2 ) • 1 NTEC RAL ( (U ( X I -P<SI (X 1 1 • *2 1 1 

IJ* 

c 


U* 

c 

UNERE U IS THE OPTIMAL SHAPE ARC PSI IS THE UcSlFEO SHAPE 

12* 

c* 


13* 


C=U. 

19* 


00 S I=1«NH 

15* 


t:C»0 

16* 

3 

CONTINUE 

17* 


CALL RONbS (START •XL*X»FQFX»HSTAR<hPIA»HRAXtLfvM*Xt AAS»K«KEV 1 

18* 

14 

SHAPtsJ. 

19^ 


CO 200 UltNM 

20* 


IF(X.ST.X2(I») GO TO IdO 

21* 


6-(X2 |l»-XLt*X* (X«X-2.«X2(l>*XL*X2 ( 1 )»*2 1 

22t 


GO TO 190 

23* 

180 

G=(X-XL1*X2 (It* (X2( II**2-2.*XL«X*X*X ) 


190 

SHAPE:SriAPt*G*&(II /(A.*XLt 

23* 

2U0 

CONT iNUe 



Pl=3. 1*159 

27* 


PSI:X*XL-X*X 

28 * 


fOFX:.>*(SttAP£-PSI t**2 

29* 


CALL R0M2 


3U* 

UU.LU.UoO TO K 

31* 

Cro*C *• A»\(S 

32* 

8l TURN 

33* 

tN J 
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Output 


THE 

L£N6TH 

OG THC 6EAH U 

tOO.uO 


THE 

UeiGHTS 

0(1 > 



1 

TO 

3 

t. OJOOO JO*OS 

(.0030000*}* 

1 .0000033*30 

THE 

VECTOR 

OF actuator positions 


1 

TO 

3 

2.5003000*01 

S*00w00304>31 

7.S030300*0( 


THE exact a >^aTRIX 


ROW 1 
ROW 2 
ROW 3 

THE e VECTOR 
I TO 3 


COL 1 
J6/3027^ Iv) 
1 .A 9 C 2 I 27 ^ 1 C 
l«OH5aJ22>iO 


1. 3737dJ9^0V 


COL 2 
X.*I 9 D 2 ; 27»13 
i* 10 ©l 3 CM 4 )C 
1*4902Un^13 


2 « 647573 U^ 0 V 


COL 3 
1 • 0450372»10 
1 . 4902 M 9410 
1«0622936*1J 


t .0737792»09 


THE exact matrix (-♦a 


ROW \ 
ROW 2 
ROW 3 


COL i 
\ *072302741G 
l.4932l27#^lj 
I* 045 a 322 ^l 0 


COL 2 
1 .49C212741C 
^.lldl 334 ^n 
1 «49C2iI9«lC 


COL 3 
1.045b322«10 
1 «49021 19^10 
l.072295b«10 


VECTOR OF optimal FORCES 
I TO 3 5*2l36d3a-02 


S.lo355C7-C2 5 *2 1 34 i> -02 * 


! 


FIRST APPRUXIMATIGM TO A MAlRlX 


ROW 1 
ROW 2 
ROW 3 


ROW I 
ROi« 2 
ROW 3 

VECTOR OF 
I TO 3 


COL 1 
1 *OS 39096^10 
1* 49^45434 I J 
1 •J 539 l 2 d«l 4 . 


COl \ 

1 . 06390 V 6 * I J 
I.HV04S43*! 0 
U 0359125 * 1 J 


COL 2 
I *490454^ ♦IC 
2 *l 07 d 222 ^l J 
l* 49 G 45 d 3 ^U 


lOL 2 
1 *49J4S^ 3 

^ .11 7o222 «1C 
U 49 j 4583 ^i J 


N*20d2374-32 


COL 3 
I *0539125 «10 
1 * 49345 o 3 ^n 
1*0539153*10 


1 *8729590407 


wOL 3 
I *05391 25M0 
1 . 49 C 45 b 34 l 0 
I *0o39l 5341 J 


4 .4323^45* 02 


^FPRCXl**ATt FCRLlS 
4 *a 3234 9 t,-j 2 


FIRST APPROXIMATION TO d VECTOR 

I TO 3 Ud7295394C9 ^*64t75874C9 

The approximate MmTRIX A^U 


t 
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THE exact cost IS 


THE APPROXIMATE COST tS .A350S»le 


POSITION 

OeSlReU SHAPE 

SHAPE 

APPROX. SHAPE 

5*00 

.47500^03 

•4C412*03 

.402C7*C3 

io«oa 

.93000403 

.79850*03 

•79459*33 

15.00 

.12750*04 

.11734*04 

. 11666 4C4 

20. OU 

•16000404 

i 15191404 

•15137*34 

25.00 

. 16750404 

•16257*04 

.18212 *C4 

30. OJ 

•21000404 

.20648*04 

.20823*04 

35.00 

.22750404 

.22919*04 

.22923*^4 

40. J'J 

.24000*04 

.24438*04 

.24470*34 

45.00 

.26750404 

.25373*04 

• 254.;tj *C4 

50.00 

.25000404 

.25692 *J4 

.25 756*04 

55.00 

.24750404 

.25373*04 

•2542b*C4 

60. 0 3 

•240004 J4 

.24438*04 

.24470*34 

65.00 

.22750*04 

.22919*04 

.22923 *C4 

70. JO 

•21330*04 

•20848*04 

.2082 3*34 

75.00 

•1675u404 

•16257*04 

. 162 U *l4 

dO.OO 

« lo0a0*D4 

.15190*04 

.15137*04 

65 . 00 

* 12 750^04 

•11734*C4 

•Il6b4*l4 

SO. JJ 

.90700*03 

. 79853*35 

.79459*3 3 

95.00 

• 4 750u* 03 

,40412*02 

.4C2l7 *C3 

1 J 0 . J J 

. 30300 

*00003 

.00030 

PONc^*t 

19d/oH6»o9PLTj 




A^lOTtP 
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B.2 The Plnned-Free Beam Control Program Listing 

I 


2* 

^ft 

4* 

7* 

bft 

n* 

n* 

12« 

15^ 

14« 

17ft 

18* 

19ft 

2D* 

21^ 

22ft 

25* 

24ft 
25* 
24* 
27* 
28* 
29* 
50* 
31* 
32* 
33* 
34* 
35* 
34* 
37* 
3h* 
31* 
4 • 
41 • 
42* 
43* 

44 • 

45 ft 
4f ft 
47ft 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


COHHON ST4RTftHST44(Hf^lllftHMAX«CRM4X«Jlt«KCY«X2tm«MM 
RC4L AUQflOftRaO) 

Rt4U YUOl) ftO4<S0) 

RC41 4R t 1 IR) «A2 (19*10) ftB! 6(1 3>ftQ CIO) 

R£*L 44(10* 10) *RH1 ( 1 0) *04 (10) *FA ( 1 1) 

RC4L 404(10*1 3)*B04ilQ)* 404(10*10) 

DtRtNSlON WOOKdOO) 

RC4L F(10)*U(90)*RS1(OJ)«X(50) 

RCAL FH12(10) 

D4T4 C/10*0«d/ 

THIS FR06R4M COMPUUS THC OPTIMAL OISCRtTL FORCCS 

FOR THC SHAPC COMTROL PROBLCH FOR THC ^PL FLCXIOLt 

eCAM* AND ORAPHS THC RCSULTIN6 SHAPC VS THC OCSIRCC SHAPC* 

PLCASC DCFINC THC FOlLOWIMO VARIABLCS* 

XL IS THF LCN3TH OF THC OCAM. 

HP IS THC NUMRCR of ACTUATORS. 

NM MUST eC SRtATCR THAM OR CQUAL TO 2. 

HP IS THC NUMRCR OF POINTS ALONG THC 8CAM AT MHICH YOU 
WISH THC graphs TO Bt PLOTTCDe 
NP IS LESS than or COUAL TO 5Q * 

XZU)* l:l*.ee*NM ARC THC ACTUATOR POSITIONS* 
hi CERTAIN X2(l) IS BfTVCCN 0. AND XL* 

Q(U* T = l*ee**NM ARC THC UCIGHTS ON THC FORCES F(l) SGUARCO 
IK THC QUADRATIC COST CRITCRION# 


PLCASC CHOOSE ONE OF THE FOLLOVING OPTIONS* 

JOPTsl ONLY EXACT OPTIMAL FORCCS WILL BE CONSlDCRCO* 

UOPTS2 C(ACT OPTIMAL FORCCS AND FIRST APPROXIMATIONS 

BASED ON EIGENFUNCTION EXPANSIONS ARE TO RC CONSIDCRCO* 


THC OCSlRtO SHAPE IS THE PARABOLA Ys LENGTH ♦ 3/4 t • X*K* 


PLEASE CHOOSE ONE OF THC FOLLOWING OPTIONS. 

KCPT DETCRMINES WHAT GRAPHS AHf GCNCRATCDe 

KOPT::! NO GRAPHS# THE OP* iMAL FORCES* SHAPE ANO COST 
WILL BE PRINTED. 

RCPTr? EXACT OPTIMAL SHAPC VS. DESIRED SHAPE. 

KOPTs' DESIRED VS. APPROXIMATE SHAPES 
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C 

K0PT>« CXACt* APPROIINAU* AliO OCSIRCO SMAPtS CN ONE GRAPH. 


C 

KOPTsS BOTH 2 AND 3* 

• i • 

C 

KOPTsG BOTH A ANO B. 


c 


53* 

c 


55 * 


KOPTsA 

5 ^* 


J0PT*2 

55 * 


ir(K0PT.GT.2l «I0 PTb2 

56 * 


NH>2 

57* 


*Lsl'»0. 

56* 


HP* 21 

59 * 


DO 1 UX.NH 

60* 


T2mtI*.S •Xt 

61* 

\ 

COHTlHUr 

63* 


HR1TC(A«2> XL 

63* 

2 

F0RHAT(////«1X.23HTHC LENGTH OF THE BCAH IS«Flt.2) 

66* 


CXLL VOUT(X2«NH«33«33HOTHE VCCTOR OF ACTUATOR POSITIONS) 

65* 


DELsXL/NP 

t6* 


NPsNP*i 

67* 

c 


66* 

V 

THESE CONSTANTS ARE NECESSARY FOR THE PE0TT1N6 SUBROUTINES. 

55* 

c 


7:* 


*LEN»8. 

7J* 


YLENsA. 

7?* 


N6S) 

73* 


NT*(> 

76* 


TlCls*** 

75* 


T!C2s*0* 

76* 


TIC5s»/» 

77* 


T1CR=«*» 

7B* 

c 


79* 

c 

These constants are necessary for the natrix inyersicn routine 

6 r * 

c 


n* 


LO»sl'' 

B2‘ 


NCB:l 

63* 


NPsl 

86* 


XJIX2U) 

85* 


MsNH-1 

86* 

c 


87* 

c 

the fOLLOHlNG VARIABLES ARE NECESSARY FOR THE UPL QUADRATURES 

8** 

c 

SUBROUTINE. 

H9* 

c 


9v * 


STARTsC. 

*»1* 


HSTAR=.01*XL 

92 « 


HNINiXL*I.E-A 

93* 


HNAX* . t5«YL 

96* 


ERRAXsl.E-9 

• 


RETS? 

96* 

c 


97* 

c 


98* 

c 

HERE THE EXACT LITTLE A HATRIX AND B VECTOR ARC COHPUTCD. 

09* 

c 


; ^ 


CALL ARAT(AI 

1 M • 


CALL FVCCte) 

1^2* 


CALL BOUT (A.NDA. VH.NH.20.23H0THC LITTLE A MATRIX) 

i03* 


CALL V0UT(H.NN.20«2)H0THC LITTLE B VECTOR) 

K6* 

c 
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105* 

c 

ns* 

c 

1*^7* 

i^e* 

no* 

111* 

in* 

75 

113* 

m* 

115* 

11b* 

C 

117* 

C 

ns* 

c 

119# 

120* 

121* 

122* 

6a 

125* 

129* 

125* 

65 

126* 

1?7* 

C 

126* 

C 

129* 

c 

15 V* 
131* 
152* 

90 

133* 

91 

134* 

155* 

95 

135* 

C 

137* 

c 

136* 

c 

139* 

19U* 

m* 

n:* 

193* 

n 

299* 

195* 

c 

195* 

1*7* 

c 

196* 

c 

199* 

15 

iSl* 

n^2* 

ID 

nr* 

C 

159 • 

155* 

155* 

C 

157* 

c 

ife* 

c 

159* 

c 

n:* 

c 

IM * 

c 


HERE THE 81$ A MATRIX AND 8 VECTOR ARE COHRUTCO. 

00 !•> I*2*MN 

niO<l>lls8<I»>K2in*BU)/Xl 
00 7S JS2«NN 

AB(I«l*J*l)sO<l>*XZ(t>*XZ(J)/(Xl-Xl)*Allt«t>-AO*U)*A2(l)/Xl«Aa«l) 
l*X7tJ)/Xl ♦Aa«lMXZ<n*X2< J)/(X1*X1 1 
COMTINUC 

CALL HOUT«AB*NOA«N*H«17*17HOTHE 81$ A MATRIX) 

CALL V0UT(e)$«<*«17«17MitTHE BTC 8 VECTOR) 

CALL V0UT«O«NM,28«2aH0FOR THIS UEICNTINB VICTOR 0) 

HERE the EXACT UEIGHTEO MATRIX A«8 IS COMPUTED* 

00 80 1:)«H 
00 80 Jsl«M 
A0tltU)iAM<l«J) 
rOMTlNUE 
00 85 lsl*M 

A0( l«l)sAQ< 1*))*0«1«1) 

COVTINUE 

CALL H0UTUQ*N0A«M«M«2X*2XH)THE MATRIX 81$ A PLUS Q) 

MOU WE SOLVE FOR THE EXACT OPTIMAL FORCES F2 TO FM* 

CALL SOR«Aa«NDA*V«81$«MOE«NB«$0$*VORK) 

CO TO 95 
UR1TE($*91) 

FORHAT<////*10X«25HMATRIX IS NEARLY SINGULAR) 
fO TO 55f 

CALL VOUT<B16«M*20«20HOTHE FORCES F2 TO FM) 

WE COMPUTE THE ENTIRE VECTOR OF OPTIMAL FORCES* 

F«l):3. 

00 I'JC 1=1«M 

Fll ) = M1)-BIC(1)*X2(1«1)/X1 
F<I*1)=RIGI 1) 

COVTINUE 

CALL VPUTtF«NN«2V*29H8THE VECTOR OF OPTIMAL FORCES) 

IF( JOPT.EO*!) GO TO ]75 
••****THE APPROXIMATIONS****** 

no 1 C 6 1 = 1*10 

0(1 )=(XL**7)*1.E>7 
CONTINUE 


V=J,927 

V2=7.;t9 


V A(0 V2 satisfy tan V = TANH V. 

THE FIRST FICENVALUF IS IV/XL)**A. 

THE SECOND eigenvalue IS (V2/XL) ** A. 

FIRST COMPUTE THE E 1 CEKFUNCT ION VALUES AT X2(l)* 
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162* 

163* 

166* 

16S* 

166* 

167* 

168* 

C 

16“^* 

113 

170* 

C 

171* 

C 

172* 

173* 

176* 

175* 

176* 

177« 

C 

178* 

123 

179* 

C 

180* 

c 

181* 

162* 

163* 

186* 

c 

185* 

186* 

187* 

188* 

189* 

198* 

191* 

15 

192* 

193* 

196* 

195* 

130 

196* 

C 

197* 

C 

198* 

199* 

2?C* 

201* 

202* 

233* 

C 

2C6* 

2L5* 

2:6* 

16C 

2;-7* 

c 

2 08 * 

c 

2C9* 
2 1 (‘ * 
211* 
212* 

c 

213* 

216* 

215* 

150 

216* 

217* 

218* 

155 


00 110 1«1*N« 

OKCtVtKZCn/Xl 

*R02«va*xzf n/xL 

f>Hl(l>>t*l. 41421 *SINU«C)*(1.969SC'>2>*(CX^URKI-CXI>< 'ARC) > 

PNl2Cl>«(l*414«)*SlNC«llfi21«<*a9lie*S)*<CXP««RC21-CXP«*ARS2)l 

CONTINUE 

THE APPROXIMATE LITTLE A MATRIX* 

DO 12U li<t«N4 
00 123 Usl.NM 

AA« 1*JUMXL**7)/CU**0))«PH1< n*PHKJ) 

AA(1« J>sAA(I«U)*l (XL**T>/(V2**an*PMI2(l)*PNI2(J) 

CONTINUE 

NOW the approximate little B VECTOR* 

DO ISO Isl«NM 
KET*5 

CALL ROMBS C START* XL«T«F0FT*HSTAR«HH1N«HNAX*ERHAX*ANS*K*KE*) 

WANTS. TS*XL«T-T«T 

AR6sv*T/XL 

Ps«>1.4142>*SINlARCI**019A9S*tCXP(ARA)»EXP«>AR6li 

FOrT«W»NT«P 

CALL R0M2 

IF K.EO.l) CO TO 10 

RA< ns(XL**3)*PHItll*ANS/<V**4) 

CONTINUE 

CALL V0UTIPHI*NN«19*lSH(tTHE PHI VEC/ORl 

CALL M0UTIAA*N0A*NN«NN«32*32H0THE APPROXIHATE LITTLE A MATRIX) 
CALL VOUT(BA«NM«S2*S2HOTHE APPROXIMATE LITTLE 9 VECTOR) 

HERE WE COMPUTE THE BIO APPROXIMATE A ANO S* 

DO 140 Isa,N4 

BBA(I-1)sPA<| )-XZU)*8A(l)/Xl 
CO 140 Js?*NM 

ABAU'l* J>l)Ea<ll*XZ< I)*XZ<U)/<X1*X1 l*AA(l« J)>AA(1*J)*XZI1)/X1 
1 -AA(1*I)*XZIJ)/X1«AA41«1)*X241)*XZ<J)/<X1*A1) 

CONTINUE 

CALL MOUT(ABA*NOA«N*N*24*24H}THE BIC APFROX A MATRIX) 

CALL VCUTIbBA*Mt24*24H0THE BIC APPROX B VECTOR* 

HERE THE APPROXIHATE HEICHTEO MATRIX BIG A « fl IS COMPUTED* 

DO 150 lsl,M 
CO 153 Jsl.M 
AOAIItJlsABAIItJ) 

continue 
DO 155 Isl.P 

AeAn«l>sAQACI*I>«0<l«l) 

CONTINUE 

CALL MOUT<AQA.N0A*M*M«31*31H3ThE APPNOX MATRIX BIG A PLUS 3) 
CALL $ORiAQA,NOA«M*8eA*N06*N8«$90*WORK) 
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219« 

160 

CALL V0UTIBeAfR«32«32H0THC APfROXlMATC FOIICCS F2 TO mi 

220* 


FAUMO* 

221* 


00 170 1»1«N 

222* 


FAIl>srAUI»0BACI)*A2tl*n/n 

225* 


FAU*nsBBAin 

224* 

170 

CONTtfUUC 

229* 


CALL V0UT«FAfNM«31OlH0THC APPAOXIMATC FORCE VECTOR F> 

226* 


GO TO 183 

22t* 

175 

WRlTt(6»17a> 

228* 


CO TO 185 

229* 

178 

FORMAT(////«3X*8HPO$mOH«4X«13HD£SlREO SHAPCf 2X 1 13HOP T I9AL SHAPO 

250* 

179 

rORHAT(////«3X«8HPO$mON«4H«lSHOCSIRCO SHAPE •2X«15H0PT15AL SHAPE* 

231* 


1 2X*13HAPPROX« SHAPEI 

232* 

180 

WRITE46*179) 

253* 

C 


234* 

c 

HERE yC COMPUTE THE SHAPES. 

235* 

c 


236* 

185 

00 2i: 1»1.N4 

237* 


X<1 l = n«l)*0EL 

239* 


T=R<1) 

239* 


PSl (1 )r.75*T*XL-T*T 

245* 


onisD* 

2*1* 


uA<n = c« 

242* 


DO 275 J=1*NM 

243* 


Z:XZ< J) 

244* 


HsT«Z«C33** XL/14? •*(Z*Z*T*n/<4«*XL)»(Z**44T**4>/t40«*XL**3) 1 

245* 


7PCT.CT.ZI 60 TO 195 

246* 


OsH-4Z*Z*T.#5*tT**3)/G.) 

247* 


GO TO 200 

24B* 

195 

G:H-(T*T*Z* •5*fZ**3)/6«) 

249* 

270 

U(1 ):Um*G*F<UI 

25r* 


1F( JOP7.eO«l> GO TO 205 

2!l* 


UA(n:UA<l)*G*fACj) 

252* 

2: 5 

CONTINUE 

253* 


IF (JOPT.EQ.n G9 TO 2C8 

254* 


yRiTE(6t2;>6i TiP$Hn«uni«uAin 

255* 

2C6 

F04MATf/,lXfF10*2t6E15.5) 

256* 


GO TO 217 

r«7. 

2‘.8 

kRl Tt C6*2 14^1 T*PSlU)*Um 

259* 

21? 

CONTINUE 

259> 

C 


26;^* 

C 

Y IS FOR SCALING PURPOSES* 

261* 

C 


262* 


NP2=2*NP 

263* 


no 219 lsl«4P 

264* 


Y41 IsPSn 11 

265t 


Y <1 .MPIiUlI ) 

265* 

215 

CONTINUE 

267 • 

C 


269* 

C 

HERE UE GENERATE THE PLOTS. 

2f 9* 

C 


27C* 


CO TC 4500.3DU*350«259.3?r.250li KOPT 

271* 

25 m 

CALL BGHPLT 

272* 


CALL PLFORM(*L INLlN^fXLEft.YLCN) 

273* 


CALL PLSCAL UfNPtNGtV tNP2«NGI 

274 • 


CALL PLA6ELi*THE FLEXIBLE BEAM CXFCR IHENT* « 28 t*lCHGTH4«6« 401 SPLACE 

275* 


1MENT*»12) 
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27i« 


CAU ^L«6Af 

27T* 


CALL PLAA|6(-2fXCCNf6« ) 

27S* 


CALL ^LCUAV(X»PS1«NP«IIT*7ICU 

2T^* 


CALL PLCURVtXfU*9P«N7«TlC2l 

200* 


CALL PLCUAVfX*UA«NP«NT«TlC31 

207* 


CALL PLTrXTU.5f.56««10f)»«49HOCS16£OI*l «S OPtlPALO) VS APPR0A15 

262* 


|ATt(/> $NAPCS«49,1) 

263* 


AO TO (271f?72*273*274«275)« n 

264* 

271 

CILL PLtCXT(3«S«7«0«*l*0**l3HTyO ACTUATOASi 13t 1 1 

265* 


60 TO 267 

266* 

272 

CALL PLTCXT<3«2*7.0«»i«0««15HTHACC ACTUATOAS f 1 50 1 

267* 


60 TO 269 

205* 

273 

CALL PLTrXTI3«3f7#0««l«0«*14HFOUA AC TUATOAS « 14^ > 

269* 


€0 TO 261 

250* 

274 

CALL PLTtXTO«3«7«a««l«0.fl4Hf IVe AC TUATOASt 14« 11 

291* 


60 TO 260 

292* 

275 

CALL rLTtXT (3«3«7«0«*lO*«13HSlX ACTUATORS «1 3«1 1 

293* 

26 0 

If(K0PT.tQ«4l GO TO 491 

294* 


CALL AOVPLT 

295* 

300 

CALL B6MPLT 

296* 


CALL PLFOA«<*HALIH»fXLtWf7LENl 

297* 


CALL PLSCAL(X»NPtMG«VfNP2«N6) 

296* 


CALL PLABCL(*THC FLtXiBLC BCAP CXPCRIHCNT ««26 •4teNGTH*«6« *2 1 SPLACC 

299* 


lPeNT*«121 

3C0* 


CALL PL6RAF 

3 Jl* 


CALL PLAX1S(*2*XLCN«0« 1 

3(2* 


CALL PLCURVIXf U«9PfNT«TlC21 

303* 


CALL PLCUAV<X*PSIfHPfNTtnCl) 

?:»4* 


CALL PLTCXT<2«5««53»«10«0.<31HOCSIR£0(*1 VS 0P7XMALI01 SHAPES, 31« 

305* 


1 1) 

306* 


GO TO (321«322f523f324f3251«R 

307* 

321 

call PLTEXT<3,3«7*0t*l*0«*13HTy0 ACTUATORSf 13U 1 

3C8* 


GO TO 33? 

309* 

322 

CALL PLTCXT<3,2«7«d«,3tO«»15HTHR£E ACTUATORS • 15 « 1 ) 

310* 


60 TO 330 

311* 

32 3 

CALL PLTEXT<3,3«7«0,,l«3#tl4NFOUR ACTUATOR S , 14, 1 1 

312* 


GO TO 531 

313* 

324 

CALL PLTtXT<3,3«7*4,,l,),,l4HFlVE AC TUA TORS, 14,1 1 

314* 


60 TO 339 

315* 

32 5 

CALL PLTEXT<3,3,7*0,«l,0»tl3HSlX ACTUAT0RS«13,1) 

316* 

33 ) 

1F<K0PT*CQ#?) CO TO 491 

317# 


CALL AOVPLT 

316* 

35 6 

CALL 86NPLT 

319* 


CALL PLFORP(*LlNLIN**XLtN,TLEN> 

320* 


CALL PLSCAL<X,NP,fHGtt,yp?,GGl 

321* 


CALL PLV8ELI*TMC FLEXIBLE BEAR E XPLR IPENT «,26 ««LENGTH*,6, *0 1 SOLACE 

322* 


1PEVT*,12) 

325* 


CALL PLGRAF 

324* 


CALL PLAXIS(*2«XLEy,n« 1 

325* 


CALL PLCURV<X,UA,NP,RTfTIC3l 

326* 


CALL PLCURV(X,PS1,WP,WT,TXC1) 

327* 


CALL PLTtXT<?*3,.56f.l J,3.,35N0tSIRE0l*) VS AP^ROXIRATE 1 /I SHADES, 

3:h* 


1 15 U) 

329* 


00 TO (371, 372f3T3*3T4, 5751,5 

330* 

371 

CALL PLUXT<3,3#7,0,,l,0,#13MTMri ACTUATORS, 13, 1 1 

331* 


CO TO 49? 

332* 

372 

CALL PLTEXTl3,2f7,C«,ltU,,15HTHREE ACTUATORStlS « 1 1 
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335 

336 
33? 
33f^ 
339 
3^ J 
36a 



GO TO 690 

ACTUATMS«14«1» 

373 

CALL PLTrXT<3.3«7«6..1«0.«X6HfOUX 
GO TO 690 

376 

CALL PLTtXT<S.3«7«0««lf 0««i6Hfl¥C 
GO TO 691 

ACTU«TMS«14«ll 

375 

CAUL PLTtXTn«3«7s0«*l«0*fl3HSIX 

«CTUATMStlS«l> 

69C 

CALL CNOPLT 


500 

STOP 

fHO 



1 

2 

3 

6 

5 
t 
? 

6 
6 

in 

n 

12 

13 

1% 

15 

U 

17 

IB 

19 

2 ) 

21 

22 

25 

29 

“ c 

« • 

2 S 
77 


SUBBOUTINC BMATU) 

COMMOM STABT«H$rAR«HMll«,HnAX«tf(NAX«i(t«Ke?fX2ll|i)«lif. 

6CAL AUOtU) 

DO 5C Jsltlir 

DO 

tXztZil} 

KJSX3< Jl 

CALI R0flBS<$?A9T«XLf K«FOfXfH$TAA«HMIMfNNAA«|Af^AA«ANS9KfKCn 

ID His X*XN <33 •« XL/160.41 xi«Xt#X*X)/M«^XU*< XI ««6«X**4|/I45««XL««3I) 

HJsX«XJ4<33.«XL/l60,4iXJ*XJ4X*X>/<6«*XL)*tX.J**64X«46)f <63«*XLt*3ll 
ir(X.€T,xn 60 TO 15 
ClrHt«(Xl4Xl4X/2#4<X**3)/6«) 
r.o TO 16 

15 C1shI-(X1«X*X/2*4<X1*43)/6.} 

16 ircx.GT.XJ) GO TO 23 
GJsMJ-Xj.XJ*X/2.*IX##3>/6. 

GO TO 21 

20 GJshj-X*X«XU/2«*<KJ**31/6. 

21 rOFXsGl^GJ 
CALI X0H2 

IFCK.CO.l) GO TO 1<> 

AU«1 )£ANS 
5 *» rOMTINUC 

DO 6>i lx2«NM 
T|=l-1 

no 6f> jsi«l I 
A(1«J)sA4Jtn 


26 

2V 

5': 


6V CONTINUr 
WLTURK 
fNO 


1 

2 
3 
6 

5 

7 

R 

9 

r 
!1 
!2 
;i 
:< 
IV 
1 1 1 
17. 
« o , 


SUBAOUTlNr BVCCCBl 

COHNON STAXT.HSTAR.HHlKtNHAXtCRHAX.XLtKeYf X2(1CI«MH 
RFAL nn o 
DO 5 > 1 = 1, N« 

z=<2<n 

CALL ROHRS (ST ART, XL fX,mFX ^H$ T AR ,hH I N ,HRA X , C R HA X ,A N S,K , KC Y) 

1*, 'i=A*Z*(53#*yL/l6T.A(Z.2^X*X,/<6.AXL)*lZ4464X446 7/<63.4XL443)> 

l/tX*GT«Z) GO TO 15 
r=G*IZ«24K««s.(x**3)/f«) 

'50 TO 20 

15 ':=G*<X.X.Z««f^ (Z*43I/6«1 

20 ^SI=.7S#X*XL*X4X 

f orxsrsMc 

CALL ROH? 

iriK^tc.n to TO 1 ' 

ii \ )=AUS 
5'’ CONTINur 

RfTURN 

nro 
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I 

Output 


[ Tht Li‘4>Tr4 OF THE OCAM IS I'JV.O*) 

( THE VECTOK af ACTUATOR POSITIONS 

j I TO 2 9.0800033*01 1.0300330*02 

' THE LITTLE A MATRIX 

COL 1 COL 2 

ROU 1 2.4299911*39 *4.1991314*39 

ROW 2 *4.1991314*09 7.14i>2920*39 

J 

THE LITTLE B VECTOR 

1 TO 2 9.4079597*00 *9.4246017*38 

THE BIS A MATRIX 

COL 1 

RCA 1 3.3462702*10 

THE BIS e VECTOR 

1 TO 1 *2*1999919*19 

FOR this weight ins VECTOR 0 

1 TO 2 0*1000030 0.0003103 

THE MATRIX 3IC A PLUS Q 

COL 1 

ROW 1 3.9462782*13 

THE FORCES F2 TO FM 

1 TO 1 *6*;961ia6*92 

THE VECTOR OF OPTIMAL FORCES 

1 TO 2 1.2192237*01 -6.8961166*12 

THE PHI vector 

1 TO 2 -l.l69':m*33 1.9992210*11 

THE APPROXIMATE LITTLE A MATRIX 

COL I COL 2 

ROU 1 2.4204669*14 *4.1487150*09 

ROU 2 -4.14(7151*19 7.1311936*09 

THE APPROXIMATE LITTLE B VECTOR 


1 TO 2 5.4616966*:* *9.3747, ’b9*38 

THE 6J6 APPRCX A MATRIX 
COL 1 

ROW 1 3*3447921*11 

THE »IG APPROX B VECTOR 
J TO 1 -2.033»:70*19 

THE APPROX MATRIX BIG A PLUS 0 

COL 1 

ROU I 3.34*7921*11 

THE «PPROX|MATC F'RCES R2 TO FM 
: TO i 
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THt APMOKINATC FOUCt VCCTOK f 


1 TO 2 

1.215 7442*^01 

•4.d767399«02 

^0$!rX0N 

OeSlRCO SHAPE 

OPTIHAL S«S*AP| 

«Q0 

.00009 

•00000 


.35009^93 

.1 047.03 

13.00 

.43939*93 

•9l472«93 

lO.O) 

.f0093^03 

.75154.03 

23.00 

.iioeo«04 

.94449.03 

25.00 

.12506«74 

.11510.94 

30. q: 

.13500^04 

.129T3.54 

35. n 

•140:0^04 

.13969.04 

40.00 

*14090#04 

.14453. 04 

45.00 

•13500.04 

.14314.04 

50.00 

.12S70«C4 

.13494.34 

55.:? 

.iiiC9*r4 

. 11935. C4 

49. o; 

*90^00.33 

• 1?fc9?3.93 

43.00 

«4SC0:«53 

.44355.93 

70.59 

.35o:o^:3 

.544bl.?3 

75.50 

.5933? 

•.39456.92 

60.09 

•.49003.93 

-.44229.03 

65. OC 

•.65732^:! 

-.91537.-‘3 

59.00 

-.I350:.r4 

•.i39:3*;4 

55.:; 


••16325.94 

103*00 

-.25COO^;^4 

-.23613.54 


I 


APPROX. SHAPe 

• 033i» 

•?99«3*)3 
•'j1224*03 
. 7*031* J2 
•94373*33 
.11*77*3* 
•1293A- 1* 
.139* :* 

•1**12*3* 
•1*274*3* 
.13*94*.'* 

• 11931**<* 
•94424*63 
•4P14**33 
•3*3«3*93 

••31575*:2 

-.*t:ir*33 
-.91274* J3 
-.13673.)* 
-.i«’7i*:* 
-.237*9*:* 
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B .3 The Sinply Supported Seen Eetimatioo Progres Liatlng 


1 » htAI. 

i* eiAL 

)• eiAL suiiiiu) 

«t* UATA *2St *ia7» / 

c 


/• 

11* 
ii'* 
15t 
190 
190 
^60 
!/• 
1 1 '* 
190 
iO 0 

^ 2 * 

25 * 

^40 

290 

d 90 

2 U 

240 

290 

il • 

%50 

i 4 m 

y 90 

yf 

yi 0 

y^0 

y90 

4i>0 

4\0 

420 

4^0 

44m 

440 

410 

4h0 

490 

900 

9X0 

920 

950 

900 

>90 

900 
91 0 
9 h0 
990 


1 004000 1NTIi0lUCT20I|«* 0444 
€ 

C THIS PfiOGHAH fCHfOAHS SHAH iStlHAtlCH fOA THC 
4 ; SIMPLY SUPPOMfCO dcAMe UPOH IHPUT OF HP JtSCM£T£ 

C ObSlPVAt ions blU Cf UlSPLACLPtAT AT POStTlOAS Aptllt 
C AN optimal cIaTImaTl Of THL SMAPi OF THC OLAM 
C UITm KiSPiCV TO THL LAIfiMlOA 

L 

C U/2U (SUM < » (I/O ill n^lNUL AAL 4 tU*F I«a^ II 

C 

C UhkMt LIAI IS THL SHAPi AND F4XI IS ThL NCISL IN Tnt 
L QYNAMU MOUtL CUCL*«AIU : F(XI« 

L 

C PLLASl uLFlNi THt f OLLOellNb YAMiAdLiS* 

C 

C XL IS The LlNi»TH 3r THi OtAMe 

C NM IS ThL NUPALN OF OKStM YAT IONS • 

L XZUlf Utfee.fNM AhL THC POSITIONS AL0N6 The ALAM At A^ICH 
C OtSLKVATlONS AML TAPLN, 

L 0(11 AML THc 00lk9ilX9t904 OF T Ht ^lUHTS ON THi OtoSiMYATIONS IN 
C Tht PlNFOMMANCC CMITLMlONeC 

C NP IS THl NiiMocN OF POINTS ON LACH CUMYt TO At PLOTfCb* 

L 

PUielAlSV 

XL = U 

NP = 3 

NUA-IC 

NUb^t 

Mbit 

NP;23 

CtL:XL/NP 

NP:.NP^ 1 

NPk:2*NP 

vO I lZim 9,4 

0( nr(XL«^n»IeL 7 
X2( :XzU90i 
X CONTlNUi 

bNirt(o«2YV) 

4.99 FOFPFTilhll 

yKlTUb«/l 

2 TOfMATl////! 

yMiT(,(o«mL 

9 KFPAT(1X«^SHThL LLNSTh Of ThL bLAP lStFU.21 

call VOiiT (X2tNM«26 «?AHjTHc OttSiMYATlON POSlTIONil 
CALL YOiil(t>*NP«II«I7H0THA OdSLMYAT I CMS I 
L 

L tOPFUlL ThL MATKIX A 
L 

uO S 1:I«NP 
tttZil I 

LIxX#(A-2**XL I 

uU 9 j:l«NM 
trx;<JI 

t2-Y*Y-2* ♦XL*Y 
Ci : X« V»V 


0\0 i: 4 zt 0 U 0 i 0 i 

420 A( I«ju i (X-XL I* (Y*ALI/( 30 e*XL«XLn« 4 4 X « « I 1/ 1 < X ••S I • ( 1 1 Wv e « 
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63* 


J <X»*3KCl*Ci/5.>««X»V/<36.»XL»XUI •((XL**7-V»*7)/7.4*L*Y»*6-)<L*»'« 

64* 


<2 « <7* (XL ••»-¥*•»)• a3.«XL*XL«C3)-(XL**4-V«*t» >*(3.ML«*3*XL«C3> 

65* 


3 ♦(Xl**3-V»«3»* <l./3.>» (5.*XL*XL»Cl**.*XL***t*C<l J - ( XL *XL -Y -Y !• ( XL * 

66* 


« C<>« (XL«*3)*C3>*(XL*YI*(XL*XL«C4n * (X/ ( 36. *X L**2 ) 1 *(Y -rL ) •( ( Y 

67* 


b **7-X*«7)/7.*)(1* ,b* <t**6-X**6l**2*(Y**5-X**5 1 * <C2*2 • *XL • XL* X* X > - 

66* 


6 «2b*( y«*4-X**4)* i 3.*XL»C2*XL*X«X) ♦ ( Y*«3- X «*31 «C2*U *X *2. *XL»XL 1 / 

6S* 


73.-( Y*Y“X*X»*,5* <XL*X*X*C2) J 

70* 

b 

CCNTlNUc 

71* 


DO 6 U2«N8 

72* 


11=1-1 

73* 


DO 6 J=> *11 

74* 


A(l«J)=A(JtIl 

7b* 

6 

COM IMUL 

76* 


CALL nOUT(*«KJA«Y1 •SM*13«t3H0T>(£ MATRIX A) 

77* 

C 


78* 

C 

HlAL JE compute A*d« 

7V* 

t 


60* 


UO 1=1 fNM 

6l* 


Ad(I)=^. 

82* 


UO 2U J = 1 »I4M 

63* 


Ab(I)=Ab(I>* AiltU)«6(J) 

64 <« 

2 J 

CONTINUE 

6b* 


CALL VOUT(AetNH*li»l!>HuT»ii ViCTOR A*in 

66* 

C 


67* 

C 

hLPt COmPUTl A^u* 

68* 

c 


69* 


DC 25 I=lt6M 

VO* 


uO <5 J=1.NM 

91* 


AC<ItJ)=A(l»j) 

92* 

2b 

CONTINUE 

93* 


DC 30 i=l*NM 

94 * 


AQ( ltI)=AC(lf I) * 0(1) 

95* 

3li 

tONTlNl’L 

96* 


CALL MOOT (AOtNUAfY M«NM« lb* IbHOTHE MATRIX A*U) 

97* 

C 


98* 

C 

hlhl y£ solve for thl optimal Shape at positions xz. 

99* 

c 


1 JO* 


CALL SGR< ACtNO AtVI t AB » NU8 • N6 * S 3 5 « 2 ORK ) 

lui* 


CALL V0UT(#»fc#NM*24f24HUOPTIPAL SHAPE PCSITICNSI 

iJc * 


bO TO 4J 

1U3* 

3 3 

UKITt(6*36) 

104 • 

36 

FORMAT (IHO. IX. 36H» ••••• MATRIX NEARLY SINGULAR ««••«•) 

lub* 


bo TO buO 

1 J6 * 

C 


1 J7* 

c 

NOW 2t COKPLTt The OPTIMAL SHAPE* 

1 

c 


lOV* 


LO Mb 1=1*KP 

UJ* 


Z( i )=( i-l )*wtu 

n )• 

4b 

COrv T IN J t 

U2* 


U U 5 J i = 1 f N M 

U 3* 


1)0 bb v-2*NP 

114* 


IM X2< i) .bf .Z( J>) bO TO bl 

1 1 b* 


x=x? ( I) 

1 A 6 « 


Y=Z( J) 

i J /• 


bC TO 32 

1 

bl 

X=Z( JJ 

11 V 


Y=X2<1) 



US# 

1^6* 

Ul# 

U9# 
UO# 
lii# 
13i»# 
lii# 
li## 
1 J># 
XSh^ 
liT# 
li«# 
li9« 
1#U# 

IHI # 

m;# 

I4i# 

1#H# 

US# 

lH6# 

U7# 

IH»# 

IMV# 

1>0# 

ISl# 

1S4»# 

IS3# 

IS## 

iss# 
1S6# 
IS7# 
is b# 
ls9# 
Ibvi# 
l6t« 

164 # 

I6i# 
164 # 
loS# 
l66# 
167 # 
1 6H# 
l6V# 
170 # 
1 M « 

1/4# 

I 73# 
1 M# 
I ;s# 
17 6« 
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S2 t2 = t#t-2.#)il#V 

C1 = X#( X«2«#XU 
CirX#X#V#t 
C#rX#*#y#Y 

$uuu ux-xi )# iv-xtui>6«#xi#)an # Mx#*n/7 4# <x #«s i« (Ci#c 2 usu 

1 ( X##31#U#C2/3.1 ♦ U#Y/ r^b4#XL#XUI • nXL##7-Y##7 )/7.#XL#r ••6-XL*#7 

2 ♦ ,2* IXI##S-Y##M#4 li«#Xt#Xl#CiU Ul##4-Y##4| I #« 3 , #XL # #3 ♦ xL #Ci ) 

3 ♦iXL««i-Y##3l#(l«/i.|#()4#XL#XUC3##«#XL####C4) -(XL^XL-'Y^Y^^tXl# 

4 C4# (Xl##il#C31#(Xl-YI#(XU#Xl#C4M # ( X/ ( 36* #Xl # #2 1 1 • ( Y - Xl 1 * ( i Y 

5 ••7-X##71/;.-XL#.S# <Y##6-X##6) ♦ *2 • ( Y *j- X ##i ) #< L2 ♦ 2 * •< L #XL • X • X > - 

6 .2S# 4Y###-X##4># U •♦XL#C2#XL#X#Xl#<Y##i-X##3l#U*IX#X-^2**xL#XtU 
7i*-( Y«Y-X#X)#«S*(XL#X#X#C2n 

SO CONUNUt 

u»mL< 6«S6> 

56 ^OHMA74///*lX«UhlHt HATf^lX SI 

Wt71TL(6»Sfl) (X2 ( U t UltMI^I 

57 |^^UK»»ATi/»lX»F lUU flutists I 
ftOH^AX (//«nXf lOcl S*5I 

00 Si 1-2»NP 

WRlTi(o*S7l 2aii(S4Jtil« 

U( ll^i* 

00 oO JiltKH 

U( 1 UJ U ) #(BC JUAs (jn#S(J#II/i( Jl 
00 CONTINUE 

WKlTECoU ) 
bK 1 1 1 (6 • 66 ) 
s>0 6S I'l 

PSH l)iXL#2 ( n-2 4 U *#2 

scu I ) n#i. 1 
sti i i#KP)=f SI n>#i. i 

WRlTL(6f67l 24 nt»SlU)#U( 17 

65 COMINUL 

6 7 F0R4AT(/tfli*2«4c:i 

66 FOHP AT (2X#oHP0S IT 10MtiX#|2HACICAL S hA Fc • iX 1 12Ht STl P # SFAPt ) 

XLLNcd* 

YILN-6* 

N0>1 

Tit*-*#* 

T (.• 

TUi = *#* 

NTU2 
NT4-0 
NT3=-l 
CALI bbNPLl 

CALL PLrO«N( •». IN.I W* fYLtN# YttNl 
CAll PlStAl 42ttaPtAb#SCl*taP2*Ai(>l 

CALL PLAbtL (*SHAPt lSTINATUN FOR Ht SIMPLY SJPPiKTci ilAM** 

1 46»*LtNbThOF Thl biAH* • 16#* U iSPLACLPt NT * 1 12 I 
CALL PLuRAf 

CALI PIAX1S(-2*XUK»U») 

LALL PLCORY<X2«A3»NRfNTitT ICII 
CALL FLlURV (2»U»NP*KT1*T IC2> 

LALL PlCUH¥(2*PSI#NP»MT2»T103> 

lALL PL TLXT (1 .4 I .S* *1 1 40»S3h 1 SUPAtc.b ShAPl <L I VS Co SL X VA T 1 CN • I 
ION AlTUAL SHAPLf >i *1 ) 
bO TO i74i»n»72#7if 74*7S1# NP 


1 774 

70 

CALL PLUX1 *CTU*TOR«U«i 1 

1 784 


CO TO »5 

1794 

71 

C*U PLTIXT (3*3«1^t.ltO.»l3HTkO ACTUAT0XS*13*n 

1 o3 4 


CO TO «0 

1814 

72 

CALL Ft TLXT(3.2*T.(it.l»0.*33HTT<Aei AC1UAT OfiSi lt> *1 1 

I8c« 


00 TO o3 

1 « 3* 

73 

CALL PLTLXli3*3f 7*ut*l»C*«l4HF0UR ACTUaT 0F5 tlNt 1 1 

184 4 


uO TO o3 

1854 

74 

CALL PEUXT i3.3t7*0»*l»U*tl4Hf IVC ACTli AT ONS »14« 11 

1864 


uO TO 83 

1874 

75 

CALL PETEXTt3«3i7*«5t*UO*tl3HSlX ACtUATOfiSt U«l 1 

1884 

63 

CALL EuOPCt 

18V4 

5uCi 

STOP 

1934 


ENO 


THE LtN&TH Of TiE bEAH IS 

THE OfcSEfiWATION POSITIONS 
I TO 3 2.3UJ0030-UI 

THE ObSERVATIONS 

I To 3 U 37:>O0oa-‘Ul 


THE HATRU A 


HOW 1 
ROW 2 
NOW 3 


thl matrix A*0 


ROW \ 
HOW 2 
ROW 3 


COL 1 
1 .‘ i6430lS - 0A 
1 * 4 *^ 02122-04 
I * 04t63l3 - J4 


LOL \ 

!• Jb330l>-04 
i . 4V02 122-04 
U 04S6313-04 


1«00 


3*000003 3-31 


2«So0003U-31 


COL 2 
I .4902I22-C4 

2 . lOt.l 33 1-34 

I •A9wkl2t>-C4 


I •Uttl}6621*U4 


COL 2 
1 * 4902122-34 
2 *1091301-04 
U49J212:>-34 


2*S432 742-31 


7*5000300-01 


I «4750333-01 


cOL 3 
l*045bT13-04 
1 *49021 25- 14 
1 *0622904-04 


7*6762742-05 


cOL 5 
I * 045aiii - J4 
1.49M2125 04 
I *0632984-34 


1*8421149-01 


0RT1^KL SMAft POSITIONS 
I TO 3 I* 8421271-01 


The Vlctoh a*o 

1 TO 3 7 *67627vi-iiD 
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Output 


THC MATRIX S 



•2b003^00 

•SjOOO^OJ 

•75030^33 


• ^3b60*U<l 

.3iV54-U4 

.2307I-C4 

«IO 

. <iob5V-0% 

•6blU9-U« 

.45591-34 


•6b3b5-C% 

.9bb6l-0b 

•670 1-C4 

• 20 

«ttb4in*0M 

.12306-01 

.66645-34 

• 2b 

• U623-03 

.14902-03 

.1C456-C3 

• 3U 

.12137-03 

•17351-03 

•11980-33 

• 3b 

•133MO-03 

.16761-03 

• 13213-4.5 


.14222-03 

•20340-03 

.14125-33 


.45 

.14744-03 

•2J821-03 

•14693-33 

• bo 

•14902-03 

•21061-03 

•149C2-C3 

• 5b 

.14693-Oi 

•20321-03 

• 14744-J 3 

• 60 

•14l2b-03 

.20046-03 

•14222-C3 

• 65 

•13213-03 

•16781-03 

•13346-03 

• 70 

•11980-03 

.17051-03 

.12137-03 

• 7b 

. Iu45b- 3 3 

•14902-03 

• 1062 3- j 3 

• 80 

•46645-04 

•1^366-03 

•66411-04 


•67021-04 

• 9>663-3« 

•68351-34 

• 90 

.45591-04 

•6bl09-04 

•46556-04 

.9b 

•23071-04 

•32959-04 

•23580-34 

1.00 

•26560-09 

•56935-09 

•10315-Lo 


PCSITION 

ACIUaL ShAPt 

tSTlM. SHAPi 


•3JjOO 

•U0U30 

.H5 

•97500-Ui 

•41587-01 

•la 

•9mUC0-31 

•81913-01 


• I;i!7504j0 

•11979400 

«20 

• lAJOO^ao 

•15416*00 

• 2!> 

•16750^QC 

•18421 400 

«3J 

•213d0^00 

•20929*00 

• 33 

•22750^C0 

•22902400 

• <I3 

• 2Au004>JU 

•24319*00 

• i»3 

*aA750^0C 

•25169*00 

• 53 

• 25000«>OU 

•25453*00 

• 55 

•29750400 

•25169*03 

• 63 

•^6 J0Q«3u 

•24323*03 

.65 

•k2750«00 

• 22 VC3 403 

• 7U 

•2100C^00 

•20932400 

.75 

• lo75J^00 

•18423*00 

• 6u 

•loOCC^vC 

•15416*00 

• 65 

•12750^00 

• 1 19 78*00 

• 90 

•9tooo*ca 

•81901-01 

.95 

•47503-31 

•41590-31 

1.00 

•7A506-U6 

•15969-U4 
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Appendix C. The Large Space Antenna Conputer Program and Output 
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OWGSriAl. PAGE IS 
gf POOR QUAtnY 


C«1 The Large Space Antenna Computer Program Listing 


3* 

7$ 

ICS 

M* 

13* 

ia» 

19* 

ea* 

?bt 

??♦ 

^9t 

30# 

3I^ 

32* 

$3* 

3af 

3S» 

3bS 

3Tt 

Jrt* 

39$ 

a\$ 

93S 

au« 

a*j# 

^s* 

aTs 

a(|f 

Sa» 

*yt* 


Dli^b^STON MDd(20)|FRtO(SO) 

0!**ib>SICM R(6R2) ,¥(81^2) f 2(9b2) iStir,5rF92) 

OtHpfjatO^ yeCT0R(<>b8b) f 0(892) fV (882) 

OlHtlMblPN 0CHK(ie)#P(l9)#KlM),A(l8,1A),8A(l«,lA) 

Ot^b>*9TUW PNtnt#|9) 

niMffigSUl^* BfcTAdO 

OlHg^^aiOb AtPHAdl), COfcPCll)# i^ORRMrift) 

V$T49(18) . AY(18)#uaT89(18) 
tfJTpr,tH iPT(i6)»jpTna) 
iKJTilitR JSC0(2428), I8tC(682) 

PQUlv^Lb^CE (lj(l),VKTnR(l))# (ytCT09(88i)#V(U) 
fOUTVALENCE (>^<l)#yfcCT0lld765)} 

C 

C 

C t8888888«f 

C 

c 

C static SH4Pt tStlHATlON ANH CH^TWOC PP A SPACf Ab-TEJ^A, 

C 

r 

r 

r 

C TmTS PRnttRANi fSTlMATtS AKjq C0*JT»<nL8 ThK STATtC OISTPRTrHN OF A 
C CaRGF space a^TENNjA# USI^g rest CCi»RDl*«ATtS» MfJOES A»-p FRHjuENCXfcS 

r SilPPtlFO HT A FIMTe t'ttPENT MOOFL, 

C 

r THF HUOFt l^CLNOES IP «IBS# TWfc Y, A^n I CnUPOl'^ATES FHW 882 
c POINTS, OR NnoES, LnrATEP tJ^- \a rO^StCOTIVF. CIRCLES^ 

r 

r IT IS ASSuHEO That the hub of TME AmTI>N.A is PIGIOLT ATTAChEO, 

r sn that these are sn rigid rodt nonrs, 

r THFPr ARE 33 FRtOoEMCIfS ANQ CORRESPOND I»^6 HCEnAUNCT IqnS 

r (HOOFS) FOR THIS HOOEL , 

r 

r This FR 06 RA*' READS FRJh thF TFMROKARV FiLk FIKFnANT, i*MlCH IS 

r CRFATtO F»(l^ THE TAPE A1960 BY RUNNING A PREUI'*INA»¥ PRf)GHAH 

c CRr,Tfn HY VEJAYARAGHaVah (VE.TAY) A(-'4H,(WEF. fO 3'i7«S0«A3R, 

r AUGUST IR, 1980, ftY VEJat ALhAR, 

c 

C The analysis oh which This PRPGHAR ys PASFO IS FOuNO Tn- 

r JPI FH 3ttT-H2, APRIL 0, 19A1, ’Ar-ApTIrjr. STATK SHAPE CO'^TSUL/ 

C PETF^^HjNATtUN ALG0RITHH8 F f)P Th(; U3F OF FpOFS SUPPLIFO HV A 

r F'lNlTt tLEME'XT OOfiFL* • PY CONUlt hFEAS, 

C 

r 

r Thf FULLONlHf, VAHlAHlES ^''I^T rE OEFT^KO,*, 

r 

t 

r Tuf vtCTOR AUPha(i), I»U^’^f are ThF CO^FMCILnIS OF THf finr#kS 

r THE actual OISTOhUD ShaFF, hhich js Tr. yf fsti'^atet anh 

r rn»t,fLTk o, 

r 

C 

r Thf OhSF^t WAT tons aNP tO'TWni PnSlTIOfA iFF dSSn’pO Tn LiU-'CAIkT. 

r ANn APPlItO TCI the same PIRLCTIOvS iT FaC^ ‘•PpT, 
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S7$ 

hjt 

•*s» 

^*t« 

67t 

f <>• 

To# 

71* 

/r* 

’i* 

7a* 

7S* 

7fc* 

77* 

7m* 

7g* 

Mr* 

Ml* 

M^* 

M.4* 

*4* 

M5* 

Mb* 

M7* 

*M* 

M9* 

^0* 

*1* 

*i?* 

gu 

94* 

9»i* 

9b* 

97* 

9k* 

99* 

19'^* 

101* 

1 

1 0 .^ * 

1 O4* 

n'i* 

< ''b* 
1 ‘‘7* 
\04* 
1 

\ i 0 * 

n 1 * 

IV* 
^ 1 5* 


t Thus I* * rnwtt IS tC t*ir v ntwi-CTlOK At 99, 

( tt !M AsSin^n YutWk IJI an 0Ha7*«WATlu»J Ur TKC AT MHUb 

r 99, COi^Vr^StUtf 
r 
c 

c t:i Tm* n* in n.iw Ap^anAiMaTlohi^i, 

c 

r Uf'l IS T»-L .%U’U>£k or AONCfS T'^ MA APPCIH)# 

c 

^ lPTfn,l»l»,«,,NPT IS tHf OP THf MjiiAL PtilNt 

r rn h,4tc»’ It*< PP«C7 IS A^«>uipr, 
r 

r iMUCATts Tnt rj9rCTTr.-^ nr Th^ tt-- P'^o, 

f J^T(il IS I#? 0*^ i, JPTfDil ‘^fANS Th^ pOKrt 1-9 TmC ( rj^^fCTTU’^ 

r 4»t7T)»<>,5 M*A\S Tnt POWCF TS I‘J Twp Y fUA£CtTfj^ nn / nt'^TCTllN, 

r 

^ MA‘Tl t,,i.,tn APPt* 4 Pi.lHCfeS# P F*r> t*' T^P T A'iP Y CTW^rilONS 

r aT - t4| U^t ThP I r AT . lr»b to, A**n «v,r 1 - Th^. Y 

r ulPpOTlOf AT bOOp. 39 , NPTbo, lPt«n*i, 19,10.39) A-iO J«»T«(l,^, . 

0 

0 

r *u *LP^p.^TS nA T;4P 9IA0NAL KHr,HTIvC; «-«ATPtt W A«C C*^nSf- HP 
C th* VAlUi PM, On l^PuT,,, 

t 99 t* tH iNttTAL VAl.Ut OP THA riAG'AL ALCOPNT* n* Tha AtUMTTNi, 

C «ATMt4 W. 

r I I THIS PWO^PaO H is ThA IWtPA.irMT ^Y »*»ich h*.* IsrPPAMAS, 

r tha rplUNi«.U‘ HY .%HirN g« is nAT£KHff,to ,,, tha shalcpsi 

r VAt.'jr Au*^ ThP ^'aTMIx a*p ts T*‘Vt PTTt^.Lh . 

c 

r 

r SAT U‘PT tHU*U To I IP NO PLOTS AMt nf.SlMLD, 


c 

r 



r 


r ThaSA VAPlAhlPS A9£ N,‘iT t'^PuT, TMf » elAtK rfeATJAU Ofc ‘-TLL 

r Hf. i»AA tKpU I J ThA MPOQMAm 

r 

r AxAfT ipSAM VATln^s OP tma: nisti>*to sha^f ,.|ll t»a co'^^riitpu 

c Af.n STOk'tn iv Tuf vrrT»J9 YSTah^ 

C VMT*« IS TA€ VPCTOP OA OPSA M V AT 1 n.«S AT THP Pot^TS TPTm, lAl.^A^t. 

r 

r OSTA5 V lUL LOnTaT:^ Tm£ nPTiMAl SHAPe ASTi:^ATtS AT Th£ P^TJ^TS 

^ tr'TfT), 

r 

r r*iAA IS 7 h£ v^ct'jm n* PSTjMATro cOPAPicipvts nr tup «onAS 
r 1‘ T^^ t ISTUMTtl S^^APp, 

r 

C p ts ^tCTf'M (;A rrTjMAl C"l’'‘TMfu Pl.wrAS, 

f rfTl IS APPtItl' TO jPTtn T* TmA OlMtCtlON .TPT(l), 
r 

c ^iTt IS Tht VfCTL.H ir O.inAL c ^KAA If TA »:Ts -»P Sl'L T J *-.r, PMCH 

r Tr-f AA'Pi iCATins CF Tl*A .jOTr*At CO»TPO| AOtfAS, *>Xli MA 

r (f tt.'TA.J HASto MV; TwA hPTtMA SWAA^t P5TI‘;ATA STHApp f>j Tuf VACTn<< 

1 Cf.AA, 
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iu> 

US* 

119* 

1?4« 

!«!>• 

I?b9 

l?H9 

I29f 

I50t 

1 S2« 
lillt 
19«9 
tS5# 

UTt 
13«» 
1 39i 
1909 

\A\$ 

\a?9 
« 41$ 

1 aa$ 
\4b9 

|9fi« 

147$ 

14»t 

149$ 

\S0b 

1S1$ 

151$ 

1^«$ 

1S5* 

\^6$ 

l^T$ 

1SS$ 

lb9$ 

1^0$ 

1*11$ 

!•>«?$ 

l*^1# 

is$$ 

1^S$ 

t**^9 

197$ 

199$ 

199 $ 

170$ 


c 

C IHO t» TMt nmhbimino Of THE Nonis, 

C NC 15 thi K}JHUb Of fLOT fOMHA^OS, 

C J5f0 !5 tht SfOUENCt OF PtOt COMNANOft* 

c 

c 

C Th9 L^5T circle H$5 OEtETCD fbO^ Th| PtOTS 

r AfCAUSf Of 590 OETi POl^iTS ON THAT cI»Clf U 80>^t 

c 

C If IT 78 r»t8!W|0 that ThF IA8T CT«Clf Hf PLOTTEO# 8PT NCt|9^8 

C IN’ TMi 09TA STATfMENT, |W TH| 8UHPmJTl^l 0*<A9, A^O 

C TMt LOaP IK^OlVliyO C09Tl9Uk. 

c 

r mp 9 rs tme headinq qh the ailc# 

C >*ODtb 15 TM£ 9UP8EW Of Nnof$«888, 

C ^F9fC 75 TH| NUP5CP Hf f H£OUt^'C If 3«50* 

C 

c 

C $»$$$«»$9«$$$«9M««9««$$9»$»$$$9$$9«tt»«$»9$t9»»t«9$9$$««»t99«$$t 

c 

c 

r DATA CA905,»« 

c 

OAtA AlfHA/lC,.iPO,tlO,,5$0,,^,,r,.«».»r^/ 

PATa ^'C#^*^i’^PT/2l90,U#l5/ 

C 

c FI$Tm CI*^ClC 

r 

OATA IPT/187* I30f I ij»l5d» 139«1U?» 14^,146 ,1^1 « 19ft, 197 •l^^f 19 If 196| 
1 I99,l7?n7b#l7ft/ 

OATA JPT/ia«1/ 

ISYOlu 

lOPTii. 

P9tUUt«iO 

»»« ,*n» 

c 

C wcin ►M'T CO«H*MO 8tOUtNr» 

c 

»Mor8,?«) uscw(n,t«i*'in 

c 

c htflj *.f otttTt tn£ ci*'cir thc •tot cor ■•*'' 08 , 
c 

nn ;>s 

J$aI*t'S(JSt.b(KK)) 
t* f js.tt.trs) to TO f<3 

<8 Cn*.-TtMJt 

c 

C otrf.f ♦TiOf akOuf»<ct isroci),,, 
c 

( r> ?• 

l9forJ)a> 

fl’.''Tj(-U» 

C 

irnnHT.ki.u go th •<« 
t*tt •tuts 
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Itl# 


CALL f>U)T(k.0««,0i<i) 

nif 


CALL MC^0f)(.00*a) 


450 

CQNTtNUC 

1?4# 

C 


1T5» 

c 




00 4 

iTt# 


cocr (IN)BO, 



«eTA(l*<)*0, 

IT9* 

4 

COMTIHUt 

1404 


00 44 jn«i»npt 

t0l4 


DCH)((JN)bo, 

l4*f 


V$TA4(J»<)aO« 

14J* 


Ay(JK')ao. 

1944 


USTA4(JN>B'\ 

1494 

44 

CONTtNUl 

1444 

C 


1474 


as 

1444 


HEAD(«S) (H0R(K),K41,<0) 

1M4 


<*eAO(«9)NOOES,NpReo 

1404 


NCMcrNaNOocats 

1414 

C 


14J4 

c 

Atnivau^ci) are the cooi>oihates of node i. 

1454 

c 


1444 


«EADraS)(xm,T(n«z(n«T»t.HooE8) 

1444 

c 


1444 

c 

SLUGSd) la THE Misa at hooe t. 

1474 

c 


1444 


HEAor«S} (aLU6a(i),xat«Htit>Ea) 

1444 


•»rA0(A5) (F«E0(I)(I"1*HFRE0) 

2004 


write ( 4f to) (HDR(K),KRt*20) 

?0S4 


to FORMAT (t Hi,//, tax, aOAA.SX, *FRtOUENClE8‘i//) 

2024 


wRITECa* t2) (rRCO(K),KRt,HFREQ) 

2014 


12 FORNiT (5X,aEt9.a«/) 

2044 


WRITFU,U)HnOE8 

2054 


14 FnRH*T(//,40X, «N0. OF H00E8 a •,!«) 

2044 


w»!TF (AilS) 

2074 


WRITE (4»ia) 

2064 

15 

roRNtT(//,<X,<Foa|TlON8 AND OlRCCTlowa OF CONTROL/OaaCRVATlOH FOIH 

2044 


1T8') 

2104 

14 

FORH*T(//,«X,'NOnt',FX»»OXRECT10M) 

2114 


00 ?0 iBlfNFT 

2124 


JaajRT (1) 

2114 


XF(ja»2) iT>ia»i4 

2144 

IT 

«RiTE(a,«a) xFTtn 

2154 


60 TO |f> 

2144 

14 

w»TTFla,«S) !FT(X) 

2174 


60 TO FO 

2144 

14 

WRTTf(fc,«a) XFT(X) 

2144 

20 

cnwTtHUi 

2204 

22 

FORHAT (/«2X, I4»4X, *«•) 

2214 

29 

FpRH*T(/,«X,Xa,4X,'t*) 

2224 

24 

FORMAT!/, ax,ia*4X,*2*) 

2214 

C 


2244 


UO ino aFai,N»i 

2254 

c 


2244 

r 


2274 

c 

WTRM9 ta a CHECK T'1 8lF ThbT THE TAPE I8 ttlXAiC REAP PRCPERLT, 
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C 

MTVNt ■ 


219* 

C 

TS tHt Hoor (IIQCnvCCTOR) NUMBCR* 


2^09 

C 

FA IS TH| FACBUCMCV. 


111* 

C 



2129 

C 



2S39 


Ae«Q(«S)AFA(FA«NTANS»(VfCrOA(K),K«l,iiir*NS) 


2999 


XFr*«CHeCR*N(.NTKiNS)60 TO 1S9 


2199 

C 



2999 

c 

A^T(T»J)* and JB|,NAr HOtOI THf VAkUE OF 

MOOl X 

29T9 

c 

AT NOOC IFT(J) XN The OXAECTXON JFT(J)« 


2199 

c 



2999 


on 3S XitrNFT 


2909 


J«TPT(X> 


29i9 


JSajFT(X) 


2929 


IF (jSof) lOiStfSE 


2939 

so 

FMT(«F,X)aU(4) 


2949 


en TO IS 


2999 

31 

FMH«r,X)aV(J) 


2999 


on TO IS 


29T9 

32 

FHKAT .X)aa(J3 


2994 

33 

CONTIMUt 


2999 

35 

CONTINUE 


2909 


00 SO laliNOOCS 


2519 

C 



29?9 

C 

M£«r Mt COPFUTt The aNOAN OXSTOATCO SMAFE, 


2939 

c 



2999 


X(t)aA(I)«AUFHA(KF)au(t) 


2999 


» (I)a»(lJ»AtFHA(AF)*V(l) 


2599 


Z(I)«Z(I)*ALFHA(KF)*A(X) 


2579 

50 

CONTINUI 


2599 

ICO 

CONTINUE 


2599 


IF(tOFT,E0,n SO TO lOS 


1909 

c 



2919 

c 

here we r»ot THE anonn distorted snare. 


2924 

c 



2934 

c 



2944 

c 

ORAN IS A subroutine created ST 6, ROoRISUCZ TO 

•tOT three 

2954 

c 

OINCNSIONAL surfaces. IT CAUS TnE susfoutinc 

TRANS, 

2994 

c 



2974 


CAU DRAa(X>T,X.4SEO«ISE0) 


2994 


CALL MCTORU.O) 


2994 


CALL FLOT(10.0*0,0»«3) 


2704 


CALL FACTUR(,00*S) 


2714 

1U5 

CONTINUE 


2724 


CALL noiiT(R«I,na,n«,nrT.U,15NOTH| haTRIX RHI) 


2734 

C 



2744 

C 



2^94 

c 

HERE r-t COMPUTE The Matrix a and tnf vector of exact ossfavations 

2794 

c 

VSTAR, 


2774 

c 



2794 

c 



2744 


00 IRO IRJ»NRT 


2904 


p(t:RU, 


2914 


00 IRO JR1«NFT 


2924 


A(I,j;ao. 


2934 

l9o 

CONTINUE 


2944 


00 lAai.NH 



ISl 


^9J• 

«n* 

«9T* 

<9«* 

aM* 

soo* 

] 0 |* 

so<« 

30S» 

30«« 

309* 

309* 

30T* 

309* 

309* 

3t0« 

311* 

ita* 

111* 

51«* 

115 * 

319* 

31?« 

319* 

119* 

lao* 

sai* 

laa* 

las* 

la** 

jas* 

399* 

lar* 

329* 

329* 

310* 

331* 

332* 

333* 

33** 

343* 

339* 

333* 

33** 

33** 

3*0* 

(«1« 


on ano i^unpt 
cocrtiiao, 

99r«9(navsr*M1399L9HA(iiO«pt4inK,n 
00 200 3«l»**2T 

aoo co9T)9Uf 

C*U •'nuT(*,N2T,*»2T,*i2T,|S,j3«0TMf 9»T9|«l «) 

C4U VOU»t*ir*9,99T,io,|OHOTM| 9|CTP« 02 0M£»V9T10M »*) 

C 

c cpmoutatiom OF THi FROoucT «(yar*K)a«y , 
c 

00 202 l•l•NFT 
or 202 J91»hFT 

9V(n***(I)**a»J) * TITaRU) 

aoa CONTINUE 

e*Lt *OUT(*»* NOT# 19,t9N0THE VeCtON *») 

c 

f. H£Rf «c 900 TFC 0 HaTRIII TU THC 9 NaTRIX* 

c 

c 

00 203 l«l»NFT 
9n,t)*9(I«I)*9R 
203 CONTTKUe 

209 F09NAT(//,2X,«Ha *,Mi,9) 

203 CALL N0o3(9,NFT,NFT,wFT,l5,t5M01 'E MAT»IX 9*») 

N9lTf(9»a09) 99 
00 2'.l> X*1»NFT 
U«TaR(I)*9*(1) 

0(1 21U J«l»MFT 
9A(I,J)«Aa,J) 

21 0 CONTINUE 
C 

C 

C MT9F *1 NOFf TO SOITC THE •V9TfH (H*A)U*aAT» , 

C 909 18 A JFL klNE99 EQUATION SOLUTION HOliTlNe. 

C 

C 

CALL *09(94, NFT,NPT,0|Ta9 , NFT , 1 , 9290, *099 ) 

CALL V0uT(U*TA9,NFT,a*»29NOVtcTn9 OF OFTIHal ESTIHaTES) 

00 215 I»1,NFT 
00 215 4«t*NFT 

OCHK fn«DCNN(l)*9(l,4)9U|TA9(JJ 
215 continue 

CALL V0uT(0CH*,9FT,19,19MOThE VfCTOfc (A*N)Uc) 

00 220 l«l,N*t 
FQNarENE0(I)*9«)999 
90 220 J»l,MFT 

Cr)(F(I)*CUeE(I)«(T9TA9(4)*UaTA9(jnBFHl(l,J)/FUN 
220 CONTINUE 

CALL tiI|.T(ALFN4,nn, 59, 4*N0THt TICtpN Of ACTUAL C«‘EFF IC ICNTS) 
CALL vn|iT(C0Ef,NN, ji.jImoTme VfCTO* Of t8TINATfO COEFFICIENTS) 
f 

C NON wf CONFUTE TN(t tSTlNATCn 9NAPE, 

r 

r 


NEnImL *5 

NEAOfNb) (N09(h),«al,20) 


152 




*C«P(«f)*iOOU»NPHO 

S4St 


NeMec*^*NOoft*s 



n*oMt)U(i)»T(n,tu)«i«t**ioort} 



n*o(«t) (Kuai(n,ui,Mooct) 



•C«0(«S) CfM0(l)i iat»NMrO) 

30# 


on tao uraiiMM 

JM» 



iMf 


ir(NCHfcK«Ne,««TRM«)oo TO Its 

330t 


00 |«0 l■l*NOOCI 

J3|» 


icn«i(x)*cQfPUP)*u(i) 

3S2t 


T(n»»ci)*coiPc«F)*v(i) 

393* 


I(t)«2(l)«C0CP(ilP)aM(|) 

33#^ 

140 

CQMttNUf 

3Mf 


tP (tOPT«CQ*l} 00 to 900 

3349 

C 


33T9 

c 

MCRC fC Ptor TMf COrtOAtfO OHiPC. 

3309 

c 


3399 


c*|.t DR«M(«»T»tfJ0CO*t0eO) 

3409 


Ctu PACTOR(UO) 

3919 


CALL PL0t(10„e.0««S) 

34{9 


C*U P«CtOH(,OOftO) 

1439 

c 


3449 

c 


1439 


09 TO 900 

3449 

*50 

H0ITC («»t91) 

34?9 


t9t0al010*l 

3449 

13t 

P0RMtT(t9R*K0N0a«rRII> 10 N0ANLV tINOUUR) 

1499 


MlTt(*»090) 

3t09 

t3l 

Pnati»T(//»i9ii»9N0HlT,) 

JM9 


irtllXO.OT.O) 00 TO aoo 

3T49 


Hain(0|i9S) 

|T19 

<33 

aONN«T(//«19«ft9HREDCPl<«E Th*. HaTRII A«R< 

3T«9 


DO oao 1«1«N0T 

3T39 


*(T,naA(l«I) ♦ M 

3?49 

<40 

CONTlhUe 

IT19 


RRalO.tPR 

3T49 


•aO,*RR 

3T99 


00 TP lOS 

3409 

C 


3419 

300 

CONTINUE 

3449 

C 

MCap NE COMPUTE ThE VECTOR 0 IN THf CONTROL PRO0LEM« 

34J9 


00 HO i«i»nn 

3449 


00 JtO 4*I»NPT 

3439 


0(J)at)(J)4PH!(I,J)*c0EP(n/(PRe0(naa{) 

3449 

310 

CONTTNUl 

3479 


CALL V0gT(0,NPT,t9*l*N0TNE VECTOR 0) 

3449 

313 

CALL P0UT(*#NaT,NRT,NPT,19«l9HOTHE MATRH ATR) 

3499 


NRTTP (A»t04) RR 

J909 


on JPO |a|»NPT 

3914 


P(t)ab(t) 

3919 


0CNH(|)80, 

3919 


PO S«c JNI.NPT 

3949 


AA(l,J)a*U«J) 

3919 

KO 

CONTlNUt 

1944 

C 


3974 

c 


3944 

c 

HERE Ni hOPC (PEHVEntLV} TO fOLVt TmE naTBIk (A*a)PaO, 
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\ 

I 


i 

i 


f 


1 

f 


i 


) 


! 

V 


nn 

00 

01 

01 

ftj 

0« 

0^ 

0 % 

or 

01 

00 

to 

tl 

12 

13 

U 

IS 

U 

IT 

10 

10 

20 

21 

22 

21 

20 

2S 

20 

27 

20 

20 

30 

31 

32 
31 
U 
35 
10 
37 
30 
30 
• 0 
•1 
«2 
«3 
00 
05 
00 
0? 
ao 

ttO 

50 

St 

5? 

5S 

So 

55 


C ono IS 0 4PL LlNOAO (QUOTION OQtuTlON OOgTXt^lt 
C 

c 

CALL S0ll(00|l«2T»O0T,rfNfT»|»0}Sp«00OlO 

cut 70lT(0,NOr^|9,{fHOVIC?0« 02 02TIMt FOOCIt) 

00 325 \0|#N2T 
on 125 J0|fN2T 
ocMo(i)oocMi((t)oo(x«j)or(j) 

S2S CONTTNUI 

CAU 7Qtt(0CH0«N2T«lOfUM0THC VCCTOr ;OoO)f) 

00 3on fifHH 
2Qa20AU(t)o0| 
no 310 J«U>*2T 

oiTAm»ocTA(nof cJ)02Hic!f j)/ra 
330 CONttHUf 

CAu voutcicTAi MMf 30f soHOMOOot cri22icicot| ronM contool oomcco 
i) 

CALL V(tor(AL2MAfl^N|30,34HOTHE vrCTOtt 02 ACTUAt Cn€22IClt<«ti) 

c 

c 

C NQ- ¥t C0M2U70 THE $HA2t AOiUOT*«tNtg 

c 

f 

00 SIS 1*1#N« 
co72n)oAL2HA(n«eiTA(n 
m coNTtNut 
OColKO OS 

OrAD(OS) (HON(lO|R«i«|0) 

0|ACCoS)NOOESfNflleO 

^eMfcR45upcoos 

oc4f)(o5) (*(n#t(ni2(i)ii»i#5uofs) 

ME AO f 05) ( 0LU6S ( n t lot I NODM) 

0EAD(05) (2000(1)1 l«liN20r.Q) 

Pn OAO a2o),nh 

0( on (05) 020,20# NitONO# (7fCTOO(tO#r.ot,WTONi) 

lF(f CKECR.NC.KTONS) 00 TO 125 

00 3at X01#NOOC8 

*(noACl)^CJE2(R2)ou(X) 

vC?)tMnoCOA2(RA)ov(X) 

xcno/rn^coEr (o2)oou) 

Sou CDNTt^U^ 

12 , '<OT,EQ^U GO To 5^0 
C 

c HfO» OLUt •«£ cO» 2 tCTrO SHOOt. 

c 

Call tOA«»(v#T,2,4$ic,Xlcu) 

Call 2actuo(uo) 

CAiL 2LU^(in,n#0,0:»S) 
rUL 2ACTOO(,0006) 
on Tn ooo 
A50 •‘OlTf (0,251) 
istroiixw^i 

fr (!8lG,6t^}0) GO TO 800 
•0tTC(0,253) 
f'n ur x»lr^2T 
An,t)lA(l#l)oO 

iifcO 



00 TO 3iS 


154 i 


COMTTNUC 

eoMTTNue 

CONTINUE 

iraoPT.10,1) 00 TO 810 
CEU PLOT(1O,#O,»000) 
CQNTthUe 
STOP 


0964 

0874 

0994 


0994 

c 

0609 

c 

0614 

U9 

0624 

«oe 

0634 

0604 

0684 

soo 

0664 

0674 

0604 

510 


<0 OP comiutigni no ozAGNoarica* 

CTPUO07 auPai9,30O 


!• 

2» 

3P 

4P 

8« 

04 

74 

04 

44 

104 

114 

124 

134 

104 

194 

164 

174 

164 

144 

204 

214 

224 

234 

204 


C 

aUBNOUTZNE 0PAH(UX»U7,UZ|a£Q»aEQl) 
c 

C ORAN IS A SUaROoriNE CPCATEO BY 6 , ROORIGUEI TO PUOT THREE 
C OINEnSIONAL auPPACES* IT CALLS Th£ auRROUTlNE TRANS, 

c 

PARAHETCR NP 4 S 02 # NC 4210 O 
REAL UI(NP),uV(NP),u 2 (NP) 
integer 3 E 0 (NC)fPLAGfSEQl(NP) 

00 10 XiilfNC 

J«|APS(SE 0 (I )3 
00 9 K 41 #NP 
IP(SCU 1 (R),EQ«J) KKRK 
8 CONTINUE 

PIAG «2 

IP(Sf 0 (Z 3 ,LT,O) PLA 0 S 3 

% m 

VaUV(AR) 

ZtuifKR) 

CALL TRaNS(K|Y#Z,XP^VP) 

CALL PLPT(AP,YP,FLAG) 

10 CONTINUE 
RETURN 
END 


14 

24 

34 

04 

94 

64 

74 

04 


SURRHUTINE TRANS{X,Y,Z#XP,YP) 

real X#Y^ZrXP»YP 

THET4B30.0 

ORa 3 * 10 l 6 /lSO,O 

XPb(X#Y)4COS(TmeT04C)A) 

4P4(I*Y)4SIN(TMETa40R)42 
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